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Propagation of TE,, Waves in Curved Wave Guides 
By W. J. ALBERSHEIM 


TEo. waves transmitted through curve wave guides lose power by conversion 
to other modes, especially to TMn. 
This power transfer to coupled modes is explained by the theory of coupled 
transmission lines. It is shown that the power interchange between coupled lines 
and their propagation constants can be derived from a single coupling dis 
criminant. 
Earlier calculations of TE, conversion loss in circular wave guide bends are 
confirmed and extended to S-shaped bends. 
Tolerance limits for random deflections from an average straight course are 
given. 
HE TE : mode of propagation in circular wave guides has great 
potential value for the transmission of wide-band signals because 
its attenuation decreases with frequency. In order to take full advantage 
of this property one must use sufficiently large wave guides to operate well 
above the cutoff of the lowest transmitted frequency. The difficulty of this 
transmission method lies in the fact that TE: is not the dominant mode 
and that energy may be lost by transfer to the many other modes capable 
of transmission in the wave guide. In an ideal wave guide, which is perfectly 
straight, perfectly circular and perfectly conducting, the propagation is 
undisturbed; but slight imperfections and especially a slight curvature of 
the wave guide axis may produce serious disturbances. 

The character of these disturbances has been investigated in several 
publications by Prof. M. Jouguet' and in unpublished work by Mr. S. O. 
Rice of the Bell Telephone Laboratories. Both Jouguet and Rice use the 
method of perturbations, which is a form of calculus invented by astronomers 
to compute the deviations from the exact elliptical orbits of the planets 
which are caused by the disturbing influences of their fellow planets. 
Although the above-mentioned authors obtained valuable results, the 
interpretation of their solutions is difficult due to this rather abstract 
mathematical formulation. To most engineers the understanding of a 
physical problem is greatly helped if it is possible to use a method of analysis 
which is elementary in character and easily interpreted in familiar physical 
terms. The familiar concept on which the present treatment will be based 
is that of coupled circuits. 

1See References:2 and 3, listed on page 7. 
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It has been stated by the earlier authors that the curvature of the wave 
guide produces a coupling between modes. Before going into a detailed 
analysis one may estimate by inspection the nature of this coupling and the 
kind of modes that are most strongly coupled to each other. Figure 1a 
shows the cross section and the longitudinal section of a straight cylindrical 
wave guide. The location of every point inside the wave guide is determined 
by three coordinates: the radial distance r from the cylinder axis; the 
azimuth angle ¢ from an arbitrary 0 line and the axial distance z from the 

















Q- CYLINDRICAL CO-ORDINATES IN STRAIGHT WAVEGUIDE 











b- TOROIDAL CO-ORDINATES IN CURVED WAVEGUIDE 
Fig. 1 


origin. If the wave guide is bent as shown on Fig. 1b, but a wave front 
at right angles to the cylinder axis is to be maintained, the waves must be 
shortened at the inside of the bend and lengthened at the outside of the 
bend. Regarding compression as a positive and expansion as a negative 
deformation, one sees that the distortion of the wave shape is proportional 
to the curvature of the wave guide multiplied by the cosine of the azimuth 
angle. It is natural to assume that the coupling between modes is propor- 
tional to this distortion. 

Now it is known that all modes of propagation in a circular wave guide 
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can be derived from functions J,(xr) cos ng. iIn these functions, is 
called the aximuthal index because it indicates the'type of symmetry around 
the circumference of the wave guide. When theée characteristic functions 
are multiplied by the distortion factor cosine ¢,' the resulting expressions 
are proportional to the sum of cosine (x + 1) ¢ and cosine (n — 1)¢. This 
means that the bending of the wave guide couples mainly those modes which 
differ by +1 in azimuth index. Since the TEs: mode has the azimuthal 
index Q, it is coupled to all modes of the type TE;» and TMi . 

In the above qualitative discussion we have claimed that coupling exists 
without defining the physical coupling parameters and their effects. We 
must now supply this definition and show that the TE»; mode is particularly 
susceptible to coupling losses. 


neice Ay —1; 
+ ” M3 % iy 
Z4 | Zy | 21 | 
z y, Zz y, Z 


m m m m s-m Ym 
Z 








1 











y 4 & 
ee ae a os 
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Fig. 2 


Our investigation is guided by S. A. Schelkunoff’s statement’ that a wave 
guide mode has the same equation of propagation as a high-pass transmission 
line. Schelkunoff further points out* that the high-pass character of circular 
wave guide modes can be interpreted as the effect of interfering plane waves 
whose directions of propagation deviate from the wave guide axis by a 
constant slanting angle. 

We therefore approach the problem of coupled wave guide modes by 
studying the behavior of two coupled transmission lines such as shown on 
Fig. 2a. Each transmission line is schematically shown as an array of small 
ladder sections. The series impedances per unit length of the lines are 2; 
and 2); their shunt admittances per unit length, y; and y.. The two 
lines are loosely coupled by small mutual series impedances per unit length 
(g,) and by small mutual shunt admittances per unit length (ym). 

A network of coupled ladder sections is more tractable than a wave guide 
structure, but still somewhat complicated. Let us therefore carry the 
analogy one step further. Figure 2b shows two resonant circuits, each 


2 Ref. 4, pp. 378 and 381 of the book. 
3 Ref. 4, p. 410 of the book. 
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consisting of single capacity C, and an impedance Z which includes an 
inductance and a damping resistance R. The resonators are coupled 
by a small mutual! inductance Z,, and by a small mutual capacity Y,, . 

The behavior of coupled resonators is very well known to radio engineers. 
They occur as tuned transformers in amplifier circuits, as band-pass filters 
and as “tank circuits” in radio transmitters. Even before the advent of 
radio, their acoustical equivalents were studied in the form of resonant 
tuning forks. The mathematical aspects of this problem were already 
clearly set forth in a paper by Wien written in 18974. He showed that the 
interaction between the free vibrations of two tuned circuits depends on 
the coupling coefficient and on the ratio of their complex resonance fre- 
quencies. The closer the two frequencies are to each other, the less coupling 
is needed to transfer energy between the two circuits. The reason is that 
the individual free vibrations of two nearly synchronous circuits remain in 
step long enough to accumulate the small energy transfer impulses of many 
vibrations. 

Now consider the two transmission lines of Fig. 2a and assume that a 
constant frequency signal is impressed upon the input of one or both of them. 
The signals are carried along the two lines as traveling waves. Again it is 
true that loosely coupled signals affect each other strongly if they remain in 
step. With traveling waves “remaining in step” means that they must 
travel with approximately equal phase velocities. We conclude that the 
phase velocities or phase constants of coupled transmission lines play a 
similar role as the resonant frequencies of coupled tuned circuits. This 
intuitive reasoning is confirmed by analysis (see Section 1 of the analytical 
part of this paper). 

We thus find that we must expect trouble for TEo: wave guide trans- 
mission if a mode with an azimuth index 1 has a propagation constant 
close to that of the TEo:. It so happens that there exists one mode, the 
TMi, which in an ideal wave guide has exactly the same propagation 
constant as the TEo:. This then should be the principal source of trouble— 
and from previous work it is known that such is the case. 

Our discussion of coupled transmission lines has shown that the interaction 
effects are functions of their relative uncoupled propagation constants and 
of the coupling coefficient. The propagation constants of the TE9, and 
TM, wave guide modes are known but their coupling coefficient remains 
to be found. 

Since the energy of the transmission modes is located in the dielectric 
inside the wave guide, we consider first the coupling between the plane 
“slant wave” groups from which the modes are built up. 


4 Reference 5. 
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As shown in the analytical part, the coupling coefficient of these slant 
wives may be defined as the energy interchanged between the modes per 
unit length of line divided by the geometric mean of the energies per unit 
length stored up in each of the modes. 

From the coupling coefficient of the slant waves the coupling coefficient 
of the wave guide modes is derived. 

On the basis of the above physical interpretation the analysis is carried 
out and the properties of TEo; propagation through curved wave guides of 
various shapes are derived in the analytical part of this paper which is 
subdivided into the following nine sections: 

Section 1 develops an approximate theory of loosely coupled, weakly 
damped circuits. The theory is first derived for coupled resonators which 
are familiar to communication engineers, and then applied in similar form 
to coupled transmission lines. It is shown that the important interaction 
properties of coupled lines are functions of a single coupling discriminant. 
The relative energy content of the two lines in each of the two possible 
coupled modes is plotted as a function of the coupling discriminant. 

Section 2 contains the field equations of a straight circular wave guide 
and their modification by a toroidal bend. 

Section 3 gives the solutions of the field equations for the uncoupled 
TEo: and TMi, modes in wave guides with infinite, and with small but 
finite conductivity. 

Section + applies the coupling theory to the TE ; and TM; modes in 
circular wave guide bends. The coupling coefficient, coupling discriminant 
and energy division between the two modes are derived as functions of the 
wave guide diameter bending radius and conductivity and of the signal 
frequency. 

Section 5 derives the critical bending radius and the attenuation of TEo: 
waves in long wave guides of constant curvature. Two numerical examples 
are given. 

Section 6 shows that in a curved section of wave guide which follows a 
long straight section or other source of pure TE ; the energy fluctuates 
back and forth between a condition of pure TE; and of predominant TMi, . 
The length and magnitude of the fluctuations are derived. 

Section 7 computes the increase in average attenuation caused by serpen- 
tine bends of regular shapes. Numerical examples are tabulated. 

Section 8 shows that the results of Section 7 can be applied to helical 
bends and to small two-dimensional random deviations from a straight 
course. 

Section 9 shows that for any given statistical distribution of random 
angular deviations the average attenuation is minimized by an optimum 
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wave guide radius for each signal wave length and by an optimum signal 
wave length for each wave guide radius. 
Numerical examples are given for sinusoidal bends. 


Summary of Results 


. The energy loss of TEo, waves in curved wave guides by conversion 


into the TM, mode is interpreted as a case of coupling between 
resonant transmission lines. 

In a pair of coupled lines the energy cannot be confined entirely to a 
single line but travels through both in one or both of two possible 
combination modes. 

All important properties of coupled circuits, including wave guide 
modes, are functions of a single discriminant. 


. When the discriminant is much smaller than one, most of the energy 


can be carried in one line or component mode. 


. When the discriminant is much larger than one, the energy flow is 


6. 


~ 
. 


10. 


11. 


nearly equally divided between the two lines or component modes. 
In wave guides of typical dimensions the coupling discriminant 
becomes one for a “critical”? bending radius greater than a mile. 
For all sharper bends, that is for most practical installations, the 
discriminant is greater than one. 

In a long wave guide section with more than critical curvature the 
average attenuation constant is the arithmetic mean between those 
of the TEo: and the TM, modes. 

If a wave guide region carrying pure TEo: is followed by a curved 
region, the energy in the curved region fluctuates back and forth 
between pure TE»; and predominant TM,,. The location of TEo; 
minima and maxima is a function of the signal frequency, the wave 
guide diameter and the total bending angle. 

For highly supercritical curvatures the bending angles at which 
minima and maxima occur are nearly independent of the curvature 
and approach the limiting values previously computed by Jouguet 
and Rice. ‘The minima approach zero. When the bending radius 
approaches or exceeds the critical value, the maxima and minima 
become shallower and their spacing is increased by a function of the 
coupling discriminant, 

For regular serpentine bends or random angular deviations from an 
average straight course which are much smaller than the first extinc- 
tion angle, the percentage increase in average attenuation is propor- 
tional to the square of the maximum deviation and to the fourth 
power of wave guide diameter and signal frequency. 

Wave guide installations of practical dimensions for frequencies now 
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attainable are tolerant to random angular deviations of the order of 
1 degree. 

12. For any expected distribution of random angular deviations there 
exists an optimum wave guide radius for each signal wave length and 
an optimum signal wave length for each wave guide radius, which 
minimize the average attenuation. 


REFERENCES 
1. Jahnke & Emde, Tables of Functions, Dover Publications, New York, 1943. 
2. M. Jouguet, Propagation dans les tujaux courbés, Comptes Rendus—Académie des 
Sciences, Paris, Feb. 18, 1946, March 4, 1946 and Jan. 6, 1947. 
3. M. Jouguet, Effets de la courbure dans un guide a section circulaire, Cables & Trans- 
mission, 1 No. 2, July 1947, pp. 133-153. 
4. S. A. Schelkunoff, Electromagnetic Waves, D. Van Nostrand Company, Inc., New 
York, 1943. 
. M. Wien, Ueber die Rueckwirkung eines resonierenden Systems, Ann. d. Physik, 
1897, Vol. 61, pp. 151-189. 
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ANALYSIS 


1. INTERACTION OF COUPLED CIRCUITS 


1.1 Free Oscillations of Coupled Resonators (Fig. 1B) 
The circuits are coupled according to the following four equations: 


a = —Zith + Zale 1.1-1 
1= YVyei + Ymee 1.1-2 
e@2 = —Zole + Zmli 1.1-3 
a oe 1.1-4 


where index ; refers to circuit 1, index 2 to circuit 2 and index ,, to the mutual 
coupling impedance and admittance. The coupled oscillations have the 
solution: 


. at , t e 
a = Ewe" a Ene’ 1.1-5 
. rat " t 
eo = E..¢€"° + Exe? 1.1-6 


In the limiting case of zero coupling (1, = 0,Z,, = 0) the obvious solution 
shows independent oscillations in the two separate circuits: 
rs , Pit = 
éw = Ailw = Ewe’ 1.1-7 
— " rot 
€x = Kolo = Exe” 1.1-8 
‘The wave impedance A, of the primary circuit is found by dividing equation 
1.1-1 by 1.1-2 
€) Zi 


me ee 
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9 om Se 
Ko = os = = 2 
: 19 / Y, 


By multiplying equation 1.1-1 by 1.1-2 one finds 


Similarly, 


—Z:\VY,=1 1.1-9 


from which one can compute the exponent ~;. In the specific circuits 
shown in Fig. 1b 


Z\ = Lips + R, and 1.1-10 
VY; = Cips 1.1-11 
From 1.1-9, 10 and 11 
° R; ® / 1 RR? 
ee ee oe E ; a .1--12 
pi 1+ Jo oL, + J V LC, aL? 1.1--1 


nd by analogy 


w co ; 1 R; as 
he eee ag 4/ LG a) 8 


In equations 1.1-7 and 1.1-8, Ey) and £29 are amplitude constants determined 
by boundary conditions. In equations 1.1-12 and 1.1-13, 6; and 6» are the 
decay or damping constants, w; and w2 the radian frequencies. 

With finite but loose coupling and small damping the circuits can oscillate 
with either or both of the two frequencies. 


pa EM Md = +05 poll — Vite!) 11-14 

pte = i Pte = +05 pil — 0/1 + 2) 11-15 

In the Jast two equations, the symbol x, defined by x = V > Pips -k, may 
1 


— pe 
be called the coupling discriminant. The first term of the product on the right 
side of this expression is the reciprocal of the fractional difference between 
the uncoupled frequencies; the second term & is the “coupling coefficient.” 

When there is only one coupling impedance, the coupling coefficient is 
usually defined as the mutual circuit impedance divided by the geometric 
mean of the separate circuit impedances. A broader definition which 
applies to all combinations of mutual impedances and admittances is 

Pr» Poy 


Vv P,P» V P,P» 1 16 








et 





CURVED WAVE GUIDES 9 


In this equation P, is the energy stored in circuit 1, Pz the energy stored 


in circuit 2 and Pj) is the energy transferred from one circuit to the other. 
One finds 


2 i2K 2 we : 

P; = psa + "; ~ : = 1K, 1.1-17 

P= 8 2 AK 1.1-18 
2 K, 2 2 es? 

Py = a me Sutg Ks wm Oe +b igi Bs 1.1-19 


; Ky 


Equations 1.1-5 and 1.1-6 contain four amplitude constants. ‘Two of these, 
for instance Fy, and EF», can be adjusted to satisfy boundary conditions. 
The other two are fixed by the equation 


EiKy _ Pe- Pi _ Pe — be _ Eb Ks 


Eig K» Pa — pr fo — fr E>, Ky 


The square root of this expression, 


Ex ,/Ki _ Ew . /Kz 


1 — Ta 
gE Vt.” mY FE, 


may be called the normalized amplitude ratio. It is a vector quantity 
denoting the amplitude ratio and phase relation of each oscillation frequency 
in the two circuits, assuming that they have been normalized to equal 
resistances by an ideal transformer. The absolute value 


| 
Ee ks ‘4 Ws 


is the ratio of the energies stored in the two circuits oscillating at frequencies 
Pa and py» respectively. 
From 1.1-14, 15, 16 and 18 
os (SE SS 
V1i+ e+ 1 
Ae=Vit¢nt—« =A 


When the indexes are left off, W 


=W 


< land| Aj, < 1 by definition. One 
sees that energy, amplitude and phase relations between the coupled 
circuits at each oscillating frequency are governed by the coupling dis- 
criminant. This also applies to the damping coefficients and frequencies 
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of the coupled oscillations. It can be shown by combining and transforming 
equations 1.1-14, 15, 19, that the coupled damping coefficients are 


6, + bo W Wa W'sa 
Oa oo - =0 ; : be — 
1 + VW " W total t HW total 
6; W +- 6b» 6} Wop Oo Wop 
Op =— — 


1 + W W total W total 


The damping constants of two coupled resonances are found by combining 
the uncoupled damping constants in the same proportion as the energies 
oscillating in the two resonators. 

The coupled frequencies are 


a aS Wee and 
“—" T=W 
om We eo Wor 
, i- w 


1.2 Forced traveling waves in coupled transmission lines (Fig. 1A). 
The two lines are coupled according to the four equations 


Te, = Sil, + Zmte 
My 


Te. = Solo + Zml1 


yi) + Vm2 


Tig = Y2€2 + Vme1 
which may be compared to the corresponding equations of section 1.1. 
There is a dimensional difference because in transmission lines the series 
impedances z are measured in ohm/meter and the shunt reactances y in 
mho/meter. I is the propagation constant of the wave traveling in the 
+s direction. Ifa sinusoidal signal with the radian frequency w is impressed 
upon the input of the lines, the coupled waves have the solution 
’ jwt—T ’ jwt—T'ps 
a= Eye” “+ Ene” ” 1.2-1 
+ _jwt—Tgs + jut—P 
& = Ene” " + Ene” 1.2-2 
For zero coupling one finds, in analogy to Section 1.1 
oe 3 . jot—T 
€i0 = Ayt, = Hye si 


. i 
ey = Kot: = Ene” aie 


) 
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and 
r, = V 21 
Tr, = V yo % 
Ey and £2 are independent integration constants. For finite but loose 


coupling and small attenuation constants one finds in analogy to 1.1-14 
and 1.1-15 


r I> T as T. Pp 
ra PES Vite = +051 - Vite) 
a a ee oe a 
ae as PO Ovi +e = +051 - Vite) 
where 
k ; 
= ‘PVs 2 
K | ia I's av ] 1 I 2 1.2-3 


is the coupling discriminant. Just as in Section 1.1 the coupling coefficient 
k is defined by the equation 
Pr» Po 


-aese 1.2-4 


VPiP, VPiP2 


P, is the energy per unit length stored in line 1; P:, the energy per unit 

length stored in line 2, and Piz = P21, the energy per unit length inter- 

changed between the lines. The waves can travel in the coupled lines with 

either or both of two transmission constants. Two of the amplitude 

vectors in equations 1.2-1 and 1.2-2, for instance Fy, and Ex, are free to 

satisfy boundary conditions; the other two are determined by the equation 
_re—-hh —Ts_ Enks 


Eta ky : ’ 
Exa Ke ls — rT, r, a r; Ex Ky 


Ex 4/® = Ew s/= = A. = ae 1 
Ei Ks Ex K, ws Ap j 


A is the normalised amplitude and phase ratio for two lines transformed to 
equal wave impedances. 


bho 
an 


Ezy Ky _ wg 
| Eta Ke " We 
~ i Vite ri _ . 
W.= = = W 
IVitett 


is the ratio of energy flow in the two lines. At the propagation constantT,, 


A, = V1 + Kn? — ~ a A 1.2-6 
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When the indexes are left off, W < 1 and | A| < 1 by definition. 
In a manner analogous to that of Section 1.1 it can be shown that the 
coupled attenuation constants are 


ay + ayW Wa Wea - 

a= —— = — .——— 1.2- 

we Tee Oe ee 
._< a» + a) HW a WW 1b + ™ WW 2b 1.2-8 


* Weta 


iT.” 


The attenuation constants of the coupled waves are found by combining 
the uncoupled attenuation constants in the same proportion as the energies 
traveling in the two lines. 

The coupled phase constants are 


Bi — WB, 
= ' é ° -9 
Ba i-Ww and 1.2-9a 
Bs — WB, a4 
= in 2-9 
Bs 1 <a W 2: b 


From equations 1.2-5 to 1.2-8 one sees that the coupled propagation 
constants are conveniently described in terms of the power ratio W. W 
itself is a known function of the complex coupling determinant «x which is 
shown on the attached Fig. 4 for the following three special cases: 

Case 1. 

The two lines have equal phase constants and different attenuation con- 
stants: B2 = Bi a Say 

k is an imaginary number. 

W changes its character abruptly at the critical coupling. 

| x critical | = 1 

For|«| <1 

W<1; a2aa; Bs=Be 
For|«x| 2 1 
W = ites Ap = Aa; Bo 2 Ba 
Case 2. 
The lines have different phase constants and equal attenuation constants. 
x is a real number 
W changes asymptotically from 
Wo = 0 to 
W, = 0.172 and to 


W. = 1 


II 
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Case 3. 

The phase and attenuation constants differ by equal amounts. As shown 
below, in section 4, this case applies to the coupling between the TEo; and 
TMi: modes in curved circular wave guides with finite conductivity. 

«? is an imaginary number. 

W changes asymptotically from 


Wo = @ to 


W, = 0.217 and 
Wex=l 
For « > 1 all three cases approach the limit 
» « Atk & - ) 
bce 9 (1 + D 1.2-10 
. a + am +k . a + a 


2. DERIVATION OF FIELD EQUATIONS 


Consider a straight circular cylinder with an inside radius such as shown 
on Fig. 1 A. Let the radial coordinate equal r, the azimuthal coordinate 
equal ¢ and the longitudinal coordinate equal z. Let the dielectric losses 
inside the cylinder be negligible. 

The field equations inside the cylinder are® 


dE, dE, 


= —iwu H, 
raps eo 
OE, o8, P 
ac ar = —iwur H, 

d(rE,) OE, _ . 
‘a =<" jour H, 
and 

OH, ” OH, wdiadlt 
ray 02 ‘ 
oH, _ 0H, sa sal 
az or 

0(rH,) _ OHr «fil, 
or 0g 


® See Ref. 4, pg. 94 of the book. 
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The natural transmission modes which satisfy these equations have the 
form 


E ae frlr) ereteo) , F atin 0 


Each of these modes conforms to the same equations as a wave traveling 
in a transmission line with an impedance and phase velocity dependent upon 
the mode. In a straight cylinder with perfectly conducting walls, there 
exists no coupling between the different modes so that any and all can 
exist without interacting. If the conductivity of the walls in a straight 
circular cylinder is finite, it produces a resistive coupling between modes of 
equal azimuthal index ( in equation 2-0), In copper tubing and at the 
frequencies now obtainable (w < 10") this coupling effect is negligible. 

A stronger coupling may be caused by deviations of the wave guide from 
the shape of a straight circular cylinder. The deformation considered 
in the present analysis consists in a circular bend of the axis, as shown 
schematically on Fig. 1b. 

In such a circular bend the longitudinal coordinate is replaced by the 
product of the bending radius R by the bending angle @: 


z= R06 
This transforms the first two component equations of curl £ into 


A(REs) dE, 


~~ ne 
dE, ARE). 
nun” So 


The variable R can be eliminated by the relation 
R= Rv) — rcos¢g 


where Rp is the bending radius of the cylinder axis. The coordinate @ can 
be replaced by a longitudinal coordinate s, measured along the cylinder 
axis. Hence, 


5s = ORo 
The progressive modes which we investigate have the approximate form 


E ~ filr) ” ied tee! err 


Hence 


0 
. Os Ro 


A ; : 
im = jn for all field components. 
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wn 


0 F 
ay «may be expressed by a prime: 
r 


oF = Ff 
or 


Thus the equations with curl E become 


me, E, sin ¢ RTE : 
J + ints In ¢ bee BPs. ie ¢ = — wp H, 
r Ro—rcosg Ry —rcos¢g 
—R TE, v E, cos ¢ . 
—— — E,- ——___—_—_ = — H 
Ry — rcos¢ . 2 sale a 


For gradual bends 


Ry>a>r 


One may therefore approximate 


o.oo i+; cos ¢ 
Ro — rcos¢ Ro 
[t is convenient to introduce the symbol 
a 
Ro 


which is proportional to the coupling coefficient. All powers of ¢ greater 
than the first will be neglected. One can now write the approximate field 
equations in the curved cylinder: 


_ + TE, + < E, sin g + cI E, cos g = —jwu H, 2-1 
Th ~~ eT - E, cos g + “ E, cosy = —joun H, 2-2 
E, + rE, — jnE, = —jourH, 2-3 

a= + TH, + “H, sing + oF : H, cos y = jue E, 2-4 
~fti, + ian “ H, cos ¢ + < H. cos ¢ = jue Ey 2-5 
H, + rH, — jnH, = jwer E, 2-6 


The coupling terms all contain the factor cos ¢ or sin gy. This means 
that every transmission mode is coupled only to modes with an azimuth 
index differing from its own by +1. 
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3. CHARACTERISTIC EQUATIONS OF TEo; AND TM,, Mopes 


The mode which one desires to propagate through the wave guide is the 


TEo: mode. 
it is characterized by the following equations: 
n= 0 
En = Ea = He i 
Ean = Eye" Ii(y) = ev Ji(y) 


ey T; 
H, ae Ji(y) 
JBon =~ 
e€ 
Ha = — 2% Joly) 
JBon 
In these equations 
7 = /'u = 377 ohm (intrinsic free space resistance) 
€ 
— 2a 
Bo = wv eu = — 
Xo 
7 ey 
Mh = Vx’ — 6 
In a perfectly conducting wave guide 
? 
x, — 3832 
a 
.O1X 
y= - = —— (cutoff factor) 
Do a 


lr, = jhV1 — v= 7B: 


In a straight wave guide with perfectly conducting walls 


3-1 
3-2 
3-3 


3-4 


3-12 


If the wave guide has the conductivity of g mho/m, its intrinsic high 


frequency impedance is 


Ze= (1+ HR = (1 + j)34d 4/4 
&Xo 
This changes x to 
IX = xo — 3 = he and 

na 
; v° (1+ j)Rit 
Tr = / ~ 
1 IB: + J/1 — P on 
* Ref. 4 pg. 83. 
t Compare Ref. 4, pg. 390. 


3-13 


3-14 


3-15 
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Due to the curvature of the guide, this desired mode is coupled to all modes 
which have the azimuthal index number 1. 

However, for low curvatures, this coupling is very loose and only causes 
appreciable effects if it can act over a great length of wave guide without 
phase interference. 

This means that the disturbing mode must have nearly the same phase 
velocity as the desired mode. It so happens that in a perfectly conducting 
circular cylinder there exists one mode, the TM, which has exactly the 
same velocity as the TEo:. Such a coincidence is called “degeneracy.” 

In the analysis of very gradual bends, only this TMi; mode need be 
considered. It is characterized in a straight guide by the following equations: 


6 = } 3-16 


Re? . ae cos (y + ¢o) = e a) cos ¢ 3-17 


Ea = 
The TMi: mode can be polarized in all directions. But since only the 
component directed toward 
¢, = 0 
is excited by the wave guide curvature, go has been omitted in the last 
term of eq. (3-17). 


dJ,(y) dj 
e2 —— 
dy 


Ey sing = @J,(y) sing, where J = a 


II 


— x2 : 
E,« Ji(y) sin ¢ 
Yr, 
Ha = 2B Ji(y) sin ¢ 
nl's 





2580 Jy) 


Hu = 
Ais nl's 


os y 
Hea = 0 
In a perfectly conducting wave guide the x defined by eq. 3-8 is 


X = x2 = xo and 


r, =; = JBi 
In a wave guide with an intrinsic impedance per 3-13, 
x = xo — o-2E ol 3-18 
nav 
(1 + J)Ri t 
T, =7 ————— 3-19 
2 jBi + +/i =~ wen 


From 3-15 and 3-19 one finds a, = agv” 3-20 
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4. INTERACTION BETWEEN TE: AND TM), Mopes 

Since the separate modes of propagation behave like traveling waves in 
transmission lines’, their interaction can be derived from the coupling 
equations derived in Section 1 of this analysis. 

The uncoupled propagation constants [; and I’; are known (equations 
3-15 and 3-19). In order to find the coupling discriminant one must derive 
the coupling coefficient from the field equations 2-1 to 2-6. The coupling 
coefficient is defined as 

k Py» Px ( ) 
= } = e eq. 1.2-4 
V P,P» V Pi P2 

In computing the coupling coefficient one may neglect the small attenua- 

tion constant. The energy stored by the TEo: wave per unit length is 


+2 
“1 


a FE 4 
P= | Pande = [iin 2ur dr 41 
0 7 ° 


This expression is not affected by the cutoff factor y because one may 
consider the field inside the guide as composed of slanting plane waves with 
the electric field strength £,. The energy stored by the TM), wave per 
unit length is 


27 pa E 27 pa 5 
P, = | —rdygdr = [ [ Honrdgdr The inter- 
0 40 7 0 0 


changed energy: 


nlx a Eb E 2x a 
Pas | | = dy dr + [ | Ho», Hn dy dr 
0 0 n “0 0 


Combining equation (4-1) with 3-3, 3-8 and 3-10 


2ra e} 


P, = A 
* ~~ 3.8322» 


3.832 : 
| yoi Ji(y) dy 
0 


From reference 1, page 140 of the book 


+3. 832 2 
2 3.832" 
i yoy Ji(y) dy = > J (3.832) 
0 2 
Hence 
0.51e; a" 
a Sle a 4.2 
n 


In a similar manner one finds 
a8 2 
O.51e5 a 
i 
2(1 — v*)n 


2 Loc. cit. 
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and 


0.51ce,e:a° 
Py = Py = — 4-4 
_ 

Substituting the values of 4-2, 4-3 and 4-4 into 1.2-4 one finds for the 
coupling coefficient between the two groups of plane waves traveling at a 
slant to the wave guide axis 

b, = VIVI=- # _ 03690 ——, 
: 3.83 Ro ' 

The coupling coefficient k between the TEo; and TM; modes which are 
the resultants of their slant wave groups is greater than &, according to the 
following reasoning: 

From 1.2-10 

B = B(1 + 0.5 &,) 
This makes the cutoff factors of the coupled modes 
D4 Vo 
= = e 
B 1 + 0.5k, 
and the coupled propagation constants 
T= pV1 — 2 = Pov (1 + 05k)? — v2 
For k, < 1 


ra avi=a(ts 


0.5k. : 
.) = Tol — 0.5k) 
1 > Te 


Hence, in view of eq. 1.2-10, the effective coupling coefficient of the wave 
guide modes is 


k k __ 0.369a an 
~ 1 RhVIi-# e 
From equations 3-15 and 3-19 
rT, —T, = es pete 4-6 
an 
Vir Pee... ea ee... 4-7 
i a Fp ee — j)RiV1 -F 
From 1.2-3, 4-5 and 4-7 one obtains the coupling discriminant 
— 0.369-29-377a° 
c= 
RoR — 7) — v’) 
Since 
R; = 0.00452 Xo": p, * 4-8 


* Loc. cit. 
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with the relative high frequency resistivity 


Pr = / p 
Peopper 


_ __ 9.65-10°(1 + fan ’* 























































































































Ro pr(1 or v’) 
Its absolute value 
5 2, —0.5 
| l 1.366-10°a 
K = —_— -— oe 
Ry p-(1 a vy) 
1.2 
1.4 CASE | 
EQUAL PHASE | 
“a a CONSTANTS y 
| wna n ae === 5 po-L-}--b =p. 
¥ en 
a See Se ee 
0.9 = | eel 
I va ae 
0.8 | Pa ly 
: ai T ore a ak 
a } \_CASE 2 
Q 0.7 1 EQUAL ATTENUATION 
te 4 Ps CONSTANTS 
a 0.6 i 
& 0.5 CRITICAL, _ - re =e 
x COUPLING * Pd 4 
& 0.4 4 f. = 
1 4 | 
- ! ye CASE 3 
y Ad ~>=-~-~-TEo, AND TM, MODES: 7 
0.2/— Vi EQUAL DIFFERENCES BETWEEN PHASE | 
‘ 1x AND ATTENUATION CONSTANTS Ez 
ay 
0.1 A | 
: | ch nicadbemindtibel 
0.1 0.2 0.4 0.6 { 2 4 6 8 10 20 40 60 100 
COUPLING DISCRIMINANT, |[K| 


Fig. 3 


As shown in equation 4-6, the differences between the propagation 
constants of the TMy; and TE»; waves are proportional to the intrinsic 
skin impedance R; which contains the factor (1 + 7). This means that the 
phase and attenuation constants of the two waves differ by equal amounts 
in accordance with “Case 3” of Section 1. For this case the power ratio 
per 1.2-20 may be written 


a l\Vitjl«P-1 
ee erie re 


The numerical value of this function is plotted on Fig. 3. It can be 
computed conveniently by means of the following auxiliary parameters: 


V1+j) «x? = p+ jq = cosh x + j sinh x 4-9 
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p= V0.5 + 0.54/1 + [x = cosh x 
V-0.5 + 0.54/1 + |x|! = sinh x 


s 
\| 


, —- | x 
W = p tanh - 
p 2 
In analogy to Case J of Section 3, the condition of | x | = 1 may be called 
critical coupling. It occurs at the critical radius of curvature 


1.366-10° a’) °° 
p(1 — 2”) 


For subcritical coupling (Ry > Ri:) W approaches 


Ro ed meters 4-10 


. ‘ K 
W suber. = t —— > 0) 
lor supercritical coupling (Ro < Rer), 
ii 
W esaperer. =] — \ — » I 
| «| 


Krom the above results, it is possible to predict the behavior of waves 
originating either as TE; or as TMi: modes, in any given wave guide 
configuration. This is done for some typical cases in the following sections. 


5. PROPAGATION IN LONG WAVEGUIDES WITH CONSTANT CURVATURE 


It has been shown in Section 3 that for each curvature there exist two 
modes of propagation. 


In one, 
W sit Pom < 1 5-1 
. P reo 
and the attenuation 
“a < QTE 5 QT™ 5-2 
In the other, 
; 1 ss 
W; == >I1 5-3 
We 
and the attenuation 
4 > oe tim a 
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In a long wave guide the “6” mode will die down due to its greater attenua- 
tion, no matter how much of it was initially present, so that one need only 
consider the ‘“‘a” mode. 

This mode has a phase velocity slightly smaller than that of the uncoupled 
TEo: wave and an attenuation nearer to that of the TE»; than the TMi, 
wave. 

The magnitude of the critical radius is illustrated by the two examples 
of Table I. 

TABLE I 
CHARACTERISTIC VALUES 


Parameter Symbol Equation Example 1 Example 2 

Wave guide radius a .0S5 m .O5 m 

Free space wave length Xo .03 m 01m 

Cutoff ratio v 3-11 . 366 £22 

Attenuation constant QTE (cu) 3-15,4-17 2.04 10 3.58 K 10° 
neper/m neper/m 

Attenuation constant apM(eu) 3-19 1.36 10% | 2.41 xX 10° 
neper/m neper/m 

Critical Radius Rerit 4-10 2.12 km. 3.44 km. 


TABLE II 
RELATIVE ATTENUATION VERSUS RADIUS OF CURVATURE 


General formulae Example 1 Example 2 

K Ro/Rer W a/ao Rokm a/ao Rokm a/ao 
0 20 0 1 x 1.00 co 1.00 
0.1 10 0.0025 1 + 0.0025 (v? — 1) 19.7 1.02 | 34.15 Oe IY 
0.2 5 0.01 1+ 0.01 a 9.85 1.06 17.08 1.66 
0.5 2 0.06 1 + 0.057 ‘i 3.84 1.38 6.83 4.45 
1 1 0.22 1+ 0.18 Js 1.97 | 2.16 3.42 | 12.94 
2 0.5 0.48 1 + 0.32 3 0.98 3.11 4k | 22.2 
5 0.2 0.75 1 + 0.43 0.38 3.83 0.68 | 26.6 
10 0.1 0.87 1 + 0.46 ee 0.20 4.05 0.34 27.9 
20 0 1.00 1+ 0.50 2 0 4.30 0 34.1 


The increase of attenuation in long wave guides with uniform curvature 
is shown on Table II, with numerical values for the same examples as in 


Table I. 


6. PROPAGATION IN A UNIFORMLY CURVED SECTION OF WAVE GUIDE 
FOLLOWING A LoNG SrrRaiGuT SEcTION. (FIG. 4) 


No matter what mixture of modes may prevail at the beginning of the 
wave guide, all modes except the TEo; die down in the long straight section 
due to their higher attenuation, so that the wave form at the beginning 
of the curved section is pure TEp: . 

Since it has been shown in Sections 1 and 4 that each of the two possible 
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modes of propagation in a curved wave guide consists of both TE»; and 
TM), waves, it follows that both modes must be superimposed in such a 
manner that at the transition point the TM components cancel each other by 
interference. 

Let the relative amplitudes of the two TE components equal a and ); 
then the corresponding amplitudes of the TM modes are aA, and DA, 


Y 





Ro 





2 s 
o—! 











Fig. 4 


respectively, where Ag and A, are the normalized voltage ratios per 1.2-5. 
At the beginning of the curved section 
a+b=1 (TE amplitude) 
aAqg + bA, = aAqg — b/Aa = O (TM amplitude) 


Hence 


6 = 


b 
a 
‘The two waves have different phase velocities and therefore interfere 
with each other. According to 1.2-9, the difference between the phase 


constants is 

) i+w 
1-W 

By means of equations 4-9, this can be transformed into 


Bo — Ba = (82 — Bi)(p + Q) 6-1 


By — Ba = (B2 — By 
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The length of one complete interference cycle is 
2r 
(p+ g\(G: — Bi) 
If both components had equal attenuation, the beats superimposed on the 
decaying envelope of the TEo: wave would correspond to an amplitude ratio 


eomax 1+W _ 
eomin 1—W nee 


However, during the progress of the mixed wave through the curved section, 
the intensity of the fluctuations is reduced by the greater attenuation 
constant of the faster and weaker “‘b” component. In one complete inter- 
ference cycle, the differential attenuation reduces the weaker component to 


Aes (a9—a@))8oe (24/p+q) 


oe = ¢ 
Ay 


Approximation for Weak Coupling 
For|«|<«1 
By — Ba = (Bz — B1)(1 + 9.5 | x |?) 


From 3-15 and 3-19 


R; 
& — Be =—vVi — 
an 


ss 4/ Ae (1 + 0.5 | « (?) radian/m 
For intermediate coupling, 6-1 may be transformed into 
By — Ba = kBif(x) 
with fi) = an a v1 + v1 + |x + v= T+ V1+ |x} 


V2 | «| 2| «| 
Approximation for Strong Coupling 
For |«|> 1 
fin = 1+ 0.125 | «| By — Ba = BBL + 0.125 | « |~4) 
Substituting the value of & from 4-13 and transforming, 


bs — Bp = 2320, 2.820 
b a hy Ro (Kx ho Ro 


The phase difference between the two components is 


(1 + 0.125 | « |“) 6-2 


2.5285. , 2.3200 . 
Wo = i (.) = M6 6-3 
. AOR 2 Ao J ; 
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where @ is the bending angle of the wave guide. The power carried by 
the TE»: wave is 
2 Me 


6-4 
2 


Pre = cos 
Minima of TE: occur when this phase difference is an odd multiple of x. 
Hence, the bending angles producing minima of TE; amplitudes are: 


1.36X9 _ (Qn + 1)2.22y 


1.360 | 6-5 
a(1 + 0.125 | x |-*) It) 


Ormnin = (2n + 1)- 
The initial fluctuation ratio approaches 


€1 max : | | 
= pta* V2 ln 
€1 min 
which is a large value tending to infinity. 
The relative attenuation of the slightly weaker component during one 
beat cycle is 
Aoy __  _—24/p+q Jeri inl 4.44 
€ = 1— 
Ao kK 


which is a small reduction tending to zero. Hence, the fluctuations persist 
through a large number of beats. The power is transformed back and 
forth between the TE»; and the TMi; modes. 

In Section 5, it was shown that in a long, uniformly curved wave guide 
the attenuation is intermediate between that of the TEo: and TMi: modes. 
But from equations 1.2-7 and 8 it follows that the two modes contribute 
to the attenuation in proportion to their relative power flow. Since at 
the beginning of the bend the power of the TMi; component is zero, 
it is to be expected that the initial rate of attenuation equals that of the 
TEo: wave alone. This is proved by differentiating with regard to s. One 
finds for all values of x that 


d T a8 Ty, 
—jae * + be” = —ay 
ds a=0 


Discussion of Results 

Equation 6-2 corresponds directly to an equation derived by S. O. Rice 
and, after allowing for the different choice of variables, to M. Jouguet’s 
equation (75)°. It differs from the results of these earlier calculations 
by the factor f(x) = 1 + 0.125 | x | which is a reminder that the simplified 
form of the equations given by the earlier authors is an extrapolation to 
infinite conductivity or infinite curvature of the wave guide. 


6 Reference 3, pg. 150 of Cables and Transmission, July 1947. 
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From equations 6-3 and 6-4 it is seen that the TEo: wave is recovered by 
bends which are an even multiple of Onin. But such bends are efficient 
transmitters of TEo: waves only over a narrow frequency range since 
Omin Varies with frequency. 

If the circular bend is followed by a long straight section, the TE; and 
TMi, components existing at the end of the bend are carried over into the 
straight section, but the TMi component dies down due to its greater 
attenuation and constitutes a total loss. 


Numerical examples for first extinction angle. 


Using the same dimensions as in Table I of Section 5, one finds from 
eq. 6-5 for: 


Example 1: @n:, = 0.816 Radians = 46.8° 
15.6° 


Example 2: @ni, = 0.272 Radians 


7. SERPENTINE BENDS 


Sections 5 and 6 dealt with bends continued with uniform curvature 
over large angles. The present section considers the small random devia- 
tions from a straight course which are unavoidable in field installations. 

Actual deviations are expected to be random both with regard to maximum 
deflection angle and to curvature; they are likely to approximate a sinusoidal 
shape. For purposes of computation, the following analysis assumes as a 
first case circular S-bends which consist of alternate regions of equal lengths 
and equal but opposite curvatures. An exaggerated schematic of such 
S-bends is shown on Fig. 5A. 

Each circular bend tends to produce a single mode with an attenuation 
per equation 1.2-7. However, the discontinuous reversals of curvature 
at the inflexion points produce mixed modes, and the initial part of each 
region reduces the amplitude of the TM components produced in the 
previous regron. 

Each region may be treated as a discrete 4-terminal section of a trans- 
mission network. Regardless of the wave composition at the input terminal, 
differential attenuation will establish in a long serpentine wave guide a 
steady state condition. In this steady state each region produces equal 
attenuation. This attenuation per region and the resulting average 
attenuation constant will now be derived. 

The TEo; and TMi waves each consist of ta” and ‘*b” components with 
separate amplitude ratios and propagation constants, as derived in Sections 
1 and 4. In the first region (between points 1 and 2 of Fig. 5A) Ro is taken 
as positive, and 


Az=zVY14+e-4 1.2-6 
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to 
~I 


In the second region, between points 2 and 3, the polarity of Ro, and 
consequently of x and the ratio of TM to TE amplitudes, are reversed. 











b- SINUSOIDAL BENDS 
Fig. 5 


Except for this change of polarity the amplitude ratios at points 1 and 2 
are equal. Introducing the symbols 


—T'sm 


g 
£m = (ga + Zo) 


ga = 43(go — a) 
one can tabulate 
point Crk Cr 


1 ad + b aA ne 7-4 


é l 
gm (! + be) gm e : ace «) 7-2 
8a 2a A 


ie ale ; 
Calling — = y, one finds 
a 


to 


1+ yee _ —A° + ye _ ba 


_ Saverage 


i+; 8-35 "sO so 
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1 i+y 
Lave rage ~— yr roe iz 3 2 ] 
averag re. + mh lo i+” with 


:, are = ue Ke —2\2 
y= x m6] + ga) +A / a + _ (1+ ga)” 


This formal solution is hard to evaluate. It can be greatly simplified 
for the subcritical and supercritical cases. 


1. Subcritical Curvature |x| &1 


re 
eae 
D'average _ I, (1 + - : a «) = Is = lr, 
“1+ ga 


For very low curvatures, the average attenuation approaches that of the 
‘a’? mode, and this in turn approaches that of the TE; wave. 
2. Supercritical Curvature. |x| > 1 


The differential attenuation constant is small compared to the differential 
phase constant. 


y | - | A | aad 1 
Substituting these values into 7-1 and 2, one finds 
y = e 0.570 m 


Expressed as a function of 6: 
ve = cosy + 7 siny with 


y = M(6 — 0.5 6,,) 7-3 
M has the value per eq. 6-18. 
The power ratio of the combined TM); and TEo: waves is 
W, = tan? y/2 
In view of equation 1.2-7 the instantaneous rate of energy loss is 
ag = Qa cos’ ¥ -f- ay sin’ ¥ = dQ) + (a2 <= a) sin’ : 7-4 
1 + 8m 1 Om 
Quveruge = — a, d, — oe ae dé 
Sm 40 m “0 
1 sin M6, 
Qaverage — Q@ + ag~- a — — 
rag 1 ( 1) (; 2M0,, 


In view of 3-20 


7 y* +1 » ¥ — 1 sin M6, 
Qaverage = M1 > 2 M?é@, 
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For small deflection angles, M0 <K 2 
y? — 1 9 2 
Qaverage = a E + — 6 M a | 
If a p% increase in attenuation is the tolerance limit 


Om & . a Substituting the value of M from 6-3, 
v 





~ 10M 2-10 
0. = 0.105-V/p vro 
~ i=? 


The maximum deflection equals 


Ae = 0.50,. Hence, in view of 3-11 


' 2 : . 
0.032V'P adians — L840 P degrees 


> 


avV/1 — v a? ll a 


3. Sinusoidal Bends with Predominantly Supercrilical Curvature. 


Ag 


Sinusoidal bends cannot be supercritically curved over their entire length 
because at the inflection points the curvature drops to zero. For sufficiently 
short bends, however, no great error is caused by treating the entire length 


as supercritical. In that case, equations 7-3 and 7-4 remain valid. 6 


takes the new value 
Ci. 6 


> sin aos Sm) 


Hence 
— MOn is a(s — Sm) 


8m 
aQ2—~ ad) +9 M6 P ms —s 
Qaverage = 1 + sin | ™ sin ») ds 
S 0 2 25m 


Sm 
For small deflection angles, M@,, << 2 


22 8 
d= ad M Om 28 3 2 a(s coi Soa 
sin ; ds 


average = On Sen 4 th 25m 
M’6;, +. 
= a1 + (a, — a) =a}/1+ 3 t. M’6,. 
4 4 : 
).026A9 ' A9N9 
= ie ve radians = a "4 P , degrees 
avi-Yr a aa 


The tolerance limit for sinusoidal deflections is 20% smaller than for 


circular S bends. 
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The effect of supercritical but shallow circular and sinusoidal S bends is 
illustrated by the following numerical examples. 


TABLE III 
INCREASE OF ATTENUATION IN S BENDS 


Maximum deflection 46 (in degrees) 


Attenuation 
Increase p% Example 1 (v = .366) Example 2 (v = .122) 
Circular Sinusoidal Circular Sinusoidal 

10 2.25 1.82 0.23 0.19 
20 3.18 2.58 0.33 0.27 
30 3.89 3.15 0.41 0.33 
40 4.50 3.64 | 0.47 0.38 
50 | 5.03 4.07 


| 0.52 0.42 


8. HeLticAL BENDS AND RANDOM Two-DIMéNSIONAL DEVIATIONS 


A helical bend may be treated as a bend which has a constant absolute 
magnitude, but a changing direction of curvature. As indicated in eq. 3-17, 
the TMi wave can be polarized in all directions. At any differential 
element of wave guide length, the TM: component polarized in the local 
bending plane is coupled to the TEo: wave; the TMi: component polarized 
at right angles is not coupled and persists unchanged. By requiring that 
the absolute magnitude of the TM),/TE;: amplitude ratio remain constant, 
a steady state solution can be found. 

Shallow helical bends of small curvatures may be treated as the super- 
position of two sinusoidal bends offset by 90° in the longitudinal direction 
and in the bending plane. The increases in attenuation due to these two 
sinusoidal bends are computed from eq. 7-5 and added. 

It is believed that random deviations from a straight course approach 
sinusoidal shape more closely than circular shape, hence equation 7-5 may 
be used to establish a tolerance limit for such random deviations. For 
quantitative results the statistical distribution of the squared deviation 
maxima must be taken into consideration. 


9. OPTIMA OF WAVE GUIDE Rapius, SIGNAL WAVE LENGTH AND 
ATTENUATION AS A FUNCTION OF ANGULAR DEVIATION 

In a straight wave guide the attenuation decreases with wave guide radius 
and signal frequency. However, the deterioration due to wave guide 
curvature increases with wave guide radius and frequency. Hence, for a 
given tolerance limit to angular deviation from the straight course there 
exists an optimum radius for each wave length and an optimum wave length 
for each radius. This will be shown for the case of uniform sinusoidal bends, 
under the simplifying assumption that the cutoff radio vy < 1. 
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Solving 7-5 for p one obtains 


p = 0.45A°a* 9-1 


where p si the percentage increase in attenuation, and A the deviation angle 
in degrees. 


Hence the average attenuation 
as = a(i + 0.01p) 9-2 


From 3-15 and 3-11 


2 


a; 
Introducing the R; value from 4-8 
a = 4.5 10-5p,\!5a-* 9-4 


where p is the high-frequency resistance of the wave guide relative to copper. 
From 9-1, 2 and 4 


ag = 4.5 10° “pA! 5a “8(1 + gra") 9-5 
with 


gq = 4.5 10°A? 
The attenuation reaches a minimum when 
fd, a) = MPa + gd? "a = minimum 
Case 1. 2 is given 
5f/5a = —3\' 5a + gh? = 0 
Gope = 1.32Agq~°-> = 5.24A-9-5 


From 9-5 
Qsom = 4a = 1.29 107"pd8A15 
Case 2. a is given 
6f/65X = 1.5d°- 5a — 2.5g\*5a = 0 


Nopt = 1.14aq’-” = ().2904aA°-5 
From 9-5 


Asopt = 1.6a@ = 1.15 10 ty a 


Numerical Example 


Let A = 0.42° 
a = 0.05 m 
A = 0.01 m 
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From Table III 
ay = 1.50a@ = 5.4 10-'p neper/m 

Case 1: X fixed at 0.01 m 
dope = 0.08 m 
Qopt = 2.76 10-59 neper/meter 

Case 2: a fixed at 0.05 m | 
Nopt = 0.0097 m 
opt = 5.36 10-5p neper/m 


Assuming sinusoidal bends with a 0.42° maximum deviation, the attenua- 
tion of centimeter waves can be reduced to one half by increasing the wave 
guide radius from 5 to 8 cm. For a 5 cm wave guide radius, 1 centimeter 
wavelength is close to the optimum. 























A New Type of High-Frequency Amplifier 
By J. R. PIERCE and W. B. HEBENSTREIT 


This paper describes a new amplifier in which use is made of an electron flow 
consisting of two streams of electrons having different average velocities. 
When the currents or charge densities of the two streams are sufficient, the 
streams interact to give an increasing wave. Conditions for an increasing 
bre and the gain of the increasing wave are evaluated for a particular geometry 
of flow. 


1. INTRODUCTION 


N CENTIMETER range amplifiers involving electromagnetic resonators 
or transmission circuits as, in klystrons and conventional traveling-wave 
tubes, it is desirable to have the electron flow very close to the metal circuit 
elements, where the radio-frequency field of the circuit is strong, in order 
to obtain satisfactory amplification. It is, however, difficult to confine the 
electron flow close to metal circuit elements without an interception of elec- 
trons, which entails both loss of efficiency and heating of the circuit elements. 
This latter may be extremely objectionable at very short wavelengths for 
which circuit elements are small and fragile. 

In this paper the writers describe a new type of amplifier. In this ampli- 
fier the gain is not obtained through the interaction of electrons with the 
field of electromagnetic resonators, helices or other circuits. Instead, an 
electron flow consisting of two streams of electrons having different average 
velocities is used. When the currents or charge densities of the two streams 
are sufficient, the streams interact so as to give an increasing wave. Electro- 
magnetic circuits may be used to impress a signal on the electron flow, or to 
produce an electromagnetic output by means of the amplified signal present 
in the electron flow. The amplification, however, takes place in the electron 
flow itself, and is the result of what may be termed an electromechanical 
. . 1,2 
interaction. 

While small magnetic fields are necessarily present because of the motions 
of the electrons, these do not play an important part in the amplification. 

‘Some electro-mechanical waves with a similar amplifying effect are described in 
“Possible Fluctuations in Electron Streams Due to Ions,” J. R. Pierce, Jour. App. Phrys., 
Vol. 19, pp. 231-236, March 1948. 

2 While this paper was in preparation a classified report by Andrew V. Haefi entitled 
‘*The Electron Wave Tube—A Novel Method of Generation and Amplification of Micro- 
wave Energy’’ was received from the Naval Research Laboratory. Dr. Haeff’s report 
(now declassified) contains a similar analysis of interaction of electron streams and in 
addition gives experimental data on the performance of amplifying tubes built in ac- 
cordance with the new principle. We understand that similar work has been done at the 
RCA Laboratories. 
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The important factors in the interaction are the electric field, which stores 
energy and acts on the electrons, and the electrons themselves. The charge 
of the electrons produces the electric field; the mass of the electrons, and 
their kinetic energy, serve much as do inductance and stored magnetic 
energy in electromagnetic propagation. 

By this sort of interaction, a traveling wave which increases as it travels, 
i.e., a traveling wave of negative attenuation, may be produced. To start 
such a wave, the electron flow may be made to pass through a resonator or a 
short length of helix excited by the input signal. Once initiated, the wave 
grows exponentially in amplitude until the electron flow is terminated or 
until non-linearities limit the amplitude. An amplified output can be ob- 
tained by allowing the electron flow to act on a resonator, helix or other 
output circuit at a point far enough removed from the input circuit to give 
the desired gain. 

There are several advantages of such an amplifier. Because the electrons 
interact with one another, the electron flow need not pass extremely close 
to complicated circuit elements. This is particularly advantageous at very 
short wavelengths. Further, if we make the distance of electron flow 
between the input and output circuits long enough, amplification can be 
obtained even though the input and output circuits have very low imped- 
ance or poor coupling to the electron flow. Even though the region of 
amplification is long, there is no need to maintain a close synchronism 
between an electron velocity and a circuit wave velocity, as there is in the 
usual traveling-wave tube. 

A companion paper by Dr. A. V. Hollenberg of these laboratories describes 
an experimental “‘double stream” amplifier tube consisting of two cathodes 
which produce concentric electron streams of somewhat different average 
velocity, and short helices serving as input and output circuits. No further 
physical description of double stream amplifiers will be given in this paper. 
Rather, a theoretical treatment of such devices will be presented. 


2. SIMPLE THEORY 


For simplicity we will assume that the flow consists of coincident streams 
of electrons of d-c. velocities #; and 2 in the x direction. It will be assumed 
that there is no electron motion normal to the x direction. The treatment 
will be a small-signal ‘or perturbation theory, in which products of a-c. 
quantities are neglected. M.K.S. units will be used. All quantities will 
be assumed to vary with time and distance exp j(wi — Bx). The wavelength 
in the stream, A, , is then related to B by 


p= 2m/de (1) 
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The following additional nomenclature will be used: 


€ dielectric constant of vacuum 
€)9 = 8.85 X 10-” farad/meter 
n charge-to-mass ratio of the electron 
n = 1.76 X 10" coulomb/kilogram 
Ji, Je d-c. current densities 
U1, Ue d-c. velocities 
Por, poz  _— d-c. charge densities 
Po = —J;/uy » Por = —J2/uz 
PL, pe a-c. charge densities 
V1, Ve a-c. velocities 
Vi, V2 d-c. voltages with respect to the cathode 
V a-c. potential 


Bi = w/m, Bz = w/uy 

Although the small-signal equations relating charge density to voltage V 
have been derived many times, it seems well to present them for the sake 
of completeness. For one stream of electrons the first-order force equation is 


lv v 1 
OF nt Ot 4, Pie, gee 


dt ot Ox Ox 
(w — Buy)o = —nBV 
—7nBV 
- 2 
(81 — B) @) 


From the conserv ation of charge we obtain to the first order 


0 ad 
a — ro (ports + pili) 
wpi = pub + mp 
ia 
(81 = B) 
Jip 
-= — 3 
p= — a2(6, — 8) 
From (2) and (3) we obtain 
nJ, B°V 
= ; - 4 
or ui(Bi — BY 4) 
We would find similarly 
nJ2B?°V 


p>  ux(B2 — 8) (S) 


It will be convenient to call the fractional velocity separation }, so that 


2(u, — ue) 
b= 6 
Uy, + UW (6) 
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It will also be convenient to define a sort of mean velocity uo 


_ 2 Ue 
lll U, + Ue 7) 
We may also let Vo be the potential drop specifying a velocity uo , so that 
uy = V2nV» (8) 
It is further convenient to define a phase constant based on uo 
Bo = = (9) 
Ny 
We see from (6), (7) and (9) that 
B: = Boll — b/2) (10) 
Be = Boll + 5/2) (11) 
We shall treat only a special case, that in which 
Ji i J: - Jo 
a ae. — 


Here Jo is a sort of mean current which, together with mo, specifies the 
ratios J;/ui and J2/us , which appear in (4) and (5). 

In terms of these new quantities, the expression for the total a-c. charge 
density p is, from (4) and (5) and (8) 


Jo8 
20 Vo 


p = pi t+ pe = 


1 1 - 


a(t —2) oT fae) oy] 


Equation (13) is a ballistical equation telling what charge density p is 
produced when the flow is bunched by a voltage V. To solve our problem, 
that is, to solve for the phase constant 8, we must associate (13) with a 
circuit equation which tells us what voltage V the charge density produces. 
We assume that the electron flow takes place in a tube too narrow to propa- 
gate a wave of the frequency considered. Further, we assume that the 
wave velocity is much smaller than the velocity of light. Under these 
circumstances the circuit problem is essentially an electrostatic problem. 
The a-c. voltage will be of the same sign as, and in phase with, the a-c. 
charge density p. In other words, the “circuit effect” is purely capacitive. 

Let us assume at first that the electron stream is very narrow compared 
with the tube through which it flows, so that V may be assumed to be con- 
stant over its cross section. We can easily obtain the relation between 
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V and p in two extreme cases. If the wavelength in the stream, A,, is 
very short (8 large), so that transverse a-c. fields are negligible, then from 
Poisson’s equation we have 


ae 
a= €0 O22 
p=ea8 V (14) 


If, on the other hand, the wavelength is long compared with the tube radius 
(8 small) so that the fields are chiefly transverse, the lines of force running 
from the beam outward to the surrounding tube, we may write 


p=CV (15) 

















ae 


Fig. 1—A “circuit” curve for a narrow electron stream in a tube. The ratio of the a-c. 
charge density p to the a-c. voltage V produced by the charge density is plotted vs. a 
parameter 8/8 , which is inversely proportional to the wavelength A, in the flow. Curve 1 
holds for very large values of 8/8» ; curve 2 holds for very small values of 8/8» , and curve 
3 over-all shows approximately how p/V varies for intermediate values of 8/8o . 

Here C is a constant expressing the capacitance per unit length between the 
region occupied by the electron flow and the tube wall. 

We see from (14) and (15) that if at some particular frequency we plot 
p’V vs. B/Bo for real values of 8, p/V will be constant for small values of 8 
and will rise as 8? for large values of 8, approximately as shown in Fig. 1. 
For another frequency, 8» would be different and, as p/V is a function of 8, 
the horizontal scale of the curve would be different. 

Now, we have assumed that the charge is produced by the action of the 
voltage, according to the ballistical equation (10). This relation is plotted 
in Fig. 2, for a relatively large value of Jo/#oVo (curve 1) and for a smaller 


value of Jo/moVo (curve 2). There are poles at 8/6) = 1+ ae and a mini- 


mum between the poles. The height of the minimum increases as Jo/moVo 
is increased. 
A circuit curve similar to that of Fig. 1 is also plotted on Fig. 2. We see 
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that for the small-current case (curve 2) there are four intersections, giving 
four real values of 8 and hence four unattenuated waves. However, for the 
larger current (curve 1) there are only two intersections and hence two 
unattenuated waves. ‘The two additional values of 8 satisfying both the 
circuit equation and the ballistical equation are complex conjugates, and 
represent waves traveling at the same speed, but with equal positive and 
negative attenuations. 














Q|> 
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Fig. 2—Curve 3 is a circuit curve similar to that of Fig. 1. Curves 1 and 2 are based 
on a ballistical equation telling how much charge density p is produced when the voltage 
V acts to bunch a flow consisting of electrons of two velocities. The abscissa, 8/Bo, 
is proportional to phase constant. Intersections of the circuit curve with a ballistical 
curve represent waves. Curve 2 is for a relatively small current. In this case inter- 
sections occur for four real values of 8, so the four waves are unattenuated. For a larger 
current (curve 1) there are two intersections (two unattenuated waves). For the other 
two waves Biscomplex. There are an increasing and a decreasing wave. 





Thus we deduce that, as the current densities in the electron streams are 
raised, a wave with negative attenuation appears for current densities 
above a certain critical value. 

We can learn a little more about these waves by assuming an approximate 
expression for the circuit curve of Fig. 1. Let us merely assume that over 
the range of interest (near 8/8) = 1) we can use 


p = ae BV (16) 


Here a’ is a factor greater than unity, which merely expresses the fact that 
the charge density corresponding to a given voltage is somewhat greater 
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than if there were field in the « direction only and equation (11) were valid. 
Combining (16) with (13) we obtain 
1 1 1 


(Oye 


Uf ws gees (18) 


2a € B80 Uo Vo 


where 


In solving (17) it is most convenient to represent 8 in terms of §» and a 
new variable 6 


B = Bo(1 + 4) (19) 
Thus, (14) becomes 
1 1 1 


—. + ous - or — 
b\' by ~ (20) 
(0-3) (+5) 


Solving for 6, we obtain 


OH OYE)T & 


The positive sign inside of the brackets always gives a real value of 6 
and hence unattenuated waves. The negative sign inside the brackets 
gives unattenuated waves for small values of U/b. However, when 


yy... 
() ae (22) 


there are two waves with a phase constant 8) and with equal and opposite 
attenuation constants. 
Suppose we let Uy be the minimum value of U for which there is gain. 
From (22), 
Um: = 87/8 (23) 


From (21) we have for the increasing wave 


AW) 


The gain in db/wavelength is 


db/wavelength = 20(27) logio e | 6 | (25) 


54.6 | 6 | 


Il 
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We see that by means of (24) and (25) we can plot db/wavelength per 
unit 6 vs. (U/U yy). This is plotted in Fig. 3. Because U? is proportional 
to current, the variable (U/U 4)* is the ratio of the actual current to the 
current which will just give an increasing wave. If we know this ratio, 
we can obtain the gain in db/wavelength by multiplying the corresponding 
ordinate from Fig. 3 by 0. 


30 


25 


20 


DB GAIN / WAVELENGTH / UNIT D 


3 


F(W/Wy)* 





1 2 4 6°8 10 20. 40 6080100 200 400 600 1000 
2 2 
(W/W )° = (U/Um) 


Fig. 3—The parameter (W/W y)? = ((U'/U y)?is proportional tocurrent. As the current 
is increased above a critical value for which (W/W yy)? = 1, there is an increasing wave of 
increasing gain. In this curve the gain per wavelength per unit }, called F(W/W y)?, is 
plotted vs. (W/W y)*. For large values of (W/W y,)?, F(W/W )? approaches 27.3 and the 
gain per wavelength approaches 27.3 b. 


We see that, as the current is increased, the gain per wavelength at first 
rises rapidly and then rises more slowly, approaching a value of 27.30 
db/wavelength for very large values of (Ul, Uy)’. 

We now have some idea of the variation of gain per wavelength with veloc- 
ity separation 6 and with current (l’/U’'y)*.. A more complete theory would 
require the evaluation of the lower limiting current for gain (or of U*,) in 
terms of physical dimensions and an investigation of the boundary condi- 
tions to show how strong an increasing wave is set up by a given input 
signal. The latter problem will not be considered in this paper; the former 
is dealt with in the third section and in the appendix. 
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3. DESIGN CURVES 

It is proposed to present in this section material for actually evaluating 
the gain of the increasing wave for a particular geometry of electron flow. 
In this section there is some repetition from earlier sections, so that the 
material presented can be used without referring unduly to section 2. In 
order to avoid confusion, much of the mathematical work on which the 
section is based has been put in the appendix. 

The flow considered is one in which electrons of two velocities, #; and ue , 
corresponding to accelerating voltages V; and V2, are intermingled, the 
corresponding current densities J; and J: being constant over the flow. 
The flow occupies a cylindrical space of radius a. It is assumed that the 
surrounding cylindrical conducting tube is so remote as to have negligible 
effect on the a-c. fields. 

It will be assumed, according to (12), that the current densities and the 
voltages V; and V2 are specified in terms of a “mean” current Jo and a 
“mean” voltage Vo corresponding to a velocity um , by 


Jy J Jo , 
yi? = hn = ys? (12a) 


The gain will depend on the beam radius, the free-space wavelength A, 
and on Jy and Vo , and on the fractional velocity separation 


2(u, — Ue) 
b a ( 1 2) (6) 
WM, + UW 
The wavelength in the beam, A, , which is associated with the voltage 
Vy is given by 
Mo V 2nVo 
ee 
c c 
dh, = 1.98 X 10 °A VV, (20) 
Here c is the velocity of light. 
A dimensionless parameter W is defined to be 


9 ‘ Ty 
Ss = 2 = ial = (27) 
w* En Uy W~ 
2 < 6 Jo 
W? = 852 X 10° — 7 (28) 


Here w, is the electron plasma frequency associated with the average space 
charge density Jo/m , and w is the radian frequency corresponding to the 
wavelength A. In (28), the constant is adjusted so that Jo is expressed in 
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amperes per square centimeter rather than in amperes per square meter, 
while f is expressed in megacycles. 

Below a minimum value of W, which will be called Wy , there is no gain. 
Wy is a function of the velocity separation b and of the ratio of the beam 
radius a to the beam wavelength, A,. A plot of (W4,/b)? as a function of 
(a/X,) is shown in Fig. 4. 

The variation of gain in the interval, Wy < Wo< «, is shown in Fig. 3 
where “Decibels gain/wavelength/unit 6” is plotted as a function of 
(W/Ww)?. This is the same curve which was derived in section 2. The 
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Fig. 4—As the ratio of beam radius a to wavelength in the beam, a, , is increased, the 
critical value of W, Wy , decreases and less current is needed in order to obtain gain. 
Here (W y/b)?, which is called H(a/d,), is plotted vs. (a/A,). 















































ratio (W/W)? is the same as the parameter (U/U )? used there, although 
Land W are not the same. 

The curve in Fig. 3 is useful in that it reduces the interdependence of a 
large number of parameters to a single curve. However, there are cases as, 
for example, when one is computing the bandwidth of an amplifier, in which 
it would be more convenient to have the curve in Fig. 3 broken up into a 
family of curves. We can do this by the following means: 

We can write the gain in db/wavelength in the form 


db/wavelength = b/(W/W 4)? (29) 
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Here F(W/W )? is the function plotted in Fig. 3. If ¢ is the total length 
of the flow, the total gain in db, G, will thus be 


G = = F(W/Wau) (30) 


We will now express (W/W )* in such a form as to indicate its dependence 
on wavelength in the beam, A,. We can write from (27) 
Wm MB (31) 
w? AZ 
Here \, is a “plasma wavelength,” defined by the relation 


Uo 


at (we/2m) 


(32) 


We further have 
Wy = BH(a/d.) (33) 


Here H(a/X,) is the function of (a/X,) which is plotted in Fig. 5. 
30 
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Fig. 5—In these curves the total gain in db, G, divided by the ratio of the length ¢ to 
the plasma wavelength X,, is plotted vs. bw/w,, which is proportional to frequency, for 
several values of the parameter (a/bA,)?. Changing 6, the velocity seperation, changes 
both the parameter and the frequency scale. 


Now, from (26), (27), and (29) we can write 


a Dh t bre a > Ne ' a ’ 1 
nei ({) (™) (): I(*) (<) wana (34) 


For a given tube the parameters (¢/A,) and (a/bd,) do not vary with fre- 
quency, while (a/A,) is proportional to frequency. Hence, we can construct 
universal frequency curves by plotting G/(¢/X,) vs. (a/X;) for various values 
of the parameter (a/bd,). It is more convenient, however, to use as an 
abscissa bA./A; = bw/we , and this has been done in Fig. 5. 
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In order to use these curves it is necessary to express the parameters 
bw /w, , and (a/bd,)? in terms of convenient physical quantities. We obtain 


b/w. = 545 X 107° OV" w/ Ih” 
Ae = 2.04 & 10° 54/9” (35) 
(a/bd.)? = 767 Io/b?V*? 


Here / is current in amperes and J, is in amperes / cm.’ 

The broadness of the frequency response curves of Fig. 5 is comparable 
to that of curves for helix-type traveling-wave tubes. 

It is interesting to note that the maximum value of G/(¢/X,) varies little 
for a considerable range of the parameter a/bd, , approaching a constant 
for large values of the parameter. This means that, with a beam of given 
length, velocity and charge density, one can obtain almost the same opti- 
mum gain over a wide range of frequencies simply by adjusting the velocity- 
separation parameter 6. 


4. CONCLUDING REMARKS 


There is a great deal of room for extension of the theory of double-stream 
amplifiers. This paper has not dealt with the setting up of the increasing 
wave, nor with other geometries than that of a cylindrical beam in a very 
remote tube, nor with the effect of physical separation of the electron streams 
of two velocities nor with streams of many velocities or streams with con- 
tinuous velocity distributions. 

This last is an interesting subject in that it may provide a means for deal- 
ing with problems of noise in multivelocity electron streams. Indeed, it 
was while attempting such a treatment that the writers were distracted by 
the idea of double-stream amplification. 


APPENDIX 
DERIVATION OF REsuLTS USED IN SECTION 3 


Consider a double-stream electron beam whose axis coincides with the 
z-axis of a system of cylindrical coordinates (r, y, z) and which is subject to 
an infinite, longitudinal, d-c. magnetic field. The radius of the beam is a 
and each of the streams is characterized by d-c. velocities, #; and “2 , which 
are vectors in the positive z direction, and d-c. space charge densities, poz 
and po2. All d-c. quantities are assumed to be independent of the coor- 
dinates and time, except, of course, for the discontinuities at the surface of 
the beam. Small a-c. disturbances are superimposed upon these d-c. quanti- 
ties and they are small enough so that their cross products can be neglected 
compared with the products of d-c. quantities and a-c. quantities. It is 
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further assumed that only those a-c. quantities are allowed which have no 


: aise = 2 se 
azimuthal variation, that is, gg 0. Fig. 6 shows the electron beam. 
(0) 


Outside the beam the appropriate Maxwell’s equations are 


1 0a k 
foie er. 5 A-1 
-- (rH) j = E, ( ) 
aH, a. 
a =—j-E A-2 
02 J No ( ) 
OE, dE ; ‘ 
2 r = A-3 
ay ae ikno He (A-3) 
where 
k=* (A-4) 
tf 
/ uo aw P 
no = NV “ = 377 ohms (A-5) 
r & 


Inside the beam, equations (A-2) and (A-3) remain the same, but instead 
of equation (A-1) we have 


1 0 ok 
1a ine . A-6 
2 (rH,) = j = E,.tat@ (A-6) 


where g; and q are the first order a-c. convection current densities of the two 
streams. These quantities can be calculated from the force equation and 
the equation for the conservation of charge. Assuming that all a-c. quanti- 
ties vary as exp j(wi — Bz), the force equation is (for stream number one, say) 


jor, — j8uyr, = —(e/m)E, (A-7) 
and the equation for the conservation of charge is 
JBporrs + 78urp1 _ +jwp1 (A-8) 
Equations (A-7) and (A-8) can be solved for 2; and p; : 
a —(e/m)E, 
aa cani (: *) (A-7a) 
m we 
J B; 
ae 
a= ( 4 1 (A-8a) 
w{1— 
1 
where : 
Bi, = = 
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Combining equations (A-7a) and (A-8a) one has 


_ Booule/m)E, 
ri = ( 4) (A-9) 
(a) — — 
Bi 
The first order a-c. convection current density is given by 
gi = ports + pitty (A-10) 


which, by combining with (A-7a) and (A-8b), becomes 


J(R/n0)( po1/meo)E, 


| aa 2 (1 rs £y (A-11) 
By 
Similarly 
j k € 
| 
n me 
at a (A-12) 
w (1 = iy 
Bo 
If we now define 
Bo = 3(81 + Ba) (A-13) 
By Be 
B= 7; B= A-14 
aus” Bo (As) 
and let 
‘ B 
Z= -15 
Bo (A-15) 
Wm="; w=. (A-16) 
é3] Ww 


where w, , the plasma-electron angular frequency given by 


2 €po 
ws = —— , ete. (A-17) 
Meo 


Equations (11) and (12) become 
—j(k/no)Wi Bi 


‘a ve 
_ —i(k/no)W2 Bz 
QQ = (Z = B,)? (A-19) 
Thus equation (A-6) becomes 
1a >. oe 
Pa ar (rH,) —_ J no Ci (A-6a) 
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where 
2 2 72 2 
Wy By Ws Bs 


ae ge ae 
1 — GB): ~ (Z— Bs)? 


(A-20) 


If we assume that the tube which surrounds the beam be taken as infinitely 
remote, the appropriate solutions outside the beam are 


Hoo = Ao Ki(yr) (A-21) 
Ba =j ~ Ao Ko(yr) (A-22) 
and inside the beam 
A, = A; h(ér) (A-23) 
En = -j VE A; Io(ér) (A-24) 
2 =— A2 — 2 ~ 2 > 
where y B k B } (A-25) 
P= yL 


The J’s and K’s in equations (A-21)-(A-24) are modified Bessel functions.’ 
At the surface of the beam (r = a), one has the following two independent 
boundary conditions 


By = Ay (A-26) 
Bus = Ba (A-25a) 
which, using equations (A-21)-(A-24), yield 
To(&a) K,(ya) » 
= —_ (A-27 
VL I(ée) Ki(ya) 
From equations (A-13), (A-14), (A-15) and (A-24) one has 
ta = ZByavW/L (A-28) 
ya = ZBoa (A-29) 
If we now define a beam wavelength, A, , by the relations 
2 
Bo = . (A-30) 


and assume for the purpose of simplifying the calculation that in the ex- 
pression for Z in (A-20) 


WiBi = WB; = W? (A-31) 


* See A Treatise on the Theory of Bessel Functions, G. N. Watson, Chapter 3. 
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We easily see that 
W? = (w/w) (A-32) 
e 
where we = —Jo/eouo (A-33) 
m 


We obtain from (A-20), (A-28), (A-29) and (A-30) 
Neon 3 / . Dad 2 
| K: (722) 1 (VE arsz) 
Xs As ) 
| ka (2292) 1, (vz 7%%”) (A-34) 


Ww" Ww 
vat eed 5 a a | 
E - By" (Z — By)? 
Equation (A-31) is equivalent to Equation (12) of the text or to 
Ji Js 


yi? a yi? (A-35) 











Fig. 6—The diameter of the electron flow considered is 2a, and the length is ¢. 


Letting ¥ = —Z and making use of the following well known relations 
between the Bessel functions 


Io(jx) => J o(x) 


(A-36) 
I,(jx) = jJi(x) 
Equation (A-34) becomes 
2 ae . 

- (20 z) r bee mu :) 
. Ng Ns 
ian 2 2 

Ko( — 7) J (vr 7) (A-37) 


NW” Ww 
~ @—B, ” zF— By 
Let the right-hand side of equation (A-37) be denoted by F(Z) and the 
middle of F2(Z). In order to find the real roots of equation (A-37) one can 
plot F; and F»2 as functions of Z on the same chart. The abscissae of the 
intersections of the two curves will then be the real roots. In Fig. 7, Fi 
is plotted as a function of Z for B, = 0.9 and By = 1.1. 
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In view of the definitions in equations (A-13) and (A-14), both B, and B, 


are uniquely defined by a single parameter, namely, the fractional velocity 
separation, b. That is 


b = 2(u; — u2)/(ui + U2) 


2(B2 st B1)/(B2 + B1) 


(A-38) 
By —_ B, 


By = 1— (6/2) 
Bz = 1+ (6/2) 


and 


(A-39) 
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Fig. 7—A curve illustrating conditions giving rise to various types of roots. 


A complete plot of #2, for any value of the parameters W and (a@/A,), 
would show that equation (A-37) has an infinite number of real solutions. 
A real solution of equation (A-37) means an unattenuated wave. Thus 
there are an infinite number of unattenuated waves possible. The waves 
which will actually be excited in any given case, however, depend upon the 
boundary conditions at the input and output of the tube. Ordinarily only 
those waves will be excited which do not have a reversal in phase of the 
longitudinal E vector, say, as r varies from 0 to a. Attention, therefore, 
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will be given only to those waves for which EF, does not change sign over a 
cross-section of the beam. By inspection of equation (A-23), it is evident 
that this requirement is automatically satisfied if L > 0. On the other 
hand, if LZ is negative, one has 


Eu, ~ Jo (v Yy om z) (A-23a) 
Thus attention will be limited to those roots which satisfy 


VY 7 7 < 2.408 (A-40) 
where 2.405 is the first zero of the Bessel function in equation (A-21a). 

Returning to Fig. 7, portions of three different /2 curves are plotted: one 
for W? = 0.01, one for W? = 0.0152 and one for W? = 0.02. All three 
curves are for (a4/X,) = 0.16. The intersections which represent roots which 
satisfy the inequality (A-40) are marked with arrows. Evidently there are 
either four real roots of this type or there are two real roots and a complex 
conjugate pair, the distinction being determined by the value of W. Thus 
there is a critical value of W? (in this case it is 0.0152) for which two of the 
real roots are identical. This identical pair is indicated by two arrows 
near the minimum of the F; curve at Z = 1. 

A pair of conjugate complex roots means that there are an increasing wave 
and a decreasing wave. Thus for each value of b and (a/X,) there is a least 
value of W? below which the tube will have no gain. 

It can be shown that the critical tangency of the F, and F; curves occurs 
at a value of Z which is less than 5? away from unity. Very little error will 
be incurred, then, by assuming that this critical point occurs at Z = 1 
if 6 is small. 

Letting Z = 1 in equation (A-37), and using equations (A-39) one has 


a Ki(2ma/d.)Jo (0/8(Wy/b)? — 1 are} ) 
sft) — 1. «fo tT Bw edt Aw Sa Us , 
sia Geen (\/8(W y/b)? — 1 2ma/d.) mn 


where W y is the critical value of W. Equation (A-41) determines (W y/6)? 
as af unction of (a@/A,). This relationship is plotted in Fig. 4. 

We will find that there will be an increasing wave in the range 
Wu < W< «. The calculation of the gain in this interval would be very 
laborious since Bessel functions of complex argument would be involved. 
However, a good approximation can be made when 6 is small. The real 
part of Z will always be near unity and the imaginary part will be found to 
be less than 6/2. Therefore one can let Z = 1 in equation (A-37) where it 
multiplies the factor (27a/A,) in the argument of the Bessel functions and 
let Z — 1 = U in the right-hand side of Equation (A-37). With these 
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assumptions Y can be determined as a function of (a/A,) and U can be deter- 
mined asa function of Y. We have from Equation (A-37) 


1 1 _1+Y 


+t Os Ww aati: 
When U = 0, W’ = Wy, = W’y , so that 
1+ Y = 8(Wy/b)? (A-42) 
and equation (A-37a) becomes 
1 1 Di tal " 
(U + 6/2)? ao (u — 4/2)" = (8/b)(Wy/W) (A-37b) 


the solution of which, for the increasing wave, is 
U = j(b/2){(1/2)(W/Wau) (1 + 8(Wy/W)? — 1) — 1]' (4-43) 
and the gain will be given by 
Gain/b = 27.3[(1/2)(W/Wau)?(V1 + 8(Wy/W)? — 1) — 1]'  (A-44) 
db/wavelength/unit 5 


“Decibels gain/wavelength/unit b” is plotted against (W/W)? in Fig. 3. 
As (W/W .)* becomes very large, the gain per wavelength approaches 
27.3 6b db. 

















Experimental Observation of Amplification by Interaction 
Between Two Electron Streams 


By A. V. HOLLENBERG 


The construction and performance of an amplifier employing the interaction 
between two streams of electrons having different average velocities are de- 
scribed. Gain of 33 db at a center frequency of 255 Mc has been observed 
with bandwidth of 110 Mc between 3 db points. 


1. INTRODUCTION 


NEW type of amplifier in which the gain is obtained by an interaction 
between streams of electrons of two or more average velocities is 
proposed in a companion paper by Pierce and Hebenstreit.!. This amplifier 
contains input and output portions in which signals are impressed on and 
extracted from the electron flow by electromagnetic circuits and a central 
portion in which gain occurs purely by interaction between streams of elec- 
trons without any circuits being involved. A small signal theory for coinci- 
dent electron streams of two d-c. velocities is presented in Pierce and Heben- 
streit’s paper. 

In this paper a description of the construction and operation of an amplifier 
of this kind will be presented. Departures of the actual conditions in the 
amplifier from the assumptions of the theory limit the expectations of quanti- 
tative agreement. It is believed, however, that the evidence for gain arising 
from the interaction between two streams of electrons is clear, and that the 
broad frequency response predicted by the theory has been confirmed. 


2. DESCRIPTION OF AMPLIFIER 


The frequency range near 200 Mc was chosen for the first experimental 
test of the proposed method of amplification for reasons of convenience. 
The theory indicates that current density requirements increase with fre- 
quency, but that these requirements become severe only at the higher micro- 
wave frequencies. Availability of circuit parts and test equipment, rather 
than anticipated difficulties at higher frequencies, led to the choice that was 
made. 

The essential features of one of the double-stream amplifier tubes which 
has been constructed and operated are represented in Fig. 1. The out- 
put helix was identical with the input helix in construction and connection 

1A New Type of High Frequency Amplifier, J. R. Pierce and W. B. Hebenstreit, this 
issue of the B. S. T. J. 
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to the coaxial line. The two identical probes p,, and p2 , extending from 
coaxial lines into the two electron streams at the beginning and end of the 
central portion of the tube between the two helices were inserted for com- 
parison of the signal amplitudes at the beginning and end of the region in 
which no circuit is present. 

A similar tube containing an output gap in place of the output helix sec- 
tion is represented in Fig. 2. 

In both cases concentric tubular electron streams originate at the ring- 
shaped emitting surfaces of the two cathodes at potentials V; and V2, pass 
through their respective control grids and through a common accelerating 
grid. An axial magnetic field of approximately 700 gauss is applied in order 
to maintain the definition of the beams. The outer and inner tubular beams 
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Fig. 2—Representation of double-stream amplifier with gap output. 


have mean radii of 0.215” and 0.170” respectively and a wall thickness of 
0.030”. in each case. 

The short sections of helix which are used for input and output are wound 
of 0.013” diameter molybdenum wire, 44 t.p.i., and mean diameter of 0.500”. 
The axial velocity of signal propagation along this helix is equal to that of 
54-volt electrons. The helix sections are each 2” long. Ceramic supporting 
rods on each helix section are sprayed with aquadag, over 14” of their 
length on the end nearest the center of the tube, for terminating purposes. 
The thickness of the spray coating increases toward the center of the tube. 
The distance between helices is 8.7”. 

The gain produced by the electronic interaction depends upon a difference 
in velocity between the two electron streams. The signal is impressed 
upon one of the streams by the helix when its velocity is that at which travel- 
ing wave amplifier interaction between the stream and the helix occurs. It 
































INTERACTION BETWEEN ELECTRON STREAMS 


is required, therefore, for this helix that one of the streams travel at a 
velocity corresponding to a potential in the neighborhood of 54 volts. 
Useful interaction occurs from 50 to 60 volts. The inner stream is adjusted 
for helix interaction in this amplifier, and the outer stream travels at a lower 
velocity to bring about the interaction between the two streams. Operation 
about a mean voltage of about 50 volts was planned in designing the ampli- 
fier, and in estimating its expected performance. The amplifier is 16 wave- 
lengths long in terms of the wavelength associated with a mean voltage of 
50 volts and a frequency of 200 Mc. Eleven of these wavelengths are in the 
center portion between the helices. 

The conditions in the amplifier tube differ from those assumed in the 
derivation of the theory of the double-stream interaction in the following 
significant ways: 

1. The beams are separated in space and not completely intermingled. 
Calculations on the effect of this separation have been made. Numeri- 
cal examples of the calculated magnitude of the effect on gain will be 
given below. 

2. Hollow tubular beams are used, instead of “‘solid” beams of uniform 
current density over their cross-sectional area. The theory indicates 
that, for the beam dimensions and currents used here, the parameters 
which depend upon beam radius and total current in the beam are 
nearly the same whether the current is concentrated in an infinitely 
thin cylindrical shell or uniformly distributed over the cross-section of a 
cylinder of the same radius. 

3. The metal wall surrounding the beams is not infinitely remote. Its 
diameter was chosen as a compromise between the requirements of 
preventing serious d-c. space charge depression of potential in the 
beam and of being far enough removed from the beam to prevent a large 
effect on the interaction due to its presence. Its proximity tends to 
increase the minimum current required for producing gain, and there- 
fore to reduce the ratio of actual to critical current on which the gain 
depends. 

4, The beams are not perfectly confined to hollow cylinders of the dimen- 
sions given. There is evidence that some spreading outside of these 
dimensions occurs. The currents reaching the collector can be meas- 
ured and these are used as “‘beam currents” in the discussion to follow 
and in comparisons between theory and experiment. Somewhat larger 
currents than these were initially launched, and the lost fraction may 
have contributed to the interaction before striking the walls. 

Although the assumptions of the theory are not fulfilled in the actual 
amplifier, estimates of its performance were first made without correction 
for the discrepancies. With voltages of 40 and 60 volts on the outer and 
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inner streams and currents of 0.5 and 0.8 milliampere, a gain of 40 db at 200 
Mc due to the double-stream amplification was predicted, with bandwidth 
to 3 db points of 90 Mc, centered about 200 Mc. A later estimate, in- 
cluding the effect of the separation of the hollow beams in space, reduced 
the 40 db figure to 23 db. 

Prediction of the performance of the device as a whole also requires an 
evaluation of the coupling of the helix sections to the electron stream. The 
length of the active portions of the helix sections was chosen to give gain of 
order unity as estimated from single-stream traveling wave amplifier con- 
siderations for the proposed operating current. 
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POWER INPUT 
DECIBELS REFERRED TO 1 MILLIWATT 


big. 3—Power output versus power input for double-stream amplifier at 200 Mc with 
heam potentials of 54 volts and 33 volts and 1.1 ma current in each beam. 


3. EXPERIMENTAL OBSERVATIONS 


The two beam potentials and the two beam currents were varied until 
good operating conditions were reached. In the amplifier tube with helix 
output, good operation with 29 db gain at 200 Me was observed at low input 
signal levels with inner beam potential 54 volts, outer beam potential 33 
volts, and 1.1 milliampere current in each beam. Measurements of output 
and input signal power are shown in Fig. 3, in which output power is 
plotted as a function of input power. The gain is seen to be constant at 
29 db from low levels up to an output power of 0.03 milliwatts, at which 
point compression sets in. Maximum power output is 0.3 milliwatts. 
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The variation of gain with frequency under very nearly the same condi- 
tions as above is shown in Fig. 4. The input signal was in the linear region 
of Fig. 3. The figure shows a bandwidth of 110 Mc between points 3 db 
down from maximum. The center of the band is at approximately 255 Mc. 

The currents required to realize the above results are some tenths of a 
milliampere higher than those used in making estimates of performance of 
the amplifier, and the voltages for best performance are lower. The use 
of the actual operating values leads to a prediction of maximum gain at 
240 Mc for the double-stream interaction in better agreement with the 


110 MC 
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FREQUENCY — MEGACYCLES 

lig. 4—Gain versus frequency for double-stream amplifier, operating in linear portion 
of Fig. 3. 
observed maximum at 255 Mc than the value of 200 Mc originally estimated 
from the lower currents and higher voltages. Estimated double-stream 
gain at 240 Mc from theory for the separated hollow beams is 40 db. 

Evidence that the gain of the amplifier resulted from the double-stream 
interaction was obtained in the following ways: by comparison of the trans- 
mission through the amplifier with both streams with that for each stream 
alone, and by comparison of signals obtained from the two identical probes 
in the electron streams at the ends of the central portion of the tube shown 
in Fig 1. 

In a typical comparison of the first type, 1 db gain for the device was 
observed with the inner stream alone. This increased to 29 db when the 
outer stream was turned on. With the outer beam alone the loss in the 
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device was very large, for the velocity of the outer beam was far from that 
at which interaction with the helices occurs. 

A signal from the probe at the end of the central portion of the tube 
23 db greater than that from the probe at the beginning of this section was 
observed in a comparison of the second type. This can probably be taken 
as a measure of the increase in signal in this portion of the tube due to the 
double-stream interaction alone, although the probe arrangement may also 
be subject to some remaining complicating effects. Overall gain for the 
device in this measurement was 32 db. Further interaction of the same 
kind occurs in the portions of the tube outside of the space between the 
probes. 

Measurements of the gain of an amplifier with helix output.as a function 
of velocity separation between the streams have been made. For fixed 
mean voltage and current, theory predicts an increase in gain from zero db 
at zero separation to a maximum and then a decrease to zero as the velocity 
separation is further increased. A maximum gain was observed experi- 
mentally as velocity separation was varied, and in the neighborhood of the 
predicted optimum value of velocity separation for the current used. 

In the amplifier tube with gap output it was possible to evaluate the a-c. 
component of current in the electron stream produced by the amplified 
signal since the impedance across the gap was known. The power output 
from this tube at saturation was 0.1 mw, a little less than the maximum 
shown in Fig. 3. For 75 ohm output impedance this power corresponds 
to 1.15 milliampere r.m.s., or about one third of the total d-c. current to the 
collector in both streams. The output power, although relatively low, is 
thus of the right order of magnitude for the currents used. 
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The Synthesis of Two-Terminal Switching Circuits 
By CLAUDE. E. SHANNON 
PART I: GENERAL THEORY 
1. INTRODUCTION 

HE theory of switching circuits may be divided into two major divi- 

sions, analysis and synthesis. The problem of analysis, determining 
the manner of operation of a given switching circuit, is comparatively 
simple. ‘The inverse problem of finding a circuit satisfying certain given 
operating conditions, and in particular the besé circuit is, in general, more 
difficult and more important from the practical standpoint. A basic part 
of the general synthesis problem is the design of a two-terminal network 
with given operating characteristics, and we shall consider some aspects of 
this problem. 

Switching circuits can be studied by means of Boolean Algebra.':? This 
is a branch of mathematics that was first investigated by George Boole in 
connection with the study of logic, and has since been applied in various 
other fields, such as an axiomatic formulation of Biology,’ the study of neural 
networks in the nervous system,‘ the analysis of insurance policies,® prob- 
ability and set theory, etc. 

Perhaps the simplest interpretation of Boolean Algebra and the one 
closest to the application to switching circuits is in terms of propositions. 
A letter X, say, in the algebra corresponds to a logical proposition. The 
sum of two letters Y + Y represents the proposition “X or Y” and the 
product XV represents the proposition “X and Y”. The symbol X’ is used 
to represent the negation of proposition X, i.e. the proposition “not X”’. 
The constants 1 and 0 represent truth and falsity respectively. Thus 
X -+ Y = 1 means X or Y¥ is true, while Y + YZ’ = 0 means X or (Y and 
the contradiction of Z) is false. 

The interpretation of Boolean Algebra in terms of switching circuits®:*:**'° 
is very similar. The symbol X in the algebra is interpreted to mean a make 
(front) contact on a relay or switch. The negation of X, written X’, 
represents a break (back) contact on the relay or switch. The constants 0 
and 1 represent closed and open circuits respectively and the combining 
operations of addition and multiplication correspond to series and parallel 
connections of the switching elements involved. These conventions are 
shown in Fig. 1. With this identification it is possible to write an algebraic 
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Fig. 1—Hindrance functions for simple circuits. 


expression corresponding to a two-terminal network. This expression will 
involve the various relays whose contacts appear in the network and will be 
called the hindrance or hindrance function of the network. The last net- 
work in Fig. 1 is a simple example. 

Boolean expressions can be manipulated in a manner very similar to 
ordinary algebraic expressions. Terms can be rearranged, multiplied out, 
factored and combined according to all the standard rules of numerical 
algebra. We have, for example, in Boolean Algebra the following identities: 


0+X=X 
0-X = 0 
1-X=X 
X+VY=Y+X 
XY = YX 


X¥+(¥ +2) =(X¥+Y)+2Z 
X(YZ) = (XY)Z 
X(¥+Z) = XY+ XZ 
The interpretation of some of these in terms of switching circuits is shown 
in Fig. 2. 
There are a number of further rules in Boolean Algebra which allow 
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Fig. 2—Interpretation of some algebraic identities. 


simplifications of expressions that are not possible in ordinary algebra. 


The more important of these are: 


X=X+X 
Ae X-X 
X +1 


X+YZ= 
X+X'= 
X-X' = 
(X+ Y)’ = 
(XY)! = 


The circuit interpretation of some of these is shown in Fig. 3. 


X+X+4X=etc. 
X-X-X = ete. 

1 

(X + Y)(X + Z) 

1 

0 

rr 

X’+Y’ 


These rules 


make the manipulation of Boolean expressions considerably simpler than 
ordinary algebra. ‘There is no need, for example, for numerical! coefficients 
or for exponents, since nX = X" = X. 

By means of Boolean Algebra it is possible to find many circuits equivalent 
in operating characteristics to a given circuit. The hindrance of the given 
circuit is written down and manipulated according to the rules. Each 
different resulting expression represents a new circuit equivalent to the given 
one. In particular, expressions may be manipulated to eliminate elements 
which are unnecessary, resulting in simple circuits. 


Any expression involving a number of variables X;, X2, «++, X, is 
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called a function of these variables and written in ordinary function notation, 
f(X1,X2,--+,X,). Thus we might have f(X, Y,Z) = X + Y'Z+ XZ’. 
In Boolean Algebra there are a number of important general theorems which 
hold for any function. It is possible to expand a function about one or more 
of its arguments as follows: 


A(X, Xe, amass » Xn) = Xi f(1, Xe, ahah » Xn) + X’f(0, X2, aite » Xn) 


This is an expansion about X,. The term f(1, X2, --+ , Xn) is the function 
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Fig. 3—Interpretation of some special Boolean identities. 
f(X%1, Xa, +++ , Xn) with 1 substituted for X, and 0 for X’, and conversely 
for the term f(0, X2,--- ,X,). Anexpansion about X,and Xgis: 
f(X1, Xo, +++ Xn) = XiX0f(1, 1, X3, °°, Xn) + X1X5f(1, 0, X3,°++,Xn) 
+ XiX2f(0, 1, Xs, -+:, Xn) + XiX2f(1, 1, Xs, -°- , Xn) 


This may be continued to give expansions about any number of variables. 
When carried out for all variables, f is written as a sum of 2” products 
each with a coefficient which does not depend on any of the variables. 
Each coefficient is therefore a constant, either 0 or 1. 


There is a similar expansion whereby f is expanded as a product: 
W(X1, X2,°++ , Xs) 
= [Xi + fO, X2,---, Xs] (Xi + f(A, Xe, +++, X,)I 
= [(X.+ X2+/(0,0, --- ,X,)] [Xi t X2+ f(0,1, ---,Xn)] 
[Xi-+ X2+ f(1, 0, +++, Xa)I(Xi + X2 + f(, 1, +++. Xn] 


= etc. 
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The following are some further identities for general functions: 
X + f(X, Y, Z,---) = X + f(0, Y, Z, ---) 
X' + f(X, Y,Z = X’+ f(1, Y, Z, ---) 
aie, ¥,.2,°« Xf(1, Y, Z, ---) 
X'f(0, Y, Z, -+-) 
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= xy Ff(11,¥Z) + X Y’#(1,0,¥%,Z2) + XY £(0,1,¥,Z) + xX’ Y’ £(0,0, Y,2) 
Fig. 4—Examples of some functional identities. 


The network interpretations of some of these identities are shown in Fig. 

4. A little thought will show that they are true, in general, for switching 
circuits. 
; The hindrance function associated with a two-terminal network describes 
the network completely from the external point of view. We can determine 
from it whether the circuit will be open or closed for any particular position 
oftherelays. This is done by giving the variables corresponding to operated 
relays the value 0 (since the make contacts of these are then closed and the 
break contacts open) and unoperated relays the value 1. For example, with 
the function f = W[X + Y(Z + X’)] suppose relays X and Y operated and 
Z and W not operated. Then /f = 1[0+ 0(1 + 1)] = Oand in this condition 
the circuit is closed. 
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A hindrance function corresponds explicitly to a series-parallel type of 
circuit, i.e. a circuit containing only series and parallel connections. This 
is because the expression is made up of sum and product operations. There 
is however, a hindrance function representing the operating characteristics 
(conditions for open or closed circuits between the two terminals) for any 
network, series-parallel or not. The hindrance for non-series-parallel net- 
works can be found by several methods of which one is indicated in Fig. 5 
for a simple bridge circuit. The hindrance is written as the product of a 
set of factors. Each factor is the series hindrance of a possible path between 
the two terminals. Further details concerning the Boolean method for 
switching circuits may be found in the references cited above. 

This paper is concerned with the problem of synthesizing a two-terminal 
circuit which represents a given hindrance function f(X1, +++ ,X,). Since 
any given function f can be realized in an unlimited number of different 





f = (W+x)(Z+S)(wW+v+8)(Z4+¥+4+x) 


Fig. 5—Hindrance of a bridge circuit. 


ways, the particular design chosen must depend upon other considerations. 
The most common of these determining criteria is that of economy of ele- 
ments, which may be of several types, for example: 

(1) We may wish to realize our function with the least total number of 
switching elements, regardless of which variables they represent. 

(2) We may wish to find the circuit using the least total number of relay 
springs. This requirement sometimes leads to a solution different 
from (1), since contiguous make and break elements may be combined 
into transfer elements so that circuits which tend to group make and 
break contacts on the same relay into pairs will be advantageous 
for (2) but not necessarily for (1). 

(3) We may wish to distribute the spring loading on all the relays or on 
some subset of the relays as evenly as possible. Thus, we might try 
to find the circuit in which the most heavily loaded relay was as 
lightly loaded as possible. More generally, we might desire a circuit 
in which the loading on the relays is of some specified sort, or as near 
as possible to this given distribution. For example, if the relay X, 
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must operate very quickly, while X_ and Xx; have no essential time 
limitations but are ordinary U-type relays, and X, is a multicontact 
relay on which many contacts are available, we would probably try 
to design a circuit for f(X1, X2, Xs, X4) in such a way as, first of all, 
to minimize the loading on X, , next to equalize the loading on X» 
and X; keeping it at the same time as low as possible, and finally 
not to load X4 any more than necessary. Problems of this sort may 
be called problems in spring-load distribution. 

Although all equivalent circuits representing a given function f which 
contain only series and parallel connections can be found with the aid of 
Boolean Algebra, the most economical circuit in any of the above senses will 
often not be of this type. The problem of synthesizing non-series-parallel 

circuits is exceedingly difficult. It is even more difficult to show that a 
' circuit found in some way is the most economical one to realize a given 
function. The difficulty springs from the large number of essentially 
different networks available and more particularly from the lack of a 
simple mathematical idiom for representing these circuits. 

We will describe a new design method whereby any function /(X,, Xo,---, 
X,) may be realized, and frequently with a considerable saving of elements 
over other methods, particularly when the number of variables n is large. 
The circuits obtained by this method will not, in general, be of the series- 
parallel type, and, in fact, they will usually not even be planar. This 
method is of interest theoretically as well as for practical design purposes, 
for it allows us to set new upper limits for certain numerical functions asso- 
ciated with relay circuits. Let us make the following definitions: 

\(n) is defined as the least number such that any function of » variables 
can be realized with not more than A(m) elements.* Thus, any function of 
n variables can be realized with \(m) elements and at least one function with 
no less. 

u(n) is defined as the least number such that given any function f of n 
variables, there is a two-terminal network having the hindrance f and using 
not more than u(7) elements on the most heavily loaded relay. 

The first part of this paper deals with the general design method and the 
behaviour of A(m). The second part is concerned with the possibility of 
various types of spring load distribution, and in the third part we will study 
certain classes of functions that are especially easy to synthesize, and give 
some miscellaneous theorems on switching networks and functions. 

2. FUNDAMENTAL DeEsIGN THEOREM 


The method of design referred to above is based on a simple theorem deal- 
ing with the interconnection of two switching networks. We shall first 


* An element means a make or break contact on one relay. A transfer element means 
a make-and-break with a common spring, and contains two elements. 
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state and prove this theorem. Suppose that M and N (Fig. 6) are two 
(n + 1) terminal networks, M having the hindrance functions U, (k = 
1, 2,--- m) between terminals a and k&, and N having the functions V; 
between band &. Further, let M be such that Uj, = 1(j,k = 1,2, --- , mn). 
We will say, in this case, that M is a disjunctive network. Under these con- 
ditions we shall prove the following: 

Theorem 1: If the corresponding terminals 1, 2,--- , n of M and N are 
connected together, then 


























n 
Us = II (Ux, + Vi) (1) 
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where Up is the hindrance from terminal a to terminal b. 
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Fig. 6—Network for general design theorem. 


Proof: It is known that the hindrance U4, may be found by taking the 
product of the hindrances of all possible paths from a to 6 along the elements 
of the network. We may divide these paths into those which cross the line 
L once, those which cross it three times, those which cross it five times, etc. 
Let the product of the hindrances in the first class be W,, in the second 
class W3;, etc. Thus 


Ua = W,-Ws3-W; --> (2) 
Now clearly 

W, = I] (Ur + Vi) 
and also 

W;=W,=-:-=1 


since each term in any of these must contain a summand of the type Uj, 
which we have assumed to be 1. Substituting in (2) we have the desired 
result. 

The method of using this theorem to synthesize networks may be roughly 
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described as follows: The function to be realized is written in the form of a 
product of the type (1) in such a way that the functions U; are the same fora 
large class of functions, the V;, determining the particular one under consider- 
ation. A basic disjunctive network M is constructed having the functions 
U;. between terminals a and k&. A network N for obtaining the functions 
V; is then found by inspection or according to certain general rules. We 
will now consider just how this can be done in various cases. 


3. DESIGN OF NETWORKS FOR GENERAL FUNCTIONS—BEHAVIOR OF A(n)- 
a. Functions of One, Two and Three Variables: 


Functions of one or two variables may be dismissed easily since the 
number of such functions is so small. Thus, with one variable X, the 
possible functions are only: 


ae Oe 

and obviously A(1) = 1, w(1) = 1. 

With two variables X and Y there are 16 possible functions: 
OX Y XY XY’ X'Y X'Y’ XY'+X’Y 
1X'Y’ X+V X+ VY X+Y X+Y’ XV4+X'yY’ 
so that A(2) = 4, u(2) = 2. 

We will next show that any function of three variables f(X, Y, Z) can be 
realized with not more than eight elements and with not more than four 


from any one relay. Any function of three variables can be expanded in a 
product as follows: 


f(X, ¥,Z) = (X + V +, 0, Z)ILX + ¥’ + f@, 1, 2)] 
[X’ + ¥ + (1,0, ZX’ + ¥' + fA, 1, Z)]. 


In the terminology of Theorem 1 we let 


U,=X+Y Vi = f(0, 0, Z) 
U,=X+Y' V2 = f(0, 1, Z) 
U;=X'+Y V3 = f(1, 0, Z) 
U, = X'+ Y' Vs = f(l, 1, Z) 


so that 
4 
Us = f(X, Y,Z) = I] (Us + Vi) 


The above U, functions are realized with the network M of Fig. 7 and it is 
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easily seen that Uj, = 1(j,# = 1,2,3,4). The problem now is to construct 
a second network N having the V, functions Vi, V2, Vs, V4. Each of 
these is a function of the one variable Z and must, therefore, be one of the 
four possible functions of one variable: 


0,1, Z, 2’. 


Consider the network JN of Fig. 8. If any of the V’’s are equal to 0, connect 
the corresponding terminals of M to the terminal of NV marked 0; if any are 
equal to Z, connect these terminals of M to the terminal of VN marked Z, 
etc. Those which are 1 are, of course, not connected to anything. It is 
clear from Theorem 1 that the network thus obtained will realize the function 
f(X, Y, Z). In many cases some of the elements will be superfluous, e.g., 
if one of the V; is equal to 1, the element of M connected to terminal ¢ can 


¥ ewes 
age 
ae 
el 
M —, 
¥ 4 
Fig. 7—Disjunctive tree with two bays. 


be eliminated. At worst M contains six elements and N contains two. 
The variable X appears twice, Y four times and Z twice. Of course, it is 
completely arbitrary which variables we call X, Y, and Z. We have thus 
proved somewhat more than we stated above, namely, 

Theorem 2: Any function of three variables may be realized using not more 
ihan 2, 2, and 4 elements from the three variables in any desired order. Thus 
(3) < 8, w(3) < 4. Further, since make and break elements appear in 
adjacent pairs we can obtain the distribution 1, 1, 2, in terms of transfer ele- 
ments. 

The theorem gives only upper limits for \(3) and yu(3). The question 
immediately arises as to whether by some other design method these limits 
could be lowered, i.e., can the < signs be replaced by < signs. It can be 
shown by a study of special cases that \(3) = 8, the function 


X@®Y@Z=X(¥YZ+ Y'2Z'’) + X' (¥Z' + Y'Z) 


requiring eight elements in its most economical realization. (3), however, 
is actually 3. 
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It seems probable that, in general, the function 
Xi @ Xo ® in ® a 


requires 4(m — 1) elements, but no proof has been found. Proving that a 
certain function cannot be realized with a small number of elements is 
somewhat like proving a number transcendental; we will show later that 
almost all* functions require a large number of elements, but it is difficult 
to show that a particular one does. 
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Fig. 8—Network giving all functions of one variable. 
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Fig. 9—Disjunctive tree with three bays. 
b. Functions of Four Variables: 


In synthesizing functions of four variables by the same method, two 
courses are open. First, we may expand the function as follows: 

i(W, X, Y,Z) = (W+X+ Y4+V,(Z))|-(W+ X + Y' + V2(Z)). 
[(W+ X'+ ¥ + V,(Z)]-(W + X’ + Y’ + V,(Z)I. 
[(W'+X+Y¥+V,(Z)]-(W’ +X + YY’ + V3(Z)). 
[W’ + X’ + Y +.V,(Z))-(W’ + X’ + Y’ + V;(Z)). 


By this expansion we would letU; = W+X+Y,U2=W+X+Y!’,::-, 
Us = W’ + X’ + Y’ and construct the M network in Fig. 9. N would 
* We use the expression “almost all” in the arithmetic sense: e.g., a property is true 


of almost all functions of ” variables if the fraction of all functions of m variables for which 
it is not true > Oasn— @. 
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again be as in Fig. 8, and by the same type of reasoning it can be seen that 
A(4) < 16. 

Using a slightly more complicated method, however, it is possible to 
reduce this limit. Let the function be expanded in the following way: 


f(W, X, Y,Z) =(W+ X + VAY, Z)|-(W + X’ + V2CY, Z)] 
(W’ + X + V3(¥, Z)]-(W’ + X’ + VAY, 2). 


We may use a network of the type of Fig. 7 for M. The V functions are 
now functions of two variables Y and Z and may be any of the 16 functions: 


(0 (Y (YZ (Y+Z AY + YZ’ 
d 4 | | | ¢ 
: oY oe oer eere | 
|Z | YZ’ V+Z 
| | | t 
|Z’ ty 4 ig + i A ; 


We have divided the functions into five groups, A, B, C, D and E£ for later ' 
reference. We are going to show that any function of four variables can 


a | 
al 





_ 


¥ z 


Fig. 10—Simplifying network. 


be realized with not more than 14 elements. This means that we must 
construct a network .V using not more than eight elements (since there are 
six in the M network) for any selection of four functions from those listed 
above. To prove this, a number of special cases must be considered and 
dealt with separately: 

(1) If all four functions are from the groups, A, B, C, and D, N will 
certainly not contain more than eight elements, since eight letters at most 
can appear in the four functions. 

(2) We assume now that just one of the functions is from group £; 
without loss of generality we may take it to be YZ’ + Y’Z, for it is the other, 
replacing Y by Y’ transforms it into this. If one or more of the remaining 
functions are from groups A or B the situation is satisfactory, for this func- 
tion need require no elements. Obviously 0 and 1 require no elements and 
Y, Y’,ZorZ’ may be “tapped off” from the circuit for YZ’ + Y’Z by writing 
itas(V + Z)(Y’-+ 2’). For example, Y’ may be obtained with the circuit 
of Fig. 10. This leaves four elements, certainly a sufficient number for 
any two functions from A, B, C, or D. 
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(3) Now, still assuming we have one function, YZ’ + Y’Z, from &£, 
suppose at least two of the remaining are from D. Using a similar “tapping 
off”? process we can save an element on each of these. For instance, if the 
functions are Y + Zand Y’ + Z’ the circuit would be as shown in Fig. 11. 

(4) Under the same assumption, then, our worst case is when two of the 
functions are from C and one from D, or all three from C. This latter case 
is satisfactory since, then, at least one of the three must be a term of 
YZ' + Y’Z and can be “tapped off.”” The former case is bad only when 
the two functions from C are YZ and Y’Z’. It may be seen that the only 


z 
—_——_—_—_— 
Zz Y 
2’ a 
oo 
2" 


Fig. 11—Simplifying network. 


Fig. 12—Simplifying network. 


essentially different choices for the function from D are Y + Zand Y’ + Z. 
That the four types of functions f resulting may be realized with 14 elements 
can be shown by writing out typical functions and reducing by Boolean 
Algebra. 

(5) We now consider the cases where two of the functions are from E. 
Using the circuit of Fig. 12, we can tap off functions or parts of functions 
from A, B or D, and it will be seen that the only difficult cases are the fol- 
lowing: (a) Two functions from C. In this case either the function f is 
symmetric in Y and Z or else both of the two functions may be obtained 
fromm the circuits for the & functions of Fig. 12. The symmetric case is 
handled in a later section. (b) One is from C, the other from D. There 
is only one unsymmetric.case. We assume the four functions are Y © Z, 
Y © 2’, YZ and Y + Z’. This gives rise to four types of functions f, 
which can all be reduced by algebraic methods. This completes the proof. 
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Theorem 3: Any function of four variables can be realized with not more 
than 14 elements. 


c. Functions of More Than Four Variables: 
Any function of five variables may be written 
(X1,---, Xs) = (Xs + filXi,-->, Xol(Xs + (Xi, +--+, Kol 


and since, as we have just shown, the two functions of four variables can 
be realized with 14 elements each, f(X;, --- Xs) can be realized with 30 
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Fig. 13—Network giving all functions of two variables. 








Now consider a function f(X;, X2,-°-:, Xn) of m variables. For 
5 <n < 13 we get the best limit by expanding about all but two variables. 


f(X1 ’ X2 pr Xn) = [X, + X» a ees + X n—-2 + Vi(X 1 »< fall 
ry! yl yt . r r 

Cit pees [Xi + Xe+---+ Xn-2 + Vs(Xn-1, Vad) (4) 
The V’s are all functions of the variables Y,_:, X, and may be obtained 
from the general NV network of Fig. 13, in which every function of two 
variables appears. This network contains 20 elements which are grouped 
into five transfer elements for one variable and five for the other.* The 
M network for (4), shown in Fig. 14, requires in general 2" ' — 2 elements. 
Thus we have: 


* Several other networks with the same property as Fig. 13 have been found, but they 
all require 20 elements. 
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Theorem 4. (n) < 2” + 18 


d. Upper Limits for \(m) with Large n. 


Of course, it is not often necessary to synthesize a function of more than 
say 10 variables, but it is of considerable theoretical interest to determine 
as closely as possible the behavior of \() for large n. 








Fig. 15—Network giving all functions of (m + 1) variables constructed from one giving 
all functions of m variables. 


We will first prove a theorem placing limits on the number of elements 
required in a network analogous to Fig. 13 but generalized for m variables. 

Theorem 5. An N network realizing all 2°” functions of m variables can 
be constructed using not more than 2-2°” elements, i.e., not more than two ele- 
ments per function. Any network with this property uses at least (3 — e€) 
elements per function for any € > 0 with n sufficiently large. 

The first part will be proved by induction. We have seen it to be true 
form = 1,2. Suppose it is true for some m with the network NV of Fig. 15. 
Any function of m + 1 variables can be written 


| i [X m1 + Fal LX m4 + fr 
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where f, and fy involve only m variables. By connecting from g to the cor- 
responding f, and f, terminals of the smaller network, as shown typically 
for gs, we see from Theorem 1 that all the g functions can be obtained. 
Among these will be the 2°” f functions and these can be obtained simply 
by connecting across to the f functions in question without any additional 
elements. Thus the entire network uses less than 
(er 
elements, since the V network by assumption uses less than 2-2°" and the 
first term in this expression is the number of added elements. 

The second statement of Theorem 7 can be proved as follows. Suppose 
we have a network, Fig. 16, with the required property. The terminals 
can be divided into three classes, those that have one or less elements di- 
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Fig. 16—Network giving all functions of m variables. 


rectly connected, those with two, and those with three or more. ‘The first 
set consists of the functions 0 and 1 and functions of the type 


(X + f) = X + fxao 


where X is some variable or primed variable. The number of such functions 
is not greater than 2m-2°” ' for there are 2m ways of selecting an ‘X” 
and then 2°” ‘ different functions fx—» of the remaining m — 1 variables. 
Hence the terminals in this class as a fraction of the total > 0 as m > 


Functions of the second class have the form 


g= (X + fi)(Y + fr) 


In case X ¥ Y’ this may be written 


x. 


XY + XY'gy-1,y~0 + X'Y gxeo,y=1 + X'Y' gy_o,y=0 


and there are not more than (2m)(2m — 2)(2°” *}’ such functions, again a 
vanishingly small fraction. In case Y = }¥’ we have the situation shown 
in Fig. 17 and the XX’ connection can never carry ground to another 
terminal since it is always open as a series combination. The inner ends 
of these elements can therefore be removed and connected to terminals 
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corresponding to functions of less than m variables according to the equation 


g = (X + fi)(X’ + fo) = (X + frxeo)(X’ + fox-1) 


if they are not already so connected. This means that all terminals of the 
second class are then connected to a vanishingly small fraction of the total 
terminals. We can then attribute twe elements each to these terminals 
and at least one and one-half each to the terminals of the third group. As 
these two groups exhaust the terminals except for a fraction which — 0 
asn— 2, the theorem follows. 

If, in synthesizing a function of ” variables, we break off the tree at the 
(n — m)th bay, the tree will contain 2" "*’ — 2 elements, and we can find 
an N network with not more than 2°” -2 elements exhibiting every function 
of the remaining m variables. Hence 
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Fig. 17—Possible situation in Fig. 16. 


for every integer m. We wish to find the integer M = M(n) minimizing 
this upper bound. 
Considering m as a continuous variable and n fixed, the function 


f(m) = 2°74 2.2 


clearly has just one minimum. This minimum must therefore lie between 
m, and m,; + 1, where 


f(m,) = f(m, + 1) 


om,+1 


i.e., ag m,+l + , ie _ ia +. o% l a, 


+1 9m 
— 2 *) 


or -_ iss Y aie | ili 


Now m, cannot be an integer since the right-hand side is a power of two 
and the second term is less than half the first. It follows that to find the 
integer M making f(M) a minimum we must take for M the least integer 
satisfying 


9M +1 


fave % 
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Thus M satisfies: 
M+1+2"¥"* >n>M+2”™ 

This gives: 
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Fig. 18—Behaviour of g(n). 


n+1 


Our upper bound for \(m) behaves something like a with a superimposed 


saw-tooth oscillation as ” varies between powers of two, due to the fact that 
m must be an integer. If we define g(n) by 


n+l 


gett "2 = g(n) = ' 


M being determined to minimize the function (i.e., M satisfying (5)), then 
g(n) varies somewhat as shown in Fig. 18 when plotted against log, n. The 
maxima occur just beyond powers of two, and closer and closer to them 
asn-—»«. Also, the saw-tooth shape becomes more and more exact. The 
sudden drops occur just after we change from one value of M to the next. 
These facts lead to the following: 

Theorem 6. (a) For all n 


)? 
Mn) < = 


(6) For almost all n 
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(c) There is an infinite sequence of n; for which 


y thle 


(ns) < a (1 + «) 
These results can be proved rigorously without much difficulty. 
e. A Lower Limit for A(m) with Large n. 


Up to now most of our work has been toward the determination of upper 
limits for \(m). We have seen that for all n 


Mn) < B. 
n 


. ; ; ’ oP 
We now ask whether this function B — is anywhere near the true value 
n 


of A(m), or may A(m) be perhaps dominated by a smaller order of infinity, 
e.g., n?. It was thought for a time, in fact, that \(m) might be limited by 
n’ for all n, arguing from the first few values: 1, 4, 8, 14. We will show that 


n 


ih a ; 
this is far from the truth, for actually — is the correct order of magni- 
n 


tude of A(m): 


ri 2” 
n 


A— <Xn) <B 
n 


for all . A closely associated question to which a partial answer will be 
given is the following: Suppose we define the “complexity” of a given func- 
tion f of m variables as the ratio of the number of elements in the most 
economical realization of f to A(m). Then any function has a complexity 
lying between 0 and 1. Are most functions simple or complex? 

Theorem 7: For all sufficiently large n, all functions of n variables excepting 


a fraction 6 require at least (1 — €) — elements, where € and 6 are arbitrarily 
n 


small positive numbers. Hence for large n 


Mn) > (1 -o= 


and almost all functions have a complexity > }(1 — €). For acertain sequence 
n; almost all functions have a complexity > }(1 — €). 

The proof of this theorem is rather interesting, for it is a pure existence 
proof. We do not show that any particular function or set of functions 


n 


; 2 ed ‘ . 
requires (1 — e) — elements, but rather that it is impossible for all functions 
n 
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to require less. This will be done by showing that there are not enough 


n 


; 2 , 
networks with less than (1 — e) _ branches to go around, i.e., to represent 


all the 2°” functions of variables, taking account, of course, of the different 
assignments of the variables to the branches of each network. This is only 
possible due to the extremely rapid increase of the function 2°". We require 
the following: 

Lemma: The number of two-terminal networks with K or less branches is 
less than (6K)*. 

Any two-terminal network with A or less branches can be constructed 
as follows: First line up the A branches as below with the two terminals 
a and 6. 


a. 1—1’ 
2-—2’ 
3—3’ 
4—4’ 
b. K—K’ 
We first connect the terminals a, 6, 1, 2, --- , A together in the desired way. 


The number of different ways we can do this is certainly limited by the num- 
ber of partitions of K + 2 which, in turn, is less than 

yk +1 
for this is the number of ways we can put one or more division marks between 
the symbols a, 1, --- , K, 6. Now, assuming a, 1, 2, --- , K, 6, intercon- 
nected in the desired manner, we can connect 1’ either to one of these ter- 
minals or to an additional junction point, i.e., 1’ has a choice of at most 


K + 3 


terminals, 2’ has a choice of at most K + 4, etc. Hence the number of 
networks is certainly less than 


2°*'(K + 3) (K + 4) (K + 5) --+ (2K + 3) 
<(6K)* K>3 


and the theorem is readily verified for A = 1, 2. 
We now return to the proof of Theorem 7. The number of functions of 


(1 — e)2” 


variables that can be realized with ~ elements is certainly less than 
1 


the number of networks we can construct with this many branches multi- 
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plied by the number of assignments of the variables to the branches, i.e., 
it is less than 


)(on 7" (1—e) (2"/n) 
H = (25) O™ | ct — ¢) =| 
n 
Hence 
2" 2” 2 
RE ee ae ee ea 
n 


= (1 — e) 2" + terms dominated by this term for large n. 


By choosing n so large that ~ 2" dominates the other terms of log H we 


dol a 


arrive at the inequality 
logs H < (1 — ©) 2” 
ee deal 
But there are S = 2°” functions of n variables and 


H 9a €,)2" 


one Se —() as no ow, 
S 22" : 
Hence almost all functions require more than (1 — €,)2" elements. 
Now, since for all x > N there is at least one function requiring more than 


n 


2" ; 
(say) ; ~ elements and since A(z) > 0 for n > 0, we can say that for all n, 
) 


M(n) > A fl 
il 


for some constant A > 0, for we need only choose A to be the minimum 
number in the finite set: 


1 \(1) (2) M3) ACN) 
>? ie , i ies a 
1 2 3 N 


Thus A(m) is of the order of magnitude of = The other parts of Theorem 


8 follow easily from what we have already shown. 

The writer is of the opinion that almost all functions have a complexity 
nearly 1,i.e., > 1— ¢€. This could be shown at least for an infinite sequence 
n; if the Lemma could be improved to show that the number of networks is 
less than (6K)"” for large K. Although several methods have been used 
in counting the networks with K branches they all give the result (6A)¥*. 
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It may be of interest to show that for large K the number of networks is 
greater than 

(6K)*"* 
This may be done by an inversion of the above argument. Let f(A) be the 


° > - ° gn . 
number of networks with K branches. Now, since there are 2” functions of 
Jn+2 
n variables and each can be realized with (1 + e) elements ( suthciently 
n 4 


large), 


i(a + «€) = alias cs" 


for n large. But assuming f(K) < (6K)*" reverses the inequality, as 
is readily verified. Also, for an infinite sequence of K, 


{(K) > (6K)*” 


Since there is no obvious reason why /(A) should be connected with powers 
of 2 it seems likely that this is true for all large K. 
We may summarize what we have proved concerning the behavior of 


; : om 
\(n) for large nas follows. (n) varies somewhat as ; if we let 
n 


n+1 
(mn) = Az 
n 
then, for large n, A, lies between } — e¢ and (2 + e), while. for an infinite 
sequence of m,} — € << An< 1+. 

We have proved, incidentally, that the new design method cannot, in a 
sense, be improved very much. With series-parallel circuits the best known 
limit* for A(n) is 

A(n) < 3.27! + 2 


n 


and almost all functions require (1 — e) : elements. We have lowered 
1 


logs 


the order of infinity, dividing by at least and possibly by n. The 


1 
logs m 


best that can be done now is to divide by a constant factor < 4, and for 


some n, < 2. The possibility of a design method which does this seems, 
however, quite unlikely. Of course, these remarks apply only to a perfectly 
general design method, i.e., one applicable to any function. Many special 
classes of functions can be realized by special methods with a great saving. 

* Mr. J. Riordan has pointed out an error in my reasoning in (6) leading to the statement 
that this limit isactually reached by the function X; ®@ X2@... @® Xn, and has shown that 
this function and its negative can be realized with about * elements. The error occurs 


in Part IV after equation 19 and lies in the assumption that the factorization given is 
the best. 
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PART II: CONTACT LOAD DISTRIBUTION 


4. FUNDAMENTAL PRINCIPLES 


We now consider the question of distributing the spring load on the relays 
as evenly as possible or, more generally, according to some preassigned 
scheme. It might be thought that an attempt to do this would usually 
result in an increase in the total number of elements over the most economi- 
cal circuit. This is by no means true; we will show that in many cases (in 
fact, for almost all functions) a great many load distributions may be ob- 
tained (including a nearly uniform distribution) while keeping the total 
number of elements at the same minimum value. Incidentally this result 
has a bearing on the behavior of y(n), for we may combine this result with 


¥ -— 





Fig. 19—Disjunctive tree with the contact distribution 1, 3, 3. 
, , , a 
preceding theorems to show that u(m) is of the order of magnitude of 7 as 
n 
n-—» © and also to get a good evaluation of u(m) for small n. 

The problem is rather interesting mathematically, for it involves additive 
number theory, a subject with few if any previous applications. Let us 
first consider a few simple cases. Suppose we are realizing a function with 
the tree of Fig. 9. The three variables appear as follows: 

ae ae appear 

2, 4, 8 times, respectively: 
or, in terms of transfer elements* 

1, 2, 4. 
Now, W, X, and Y may be interchanged in any way without altering the 
operation of the tree. Also we can interchange X and Y in the lower branch 
of the tree only without altering its operation. This would give the dis- 
tribution (Fig. 19) 

1, 3,3 


* In this section we shall always speak in terms of transfer elements. 
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A tree with four bays can be constructed with any of the following dis- 
tributions 
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= 
w% 
< 


“ 
~~ 
~ 


Ct 


. 

- 
oe 
es 
—= — 
-. ~} 
es eS 

w > 


) 


J 
- 


+ 
. 
~~ 
~ 

- 
- 


- 
~ 
~ 
~~ 


wWwnnd bd 
WN 
Daan TteceOoOa 7100 N 


Ww 
ue 
wrn bd 


ee et 
— hh 

Ww B= W DY ND bO HO 
mw wd wo - — —- 
— eG OND 
Ww We DO 


a 
w 


~ 


++4+4+4++44 


on 
wn 
oo 


and the variables may be interchanged in any manner. The ‘‘sums”’ on the 
right show how these distributions are obtained. The first set of numbers 
represents the upper half of the tree and the second set the lower half. They 
are all reduced to the sum of sets 1, 2, 4 or 1, 3, 3 in some order, and these 
sets are obtainable for trees with 3 bays as we already noted. In general it is 
clear that if we can obtain the distributions 


@1 , Go, G3, °°* , an 
bi, b2, bs, +--+, dn 
for a tree with m bays then we can obtain the distribution 
1, a1 + 61, a2 + bo,°-- ,an+ by 


for a tree with x + 1 bays. 

Now note that all the distributions shown have the following property: 
any one may be obtained from the first, 1, 2, 4, 8, by moving one or more 
units from a larger number to a smaller number, or by a succession of such 
operations, without moving any units to the number 1. Thus 1, 3, 3, 8 is 
obtained by moving a unit from 4 to 2._ The set 1, 4, 5, 5 is obtained by 
first moving two units from the 8 to the 2, then one unit to the 4. Further- 
more, every set that may be obtained from the set 1, 2, 4, 8 by this process 
appears as a possible distribution. This operation is somewhat analogous 
to heat flow—heat can only flow from a hotter body to a cooler one just as 
units can only be transferred from higher numbers to lower ones in the above. 

These considerations suggest that a disjunctive tree with » bays can be 


constructed with any load distribution obtained by such a flow from the 
initial distribution 


a ee i 


We will now show that this is actually the case. 
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First let us make the following definition: The symbol (a; , a2, --- , dn) 
represents any set of numbers 0,, d2, --- , 6, that may be obtained from 
the set a; , @2, --- , a, by the following operations: 

1. Interchange of letters. 

2. A flow from a larger number to a smaller one, no flow, however, being 
allowed to the number 1. Thus we would write 


1, 2, 4, 8 = (1, 2, 4, 8) 
4, 4,1, 6 = (1, 2, 4, 8) 
1, 3, 10, 3, 10 = (1, 2, 4, 8, 12) 
but 2, 2 ¥ (1,3). It is possible to put the conditions that 
bi, be, +++ 5 bn = (GQ), G2, -** , an) (6) 


into a more mathematical form. Let the a; and the 6; be arranged as non- 
decreasing sequences. Then a necessary and sufficient condition for the 
relation (6) is that 


(1) Yb > Da 


i=l 


(2) >> 6; = Do a;, and 
1 1 


(3) There are the same number of 1’s among the a; as among the 0;. 
The necessity of (2) and (3) is obvious. (1) follows from the fact that if 
a; is non-decreasing, flow can only occur toward the left in the sequence 


Q,, 42,43,°** , an 


8 


and the sum >, a; can only increase. Also it is easy to see the sufficiency of 
I 


the condition, for if b; , b2 , --- , b, satisfies (1), (2), and (3) we can get the 
b; by first bringing a; up to 4; by a flow from the a; as close as possible to 
a, (keeping the “entropy” low by a flow between elements of nearly the 
same value), then bringing a2 up to d2 (if necessary) etc. The details 
are fairly obvious. 

Additive number theory, or the problem of decomposing a number into the 
sum of numbers satisfying certain conditions, (in our case this definition is 
generalized to “‘sets of numbers’’) enters through the following Lemma: 

Lemma: If a1, @2, -+- ,@n = (2,4, 8, +--+ , 2”) 
then we can decompose the a; into the sum of two sets 


@a,=)b;+¢; 
such that 
b,, bo, nee » = (1, 2, 4, eee i 
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and 
Cop Mne**oGa@ CG, 2.4 °° 
We may assume the a; arranged in a non-decreasing sequence, a; < a3 < 
ds S +++ < ay. In case a; = 2 the proof is easy. We have 
L243, 9" 
ES ae 





2,4,8,-+»,2* 
and a flow has occurred in the set 
436 ---.3 


to give d2, ds, °-:,@,. Now any permissible flow in C corresponds to a 
permissible flow in either A or B since if 


¢ = a,+6;> ¢ = a, t+ bd; 
then either a,>a; or b>; 


Thus at each flow in the sum we can make a corresponding flow in one or the 
other of the summands to keep the addition true. 
Now suppose a; > 2. Since the a; are non-decreasing 


in tha <2 —- 2-4 Ss ~~ 2-3 


Hence 
a+l = 
2 ~- po, 1 < et 


ls ON aa 
” n—1 


the last inequality being obvious for m > 5 and readily verified for n < 5. 
This shows that (a; — 1) and (az — 1) lie between some powers of two in the 
set 

1, 2,4, jnm ,2~ 
Suppose 

2e-1 < (a, — 1) < 2° 


2-1 < (a, — 1) < a gqiips(n— 1). 


Allow a flow between 2% and 2%" until one of them reaches (a; — 1), the 
other (say) R; similarly for (ag — 1) the other reaching S. As the start 
toward our decomposition, then, we have the sets (after interchanges) 


ee. ote 
Finer FOS FP ne 
’ eer 





2 
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We must now adjust the values to the right of L — L to the values 
dg,@4,-*-,@,. Let us denote the sequence 


4, 8, +++, 204, (20-1 + R), 3-29, 3-2e+, ... (27 + S$), QPt?,--- 20 


by “1, #2, °** , Mn-2. Now since all the rows in the above addition are 
non-decreasing to the right of L — L, and no 1’s appear, we will have proved 
the lemma if we can show that 

4 t+3 

2. te S 2h i= 1,2,---,(m — 2) 

i=l t=3 
since we have shown this to be a sufficient condition that 

a, 44,°°* » an = (M1, Hn, es » Un—2) 


and the decomposition proof we used for the first part will work. 
i<q-— 2,i.e., before the term (2¢' + R) 


D ws = 4(2" — 1) 


i=] 
$+3 
> a; > tag > 127" > i" 
3 
qsPp 


i+3 


Lu S Lia 


Next, for (¢g — 1) <i < (p — 3), i.e., before the term (2? + S) 


D ws = 427" — 1) + RF 3-22" — 1) 
1 


< 3-7" —~4<¢ 3-7" ~—§ 
since 
R< 29 


also again 
i+3 


> a > iv" 
3 


so that in this interval we also have the desired inequality. Finally for the 
last interval, 
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Dw = 2 - a — a <2 -—a—ay—-2 


and 
i+3 i+3 a 
> a; = oa; — a, — a > 2°" —a, — a — 2 
3 1 

since 


@,,42,°°', an = (2, 4, 8, +--+ , 2”) 


This proves the lemma. 


5. Tue DisyuncTIvE TREE 
It is now easy to prove the following: 
Theorem 8: A disjunctive tree of n bays can be constructed with any dis- 
tribution 
Gi Bay? 5 Oa (1, 2, 4, ara ="). 


We may prove this !by induction. We have seen it to be true for n = 


2,3,4. Assuming it for 7, it must be true form + 1 since the Lemma shows 
that any 
a Mia Tia (2, 4, 8, mies ty! , 3°) 
can be decomposed into a sum which, by assumption, can be realized for the 
two branches of the tree. 
It is clear that among the possible distributions 


—1 
(1, 2, 4, paige ) 
for the tree, an ‘“‘almost uniform” one can be found for all the variables but 


one. That is, we can distribute the load on (nm — 1) of them uniformly 
except at worst for one element. We get, in fact, for 


n= 1 1 

n=2 1,2 

n= 3 was 

n= 4 4; 3,3, 

n= 5§ rer S 

n=6 L, 33, 33, 12, 88:23 
n=7 3, Zh, 28, al, 20, Oey 28 
etc. 


as nearly uniform distributions. 


6. OTHER DISTRIBUTION PROBLEMS 


Now let us consider the problem of load distribution in series-paralle] 
circuits. We shall prove the following: 



































Sa S 
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Theorem 9: Any function f(X,, X2, +--+ , Xn) may be realized with a 


-) 


series-parallel circuit with the following distribution: 
(1, 2 4, ee ar*), Qn-2 


in terms of transfer elements. 
This we prove by induction. It is true for m = 3, since any function of 
three variables can be realized as follows: 


AX, Y,Z) =(X + fa (Y, DIX’ + fe (Y, Z)I 


and fi(Y, Z) and fo(Y, Z) can each be realized with one transfer on Y and 
one on Z. Thus f(X, Y, Z) can be realized with the distribution 1, 2, 2. 
Now assuming the theorem true for (x — 1) we have 


flv, Xe, <°*» Ee) = ve ¥ fe he: , Mead 
[Xe + f£AX1, Xe, --> , Xa-ad] 
and ) 


A 4, 8, tee Qn-3 
4 8, ara ——_ 


dN bo 


’ bd 


4, 8, 16,--- , 2"? 
A simple application of the Lemma thus gives the desired result. Many 
distributions beside those given by Theorem 9 are possible but no simple 
criterion has yet been found for describing them. We cannot say any 
distribution 


(1, 2,4 &*- 2n-2) 


In—2 
9 = , 


(at least from our analysis) since for example 
3, 6, 6, 7 = (2, 4, 8, 8) 
cannot be decomposed into two sets 
@, , G2, 3, a4 = (1, 2, 4, 4) 
and 
bi, bo, bs, bs = (A, 2, 4, 4) 


It appears, however, that the almost uniform case is admissible. 

As a final example in load distribution we will consider the case of a net- 
work in which a number of trees in the same variables are to be realized. 
A large number of such cases will be found later. The following is fairly 
obvious from what we have already proved. 
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Theorem 10: It is possible to construct m different trees in the same n variables 
with the following distribution: 


1, 0g,°°* » an = (m, 2m, 4m, eee , 2"-!m) 


It is interesting to note that under these conditions the bothersome 1 disap- 
pears form > 1. Wecan equalize the load on all » of the variables, not just 
n — 1 of them, to within, at worst, one transfer element. 

7, THe FuNcTION y(n) 


We are now in a position to study the behavior of the function u(m). 
This will be done in conjunction with a treatment of the load distributions 
possible for the general function of m variables. We have already shown 
that any function of three variables can be realized with the distribution 


14,2 
in terms of transfer elements, and, consequently yu(3) < 4. 
Any function of four variables can be realized with the distribution 
1, 1, (2, 4) 
Hence u(4) < 6. For five variables we can get the distribution 
1, 1, (2, 4, 8) 
or alternatively 
1, 5, 5, (2, 4) 
so that u(5) < 10. With six variables we can get 
1, 5, 5, (2, 4, 8) and (6) < 10 
for seven, 
1, 5, 5, (2, 4, 8, 16) and u(7) < 16 
etc. Also, since we can distribute uniformly on all the variables in a tree 
except one, it is possible to give a theorem analogous to Theorem 7 for the 


function u(n): 
Theorem 11: For all n 


lid 


n> 


u(n) 


lA 


For almost all n 
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AP 
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For an infinite number of n; , 


n+l 


sens S* 


n? 


The proof is direct and will be omitted. 


PART III: SPECIAL FUNCTIONS 


8. FUNCTIONAL RELATIGNS 


We have seen that almost all functions require the order of 


get 
n2 


elements per relay for their realization. Yet a little experience with the 
circuits encountered in practice shows that this figure is much too large. 
In a sender, for example, where many functions are realized, some of them 
involving a large number of variables, the relays carry an average of perhaps 
7 or 8 contacts. In fact, almost all relays encountered in practice have less 
than 20 elements. What is the reason for this paradox? The answer, of 
course, is that the functions encountered in practice are far from being a 
random selection. Again we have an analogue with transcendental numbers 
—although almost all numbers are transcendental, the chance of first en- 
countering a transcendental number on opening a mathematics book at 
random is certainly much less than 1. The functions actually encountered 
are simpler than the general run of Boolean functions for at least two major 
reasons: 

(1) A circuit designer has considerable freedom in the choice of functions 
to be realized in a given design problem, and can often choose fairly simple 
ones. For example, in designing translation circuits for telephone work it is 
common to use additive codes and also codes in which the same number of 
relays are operated for each possible digit. The fundamental logical simplic- 
ity of these codes reflects in a simplicity of the circuits necessary to handle 
them. 

(2) Most of the things required of relay circuits are of a logically simple 
nature. The most important aspect of this simplicity is that most circuits 
can be broken down into a large number of small circuits. In place of 
realizing a function of a large number of variables, we realize many functions, 
each of a small number of variables, and then perhaps some function of these 
functions. ‘To get an idea of the effectiveness of this consider the following 
example: Suppose we are to realize a function 


f(X1, X2, Br a Xen) 
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of 2n variables. The best limit we can put on the total number of elements 
q2n +1 


necessary is about 7 . However, if we know that f is a function of two 


~~ 


functions f; and fz, each involving only » of the variables, i.e. if 
f g(hi , fa) 
frm fii, B+, Xe) 
h= fl Xn41 op, Ore Pee Xon) 


then we can realize f with about 


an 1 
‘.* 

n 
a) 
b hag pent : : y 
elements, a much lower order of infinity than in If gis one of the simpler | 
n ‘ 
} 
functions of two variables; for example if g(f1 , f2) = fi + fe, or in any case j 
at the cost of two additional relays, we can do still better and realize f with { 


Jnt+l j 
about 2 elements. In general, the more we can decompose a synthesis 
n 


problem into a combination of simple problems, the simpler the final circuits. 
The significant point here is that, due to the fact that f satisfies a certain 
functional relation 


f - g(hi » te), 


we can find a simple circuit for it compared to the average function of the 
same number of variables. 

This type of functional relation may be called functional separability. It 
is often easily detected in the circuit requirements and can always be used 
to reduce the limits on the number of elements required. We will now show t 
that most functions are not functionally separable. 

Theorem 12: The fraciion of all functions of n variables that can be written 
in the form 4 





f = g(h(Xi +++ X.), Xen, +++ » Xn) 
where 1 << s <n -- I approaches zero as n approaches ~. 
. me ; 
We can select the s variables to appear in / in( ways; the function h 
P ’ 


Q8 ee ae gn-satl ey spe, 0 . ° 

then has 2° possibilities and g has 2 possibilities, since it has m — s+ 1 
arguments. The total number of functionally separable functions is there- 
fore dominated by 
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n—2 n tale 
> ( ) 9" 9" * 
§ 


s=2 


and the ratio of this to 2” > 0asn— ~. 


Xs+1 





en 
SS Lhe } oo 
QR _ 


Fig. 21—Use of separability of two sets of variables. 


In case such a functional separability occurs, the general design method 
described above can be used to advantage in many cases. This is typified 
by the circuit of Fig. 20. If the separability is more extensive, e.g. 


f = gay X, ase Sieh Nhe X 41 rs Jesh. X 141 Fide Ai) 


the circuit of Fig. 21 can be used, using for ‘‘/2” either 4; or hy , whichever 
requires the least number of elements for realization together with its 
negative. 

We will now consider a second type of functional relation which often 
occurs in practice and aids in economical realization. This type of relation 
may be called group invariance and a special case of it, functions symmetric 
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in all variables, has been considered in (6). A function f(X1, +--+, Xa) 
will be said to be symmetric in XY, , X¢ if it satisfies the relation 


AX, Me, «+, Re ls Be +, BD. 


It is symmetric in XY, and X’, if it satisfies the equation 


f(Xi, Xe,-°+, Xe) = f(Xa, Xi, Xs, °° Xs) 


These also are special cases of the type of functional relationships we will 

consider. Let us denote by 

Noo +++ o = I the operation of leaving the variables in a function as they 
are, 

Nico *** the operation of negating the first variable (i.e. the one occupy- 
ing the first position), 

Now +++ that of negating the second variable, 

Niio+++ that of negating the first two, etc. 

So that Nioif(X, Y, Z) = f(X’YZ’) etc. 

The symbols Ni form an abelian group, with the important property that 
each element is its own inverse; NiNi = J The product of two elements 
may be easily found — if V; VN; = N;, , & is the number found by adding 7 
and 7 as though they were numbers in the base two but without carrying. 

Note that there are 2" elements to this “negating” group. Now let 

= ] = the operation of leaving the variables of a function in the 
same order 
be that of interchanging the first two variables 
that of inverting the order of the first three, etc. 


S32 f(X, Y; Z) i K(Z, > Y) 
Ss2f(Z, X,Y) = Stu f(X, Y, Z) = f(Y, Z, X) 


etc. The S; also form a group, the famous “substitution” or “symmetric” 
group. It isofordern!. Itdoes not, however, have the simple properties 
of the negating group—it is not abelean (m > 2) nor does it have the self 
inverse property.* The negating group is not cyclic if m > 2, the symmetric 
group is not if n > 3. 

The outer product of these two groups forms a group G whose general 
element is of the form VS; and since i may assume 2" values and j, 2! values, 
the order of G is 2"n! 

It is easily seen that S; NV; = N,5S;, where k may be obtained by per- 


* This is redundant; the self inverse property implies commutativity for if XX = / 
then X¥ = (XY) = YX" = YX. 
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forming on 7, considered as an ordered sequence of zero’s and one’s, the 
bd , 
permutation S;. Thus 


Sasa N i100 be N 1010 Soai4 . 


By this rule any product such as N;S; Vi Ni Sm N.S, can be reduced to the 
form 


NN; +++ NaSpSq°*+ St 


and this can then be reduced to the standard form NS; . 
A function f will be said to have a non-trivial group invariance if there are 
elements V,S; of G other than J such that identically in all variables 


N,S; f = f. 


It is evident that the set of all such elements, V;S; , for a given function, 
forms a subgroup G, of G, since the product of two such elements is an ele- 
ment, the inverse of such an element is an element, and all functions are 
invariant under J. 

A group operator leaving a function f invariant implies certain equalities 
among the terms appearing in the expanded form of f. To show this, 
consider a fixed N,S;, which changes in some way the variables (say) 
Xi, X2,-°-:, X,. Let the function /(Xi,---, X,) be expanded 
about X,,--:, X,: 


f=(Mit Xat e+ Xt filXar, ++, Xe) 
[Xi + Xe + ++) +X, + fol(Xri, +++, Xa) 


2 ee ce eee oe oe) 


If f satishes V;S,f = f we will show that there are at least 42" equalities 
between the functions fi, f2,---, for. Thus the number of functions 
satisfying this relation is 

< ili oat i 9 


since each independent f; can be any of just 2” ‘ functions, and there are 
at most } 2’ independent ones. Suppose V;5; changes 


} ae mere 2 A 


* a 


,* 
Rus Maa” ** Me B 


r 


P e ° ,* , + 
where the *’s may be either primes or non primes, but no X,; = X;. Give 
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X, the value 0. This fixes some element in B namely, X,; where a; = 1. 
There are two cases: 
(1) If this element is the first term, a; = 1, then we have 


| ne 3 
( ee 7 


Letting X2,--- , X, range through their 2” * possible sets of values gives 
2""* equalities between different functions of the set f; since these are 
really 


S(X1, Xa, 00+ Xe, Xegry t+» Xn) 


with X¥,, X2,--- , X, fixed at a definite set of values. 

(2) If the element in question is another term, say X,, , we then give X»2 
in line A the opposite value, X2 = zy = (Xs )’. Now proceeding as 
before with the remaining r — 2 variables we establish 2” * equalities between 
the f;. 


Now there are not mo an 2"! relatio 
N tl t re* than 2"! relations 
N iS; f = f 


of the group invariant type that a function could satisfy, so that the number 
of functions satisfying any non-trivial relation 
< 2"! 2” 

Since 
— 0 asn—> &% 
we have: 

Theorem 13: Almost all functions have no non-trivial group invariance. 

It appears from Theorems 12 and 13 and from other results that almost all 
functions are of an extremely chaotic nature, exhibiting no symmetries or 
functional relations of any kind. This result might be anticipated from the 
fact that such relations generally lead to a considerable reduction in the 
number of elements required, and we have seen that almost all functions are 
fairly high in ‘“‘complexity”’. 

If we are synthesizing a function by the disjunctive tree method and the 
function has a group invariance involving the variables 


X1,X2, >To » ae 
at least 2’ “ of the terminals in the corresponding tree can be connected to 
* Our factor is really less than this because, first, we must exclude N;.S; = J; and second, 


except for self inverse elements, one relation of this type implies others, viz. the powers 
(Ni S;)?°f = f. 
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other ones, since at least this many equalities exist between the functions to 
bé joined to these terminals. This will, in general, produce a considerable 
reduction in the contact requirements on the remaining variables. Also an 
economy can usually be achieved in the M network. In order to apply this 





~ INDEPENDENT > 
x,Y 
y’ 
(3) (4) 


Fig. 22—Networks for group invariance in two variables. 


Fig. 23—Networks for group invariance in three variables. 


method of design, however, it is essential that we have a method of deter- 
mining which, if any, of the V;S; leave a function unchanged. The 
following theorem, although not all that might be hoped for, shows that we 
don’t need to evaluate V;S;f for all V;S; but only the V;f and S;f. 

Theorem 14: A necessary and sufficient condition that N;S;f = f is that 
Ni f = S;f. 


This follows immediately from the self inverse property of the V;. Of 
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course, group invariance can often be recognized directly from circuit re- 
quirements in a design problem. 

Tables I and II have been constructed for cases where a relation exists 
involving two or three variables. To illustrate their use, suppose we have 
a function such that 


Nii Ssia f - f 


Yo i i 


=] 


Fig. 24—M network for partially symmetric functions. 


The corresponding entry Z’Y’X in the group table refers us to circuit 9 of 
Fig. 23. The asterisk shows that the circuit may be used directly; if there 
is no asterisk an interchange of variables is required. We expand f about 
X, Y, Z and only two different functions will appear in the factors. These 
two functions are realized with two trees extending from the terminals of the 
network 9. Any such function f can be realized with (using just one variable 
in the V network) 


9 + 2(2"* — 2) 42 


= get + 7 elements, 
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a much better limit than the corresponding 


2-1 4 18 
for the general function. 


TABLE I 


Group INVARIANCE INVOLVING Two VARIABLES (SUPERSCRIPTS REFER TO FiG. 22) 
Si Sor 

No (x y) (y x)* 

Noi (x y’)** (y x’)3* 

Ni (x’y)? (y’x)** 

Nu (x’y’)4 (y’x’)! 

TABLE II 

Group INVARIANCE INVOLVING THREE VARIABLES (SUPERSCRIPTS REFER TO Fic. 23) 
Nui wie XZY! | 2. ¢ ag VZX 2* aX Y * ZYX : 
Nor XYZ'* D Ph hie Faz * ¥2x" re ie LYX"* 
Nic 2 Ee" ALY” YA’2 ¥2 2 * aax* ZY'X? 
Noun XY'Z’® Aa Yaa * ee te Zan * ara * 
Niow X’YZ3 Aen ¥’ XZ * ¥.25 ° Zar * £Yx* 
Nio A’¥Z * X’ZY'* yAze = 2" £xAY* y i Ble. 
Ni X'Y'Z* be a | oe. a ya." aay? ZY? 
Nia a Y's aay As Tou dl, tl) ii ois 


9, PARTIALLY SYMMETRIC FUNCTIONS 
We will say that a function is “partially symmetric” or “symmetric in a 
certain set of variables” if these variables may be interchanged at will 
without altering the function. Thus 
NYVZW + (XVY’ + X’YW + WZ’ 
is symmetric in Y and ¥. Partial symmetry is evidently a special case of 
the general group invariance we have been considering. It is known that 
any function symmetric in all variables can be realized with not more than 
a P e ° 6 ° P P 
n? elements, where 1 is the number of variables.’ In this section we will 
improve and-generalize this result. 
Theorem 15: Any function f(X,, Xo, +++, Xm, V1, Yo,+++ Vn) sym- 
metric in Ny, X2,-+++, Xm can be written 


Fea 4 BGs *** 5 Mea 3 Sst" 


si [So(X1 ’ X2 oo oy X m) + FAV 5) Y, Ri oes )}. 
[SCX ’ X» en Au) + mee , Y» ee ae Fe )]. 


[Sm(X1, X2,°°*, Xm) + fm(¥1, Ye,---, Val (6) 


where 


Jal2 ss, Va5-°> 5 Fa) 
= f(0,0,--- , 0, See | Tas 2ee*** 3a 
k0’s (m — k) 1’s 
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and S;(X,, X2, +++, Xm) is the symmetric function of X,, X2,°+:, Xn 
with k for its only a-number. 

This theorem follows from the fact that since f is symmetric in XY, , X2, 

- , Y,, the value of f depends only on the number of X’s that are zero and 
the values of the ]’’s. If exactly A of the X’s are zero the value of / is 
therefore fx , but the right-hand side of (6) reduces to fx in this case, since 
then S;(X,, X2,--: ,Xm) = 1,7 ¥ K,and Sx = 0. 

The expansion (6) is of a form suitable for our design method. We can 
realize the disjunctive functions Sx(X,, Yo, +--+ , Y,) with the symmetric 
function lattice and continue with the general tree network as in Fig. 24, 
one tree from each level of the symmetric function network. Stopping the 
trees at J’; , it is clear that the entire network is disjunctive and a second 
application of Theorem 1 allows us to complete the function f with two ele- 


ments from }’,,. Thus we have 
Theorem 16. Any function of m + n variables symmetric in m of them can 


be realized with not more than the smaller of 
(m + 1)(AGiz) + m) or (m 4- 1)(2" + m — 2) + 2 


elements. In particular a function of n variables symmetric in n — 2 or more 
of them can be realized with not more than 


n—n+2 


elements. 

If the function is symmetric in Y,, V2, ---,X,,andalsoinY,,V2,---, 
)’,, and not in Z;, Z2,--- , Z, it may be realized by the same method, 
using symmetric function networks in place of trees for the V’ variables. 
It should be expanded first about the Y’s (assuming m < r) then about the 
Ys and finally the Z’s. The Z part will be a set of (m + 1)(r + 1) trees. 


REFERENCES 
. G. Birkhoff and S. MacLane, “A Survey of Modern Algebra,” Macmillan, 1941. 
. L. Couturat, “The Algebra of Logic,’’ Open Court, 1914. 
. J. H. Woodger, ‘‘The Axiomatic Method in Biology,’’ Cambridge, 1937. 
. W.S. McCulloch and W. Pitts, ““A Logical Calculus of the Ideas Immanent in Nervous 
Activity,” Bull. Math. Biophysics, V. 5, p. 115, 1943. 
I. C. Berkeley, “Boolean Algebra and Applications to Insurance,” Record (American 
Institute of Actuaries), V. 26, p. 373, 1947. 
. C. E. Shannon, ‘A Symbolic Analysis of Relay and Switching Circuits,’ Trans. A. 
fl. EE. V. 31,3. 113; 1988: 
. J. Riordan and C. E. Shannon, ““The Number of Two-Terminal Series Parallel Net- 
works,” Journal of Mathematics and Physics, V. 21, No. 2, p. 83, 1942. 
8. A. Nakashima, Various papers in Vippon Electrical Communication Engineering, April, 
Sept., Nov., Dec., 1938. F 
. H. Piesch, Papers in Archiv. from Electrotechnik XXXII, p. 692 and p. 733, 1939. 
. G. A. Montgomerie, ‘Sketch for an Algebra of Relay and Contactor Circuits,” Jour. 
I. of E. E., V.95, Part IT, No. 36, July 1948, p. 303. 
. G. Polya, “Sur Les Types des Propositions Composées,” Journal of Symbolic Logic, 
V.5, No. 3, p. 98, 1940. 





A Method of Measuring Phase at Microwave Frequencies 


By SLOAN D. ROBERTSON 
A method of measuring microwave phase differences is described in which it 
is unnecessary to compensate for amplitude inequalities between the signals 
whose phases are being compared. The apparatus described is alsé suited for 
the measurement of the magnitude of a transfer impedance as well as the phase. 
YITH the increasing interest in wide-band amplifiers and circuits for 
microwave communication systems the measurement of the transfer 
phases of such components has become a necessary procedure. A commonly 
used technique for measuring phase at microwave frequencies is to sample 
the signal at the input and output of the device to be measured and to obtain 
a nul! balance between the two signals by varying the phase of one signal by a 
known amount. If the two samples are not of nearly equal amplitudes, it is 
necessary to attenuate the larger one with an attenuator of known phase 
shift. The latter operation presents difficulties. 

A method of phase measurement has been developed which overcomes 
these difficulties by permitting measurements to be made with samples of 
unequal amplitudes. The method uses the homodyne detection principle 
and operates in the following manner: The output energy of a signal oscil- 
lator is divided into two portions. One portion is applied to a balanced 
modulator where it is modulated by an audio-frequency signal. The sup- 
pressed-carrier, double-sideband signal from the modulator is applied to the 
device to be measured. As before, means are available for sampling the 
signal at both the input and output of the device. The other portion of the 
oscillator power is fed through a calibrated phase shifter and is applied to a 
crystal detector in the manner of a local oscillator in a double-detection 


receiver. The signal samples are then alternately applied to the crystal 
detector where they are demodulated by the action of the homodyne carrier. 


In each case the phase shifter is adjusted so that the audio signal is a mini- 
mum in the detector output. This occurs when the phase of the homodyne 
carrier is in quadrature with the signal sidebands. The difference in phase 
between the two.adjustments of the phase shifter is equal to the phase dif- 
ference between the two samples. 

Figure 1 shows the apparatus used for measuring phase in this manner. 
Radio frequency power from a suitable oscillator is applied to the H-plane 
branch of an hybrid junction! where it divides and emerges in equal portions 

LW. A. Tyrrell, “Hybrid Circuits for Microwaves,” Proc. 1. R. E., Vol 
1294-1306; November 1947. 
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from the two lateral branches. The portion applied to the calibrated 
variable phase shifter at the top of the figure becomes the homodyne carrier. 
The remaining portion is applied to a balanced crystal modulator’ through 
a second variable phase shifter which need not be calibrated. The latter 
was introduced in order that the phase of any modulated power retlected 
due to an imperfect balance in the modulator could be shifted so that it 
would be in quadrature with the homodyne carrier and would, therefore, 
not produce an audible signal in the detector. 

The portion of the power which enters the modulator is modulated by a 
signal derived from an audio-frequency oscillator. The suppressed-carrier, 
double-sideband signal which leaves the modulator is applied, after a certain 
amount of attenuation, to the input of the device to be measured. Probes 
are provided at the input and output of the latter for sampling the signal. 
Provision is made for connecting either probe to a crystal detector of the type 
used for detecting an amplitude-modulated signal. 

The homodyne carrier emerging from the calibrated phase shifter is 
attenuated to a level of about one milliwatt and is applied to the crystal 
detector. The output of the detector is connected to an audio-frequency 
amplifier terminated by a pair of headphones or an output meter. An 
attenuator may be placed between the amplifier and the detector as an aid 
in measuring the magnitude of a transfer impedance. 

The procedure for adjusting the apparatus and measuring phase is as 
follows: 

With both sampling probes disconnected from the detector the variable 
phase shifter between the oscillator and modulator is adjusted until the 
output of the detector is zero. This balances out the effect of any signal 
reflected by the modulator. The input probe is then connected to the 
detector and the calibrated phase shifter is adjusted until the signal disap- 
pears in the audio output. When this occurs the homodyne carrier is in 
quadrature with the signal sidebands, and the resultant signal applied to the 
detector is equivalent to a phase-modulated wave having a low modulation 
index, and consequently is not demodulated by a detector of the type used 
here. 

The input probe is then disconnected from the detector and the output 
probe connected. The phase shifter is again adjusted for a null in the audio 
output. The difference in phase between the two adjustments of the phase 


shifter is equal to the phase shift between the input and output of the 


device. If the probes are not located exactly at the input and output termi- 
nals of the unknown it may be necessary to make a correction in the meas- 


2 C. F. Edwards, ‘‘Microwave Converters,” Proc. I. R. E., Vol. 35, No. 11, pp. 1181 
1191; November 1947. 
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ured phase by allowing for the known phase shift in the line between the 
probes and the actual terminals of the unknown. 

So much for the general method. Certain precautions are necessary in 
order to avoid errors in measurement. In practice the carrier is not com- 
pletely suppressed in the output of the balanced modulator. It may be at a 
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Fig. 1—Schematic circuit for microwave phase measurement. 
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Fig. 2—Vector diagram of balanced condition with the resultant carrier in quadrature 
with the signal sideband. 


level of the order of 10 to 20 decibels below the sidebands. Since the 
residual carrier will be added to the homodyne carrier in the detector, and 
since the null adjustment will be reached when the resultant carrier is in 
quadrature with the sidebands, it is desirable that the residual carrier be low 
in level compared with the homodyne carrier. The error in phase Ad 
introduced by the residual carrier is shown in the vector diagram of Fig. 2. 
A difference in level of about 40 decibels between the homodyne and residual 
carriers will give an error of not more than half a degree in phase. The 
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homodyne method of detection has all the conversion efficiency of the usual 
double-detection arrangements and, in addition, has the advantage in this 
particular application of having a very low noise level due to the relatively 
narrow band required for the audio signals. The 40-decibel level difference 
mentioned above is accordingly not a serious handicap. 

Other precautions must be observed. The homodyne carrier can be 
brought in quadrature with the signal for two different phases 180° apart. 
This is illustrated in Fig. 2. In many applications, where only the variation 











Fig. 3—Variable phase shifter using a polystyrene vane. 


in phase difference is of importance, this uncertainty of 180° can be ignored. 
The correct setting of the homodyne carrier phase can, however, be deter- 
mined very easily. Assume that the input probe is connected to the receiver 
and that the phase has been adjusted for a balance. Then disconnect the 
audio frequency drive from one of the crystals in the balanced modulator. 
The residual carrier will now no longer be suppressed and the error angle 
Ag of Fig. 2 will become larger. Whether the homodyne carrier is lagging or 
leading the signal carrier can be determined by observing whether more or 
less phase shift, respectively, must be introduced to restore balance. A 
similar test performed with the output probe will indicate whether or not it 
is necessary to add 180° to the measured phase difference. If either probe 
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test indicates a lead, whereas the other probe indicates a lag, then the addi- 
tion of 180° is indicated. 

In microwave circuits it frequently happens that the transfer phase 
varies quite rapidly with the frequency, particularly if some part of the circuit 
is at or near resonance. In measuring the phase characteristics of a circuit 
of this type over a band of frequencies it is necessary, therefore, to take the 
points of measurement close enough together to avoid phase errors corres- 
ponding to multiples of 360°. 

When a balance has been established so that the signal is minimized in the 
detector output, one may observe the presence of the second harmonic of the 
audio tone. This harmonic is a distortion term generated in the detector. 
If it is objectionable, it can be eliminated either by a low-pass filter in the 


audio output, or by using a balanced detector. 


In measuring transfer impedances it is desirable to know the ratio of the 
magnitudes of an output voltage and an input voltage as well as the phase 
difference. The equipment described here can be used for measuring ampli- 
tudes by adjusting the phase shifter for a maximum signal in the audio 
output. Maximum signal levels can then be compared with the aid of an 
audio-frequency attenuator and output meter connected as shown in Fig. 1. 

The apparatus was assembled with standard 4000-megacycle waveguide 
components. A satisfactory phase shifter was made of an ordinary vane- 
type variable attenuator by replacing the resistance strip with a vane of 
quarter-inch thick polystyrene six inches in length. This phase shifter gave 
a total shift of about 100°. Constructional details of this phase shifter are 
shown in Fig. 3. Other phase shifters could have been used with equally 
satisfactory results. It is desirable, however, that the phase shifter be 
impedance matched to the line in which it is located in order that reaction 
back on the oscillator shall be a minimum. In the shifter of Fig. 3 the ends 
of the polystyrene vane have been tapered two inches at each end to accom- 
plish this result. 

The phase shifter can be readily calibrated by using a standing wave 
detector fitted with a sliding probe as a standard of phase. The standing 
wave detector is terminated on one end and connected to the modulated 
signal source on the other. The signal picked up by the sliding probe is 
applied to the crystal detector. Knowing the guide wavelength in the 
standing wave detector, known phase shifts can be introduced by sliding the 
probe along the guide. By adjusting the phase shifter in the homodyne 
carrier path for balance, calibration points can be established. 

The measuring procedure described above has been tested experimentally 
at 4000 megacycles with very satisfactory results. With ordinary care it was 
possible to measure phase differences with an accuracy of better than half a 
degree. 





Reflection from Corners in Rectangular Wave Guides— 
Conformal Transformation* 


By S. O. RICE 


A conformal transformation method is used to obtain approximate expressions 
for the reflection coefficients of sharp corners in rectangular wave guides. The 
transformation carries the bent guide over into a straight guide filled with a non- 
uniform medium. The reflection coefficient of the transformed system can be 
expressed in terms of the solution of an integral equation which may be solved 
approximately by successive substitutions. When the corner angle is small and 
the corner is not truncated the required integrations may be performed and an 
explicit expression obtained for the reflection coefficient. Although applied here 
only to corners, the method has an additional interest in that it is applicable to 
other types of irregularities in rectangular wave guides. 


INTRODUCTION 


HE propagation of electromagnetic waves around a rectangular corner 

has been studied in two recent papers, one by Poritsky and Blewett! 
and the other by Miles*. Poritsky and Blewett make use of Schwarz’ 
“alternating procedure” in which a sequence of approximations is obtained 
by going back and forth between two overlapping regions. Miles derives 
an equivalent circuit by using solutions of the wave equation in rectangular 
coordinates. Several papers giving experimental results have been pub- 
lished. Of these, we mention one due to Elson’ who gives values of reflection 
coefficients for various types of corners. 

Here we shall deal with the more general type of corner shown in Fig. 1 
by transforming, conformally, the bent guide (in which the propagation 
“constant” of the dielectric is constant) into a straight guide in which the 
propagation “constant” is a function of position—its greatest deviation 
from the original value being in the vicinity of points corresponding to the 
corner. This type of corner has been chosen for our example because it 
possesses a number of features common to problems which may be treated 
by the transformation method. 

The essentials of the procedure used are due to Routh* who studied 
the vibration of a membrane of irregular shape by transforming it into a 
rectangle. After the transformation the density (analogous to the propaga- 
tion constant in the guide) was no longer constant but this disadvantage 
was more than offset by the simplification in shape. 

Until this paper was presented at the Symposium I was unaware of any 

* Presented at the Second Symposium on Applied Mathematics, Cambridge, Mass. , 


July 29, 1948. 
1 See list of references at end of paper. 
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other wave guide work based on conformal transformations (as described 
above) except that of Krasnooshkin®. At the meeting I learned that 
the transformation method had also been discovered (but not yet published) 
by Levine and by Piloty independently of each other. Levine has studied 
the same corner, see Fig. 1, as is done here. However, his method of 
approach is quite different in that he obtains expressions ior the elements 
in the equivalent pi network representing the corner, whereas here the 
reflection coefficient is considered directly. This is discussed in more 
detail at the beginning of Section 6. Piloty’s work is closely related to the 
material presented in a companion paper’ and is discussed in its introduction. 

In this paper the partial differential equation resulting from the trans- 
formation, together with the boundary conditions, is converted into a rather 
complicated integral equation. Numerical work indicates that satisfactory 
values of the reflection coefficient, in which we are primarily interested, 
may be obtained by solving this integral equation by the method of succes- 
sive substitutions. However, the question of convergence is not investigated. 

Although they are here applied only to corners, the equations of Sections 3, 
4 and 5 are quite general. In order to test their generality they were used 
to check the expression’ for the reflection coefficient of a gentle circular 
bend in a rectangular wave guide, E being in the plane of the bend. The 
work has been omitted because of its length. It was found that the essential 
parts of the transformation may be obtained by regarding the inner and 
outer walls of the guide system as the two plates of a condenser, solving the 
corresponding electrostatic problem (using series of the Fourier type), and 
utilizing the relation between two-dimensional potentials and the theory of 
conformal mapping. 

When the angle of the corner is small we may obtain the series (7-5) 
and (7-11) for the reflection coefficients corresponding to simple (i.e. not 
truncated) EK and H corners, respectively (a corner having the electric 
intensity £ in the plane of the bend will be called an £ corner or an electric 
corner. H corners are defined in a similar manner). When the angle of 
the general E corner shown in Fig. 1 is small we may use the series (7-18). 

The series (7-5) and (7-11) giving the reflection from small angle corners 
are related to the series giving the reflection coefficients for gentle circular 
bends. In fact, if the radii of curvature of the latter be held constant 
while the angle of bend is made small, the series for the circular bends 
reduce to those for the corners. 

As for the limitations of the method, note first that it can be used only 
for wave guide systems in which the dimension normal to the plane of 
transformation is constant throughout. Moreover, the integral equations 
of the present paper, except for the work of Appendix III, are derived 
on the assumption that the dimensions of the guide approach constant 
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values at minus infinity and the same values at plus infinity. When this 
assumption is not met, a conformal transformation may still be used to 
carry the system into a straight guide. However, there appears to be some 
doubt as to the best way of dealing with the resulting partial differentia! 
equation. One method, discussed in the companion paper’, leads to an 
infinite set of ordinary linear differential equations of the second order. 
Again, possibly the Green’s functions appearing in Sections 3 and 5 may be 
replaced by suitable approximations. 





















1. Representation of Field for Corner or Bend in Rectangular Guide 


Quite often waves in rectangular wave guides are classed as ‘transverse 
electric” or “transverse magnetic”. However, for our purposes it is 
more convenient to class them as “electrically oriented” or “‘magnetically 
oriented” waves.” Thus, the electric and magnetic intensities are obtained 
by multiplying 
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by e““‘ and taking the real part. Here w, u, and ¢ are the radian frequency, 
the permeability of the medium filling the guide (u = 1.257 X 10° henries 
per meter for air), and the dielectric constant of the same (€ = 8.854 & 10-" 
farads per meter for air), respectively. x, y, and ¢ constitute a right-handed 
set of rectangular coordinates in which the ¢ axis is normal to the plane of 
the bend. Equations (1-1) may be verified by substituting them in 
Maxwell’s equations. 
The potentials A and B satisfy the wave equation 
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where Ao is the wave length in free space corresponding to the radian 
frequency w. 

When the electric vector lies in the plane of the bend, as shown in Fig. 1, 
and the incident wave contains only the dominant mode we set 




















A = 0, B = Qsin (r¢/a) (1-3) 
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where a is the wide dimension of the rectangular cross-section, the guide walls 
normal to the ¢ axis are at ¢ = O and ¢ = a, and Q isa function of « and y 
such that 
vO , FQ 
: (1-4) 
Tio = i2@do (1 — Apa */4).” 
The guide walls are assumed to be perfect conductors and hence the tan- 
gential component of E must vanish at the walls. This requires the normal 
derivative of Q to vanish at those walls which are perpendicular to the 
plane of the bend: 
aQ 
—< = Q, 1-5 
on Os) 
When the magnetic vector lies in the plane of the bend and the incident 
wave consists of the dominant mode, we set 


A = P, B=0 (1-6) 
where P is a function of 1 and y such that 


oP . #P 
+ mies 


axt t 9y2 TaP = ¢, = (1-7) 


and 
P=0 (1-8) 


at the walls perpendicular to the plane of the bend. In this case the guide 
walls parallel to the plane of the bend are at ¢ = O and ¢ = 0. 


2. Electric Vector in Plane of Bend 


Figure 1 shows a section of the bend taken parallel to the electric vector. 
b is the narrow dimension of the guide. Let the frequency and the wide 
dimension a of the guide (measured normal to the plane of Fig. 1) be such 
that only the dominant mode is freely propagated. The position of any 
point in this section is specified by the complex number z = x + zy where 
the origin and the orientation of the axes have been chosen somewhat 
arbitrarily. 
The constant & and related propagation constants which appear in the 
formulas dealing with Q and electric bends are given by 
k = (2b/do) (1 — (Ao/2a)?]'" = —iDob/x 
2 


Yn = m — k?*; man. O, 1, 2, ++ ; Yo = tk 


Ao = free space wavelength 
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Since, by assumption, only the dominant mode is freely propagated, k and 
Ym for m > O are real and positive. 

We imagine an incident wave of unit amplitude coming down from 25 
in the upper left portion of Fig. 1. What are the amplitudes of the reflected 
wave traveling back toward z; and the transmitted wave traveling outward 
to the right towards z3? Our task is to find a Q(x, y), satisfying the wave 
equation (1-4) and the boundary condition (1-5), which represents a 
disturbance of the assumed type. 


Z5=0 et (1-2a) 





Mt - ----- 0 -----> 








Ze 
V=t 
Fic. 2 


Tt 
T 
y 





The first step is to find the conformal transformation 
s=xt+iy = f(vt+ i) = fw) (2-2) 
which carries the bent guide (shown in Fig. 1) in the (x, y) plane over into 
the straight guide (shown in Fig. 2) in the (v, 6) plane. This may be done 
by the Schwarz-Christoffel method discussed in Appendix I. This trans- 
formation carries the wave equation (1-4) and the boundary condition 

(1-5) into 
a0 


agp + (LL + ale, RQ = 0, (2-3) 


oC = Oat = Oando =x (2-4) 
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where the upper and lower guide walls are carried into 9 = 0 and @ = -, 

respectively, and g(v, 8) is given by 

1+ g(v, 0) = |f'(v + 10) Po /0 (2-5) 
[ch v + cos 6]°* 

[ch(v — 1) — cos 6]*{ch(v + 1) — cos 6] 

Here ch denotes the hyperbolic cosine, f’(v + 10) denotes the first deriva- 
tive of f(w), and from Appendix I, 27a is the total angle of the bend. /isa 
parameter which depends upon o and the ratio d/dy where d = | 24 — zo | 
and do = | 24 — 2|in Fig. 1. A table giving values of ¢ for a 90° bend 
(a == 1/4) appears in Appendix I. 

That the propagation constant is no longer uniform in the transformed 
guide shows up through the fact that the coefficient of k°Q in (2-3) is now a 
function of the coordinates (v, @). g(v, @) measures the deviation of the 
propagation constant from its value at v = —*. For example, if we 
consider a wave front coming down from 25 we expect it to get past zs before 
it reaches 29. In Fig. 2 the same wave front is tilted forward corresponding 
to a high phase-velocity (or small propagation constant) at z, where »v = 0 
and 6 = 7. This is in line with the fact that the coefficient of &?Q in (2-3) 
vanishes at 24 by virtue of (2-6). Similar considerations hold at 2; and 2. 

What is our reflection problem in terms of the transformed guide? In 


g(v, 6) = —1 (2-46) 


addition to satisfying the two equations (2-3) and (2-4) Q must behave 
properly at infinity. For large negative values of », Q must represent an 
incident wave plus a reflected wave. The incident wave is of unit amplitude 
and the reflected wave is of the, as yet, unknown value Re. For large 
positive values of v Q must represent an outgoing wave. Thus Q must 
also satisfy the two equations 


Q = en ike + Rge***, (2-7) 
Q = Tre-* , ww (2-8) 


where the subscript £ appears on the “reflection coefficient” Rg and the 
“transmission coefficient” T to indicate that here we are dealing with an 
electric corner. 

Our problem is now to take the four equations (2-3, 4, 7, 8) and somehow 
or other obtain the value of Rg. We are not so much interested in T, 
because it does not have the practical importance of the reflection coefficient. 
There are at least two different ways we may proceed from here. One 
is to transform the differential equation plus the boundary conditions 
into an integral equation which may be solved approximately by iteration. 
Another way is to assume Q to be a Fourier cosine series in # whose co- 
efficients are functions of v. Substitution of the assumed series in the 
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differential equation (2-3) gives rise to a set of ordinary differential equations 
having v as the independent variable and the coefficients as the dependent 
variables. The integral equation method is used in this paper. The 
second method is discussed in the companion paper. 


3. Conversion of Differential Equation into an Integral Equation 


The differential equation (2-3) may be converted into an integral equation 
by using the appropriate Green’s function in the conventional manner. 
The only modifications necessary are essentially those given by Poritsky 
and Blewett! in a similar procedure. 

The.conversion is based upon Green’s theorem in the form 


j AG dO 
[(0%8 -« “= 4 _ ds = - flo G — GV’ 0) dv dé (3-1) 


where the integration on the right extends over the rectangular region? < v <2, 
0 < 6 < a (inside the straight guide associated with (v, 9), i.e. the guide of 
Fig. 2) except for a very small circle surrounding the point (v, 40). 
G = G(v, 4 32, 8) is the Green’s function corresponding to 

ch > =  +RV (3-2) 
in the region — x <v< 2,0 < @ < @ subject to the boundary condition 
0V/dn = O on the walls (0V/06 = 0 at 6 = O and 6 = mr). G becomes 
infinite as —log r when r — 0, r being the distance between the variable 
point (v, @) and the fixed point (v , 4). Poritsky and Blewett* have shown 
that, in the notation (2-1), 


a 


, ~1 —|v—vol+¥m 
G = > enyn' cos my cos moe |” 8” (3-3) 


m=0 
q=1,¢6, 22 for m=1,2,3- 


Equation (3-1) leads to 


a: 00 Vr QO 
2O(v% , 80) +[ | -0° +6 eI dd + | {o% es 4 d6 
Ov Ov ve 
e | do | d0 g(v, 0=)OG 
v1 0 


from which the required integral equation for Q is found to be 
O(v0 ) 6) _ 


(3-4) 


a Pa x ‘eo 
¢ wre + | av [ d8 g(x, 6) O(®, 8) Sen Yn cos mB cos mB e °°" (3.5) 
“7 J/—x 


m= 


* We have replaced their i by —i since here we assume the time to enter through the 
factor e'** instead of e *#¢. 





CONFORMAL TRANSFORMATION 111 


where Ym is given by (2-1). The term e~‘**o comes from the first integral 
on the left side of (3-4) as 7; + —*x. Equation (3-5) is a general equation 
which may be applied to a number of wave guide problems by choosing a 
suitable function g(v, 8). For the corner of Fig. 1 g(v, 0) is given by (2-6). 

If g(v, 6) approaches zero when |v, becomes large, as it does for the 
corner, expressions for the reflection coefficient Ry and the amplitude Ty 
of the transmitted wave may be obtained by letting 7% — + in (3-5). 
For very large values of | v9 | the contributions of all the terms in the summa- 
tion except the first (m = 0) vanish. Comparison of the resulting expression 
for Q(v , 00) with the limiting forms (2-7) and (2-8) defining Rz and Ty gives 


Rg = — = | av | d6 g(v, 8) Or, 6)e * (3-6) 
2r — 90 0 


Ty, =1- =} ac | dO g(v, 8) Q(z, ae” (3-7) 
Qn — 00 0 


Since the integrands involve the as yet unknown ((v, @) these expressions 
are not immediately applicable. In fact, if we knew Q(v, @) it would not be 
necessary to use these integrals for Ry and T,—we could simply let »— + « 
and use (2-7) and (2-8). Nevertheless, (3-6) and (3-7) are useful in obtain- 
ing approximations to Ry and J, when approximations to Q are known. 

In Appendix IV it is shown that Rg is the stationary value, with respect to 
variations of the function Q, of an expression made up of integrals containing 
Q in their integrands. From the integral equation it follows that when 
k — 0), i.e., when the frequency decreases toward the cut-off frequency of the 
dominant mode, Q becomes approximately exp (—ikv). Furthermore, Re 
approaches zero. This is in contrast to the apparent behavior of Ry which, 
according to the discussion given in Section 5, may possibly approach —1 
under the same circumstances. Thus reflections from the two types of 
corners, or more generally, irregularities in the £ plane and in the H plane, 
appear to behave quite differently as the cut-off frequency is approached. 


Ry and T, are not independent. Since the energy in the incident wave is 
equal to the sum of the energies in the reflected and transmitted waves we 
expect 


RG + TrT> = 1, (3-8) 


where the asterisk denotes the conjugate complex quantity. In addition, 
there is a relation between RX» and Ty which for a symmetrical irregularity, 
i.e. for g(v, @) an even function of v, states that the phase of Re differs from 
that Ty by +7/2. In this special case T, is determined to within a plus or 
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minus sign when Regis given. These relations may be proved by substituting 
various solutions of equation (2-3) for Q and ( in the equation 


a0 -00] _[,00_ p00 
E ? e| — E av 0 el (3-9) 


ry 

where 2; and 22 are large enough (2; negative and v2 positive) to ensure that 
Q and OC have reduced to exponential functions of ». Equation (3-9) follows 
from Green’s theorem. When Q is taken to be the solution for which (2-7) 
and (2-8) holds and 0 its conjugate complex Q*, equation (3-8) isobtained. 
Keeping the same solution for Q but now letting O denote the solution 
corresponding to an incident wave of unit amplitude coming in from the 
right: 


se | 


O, = ef 4 Rye te, 
QO; T,e**" ‘ 7— — @ 


gives T = T, where we have dropped the subscript £ and have assumed 
that g(v, 6) may be unsymmetrical. Taking ( to be or gives 


RT + iT =0 


which is the relation sought. In the symmetrical case R = R1, R/T + R*/T* 
is zero and hence R/T is purely imaginary as was mentioned above. The 
same relations hold for Ry and Ty. ‘These results are special cases of a 
more general result which states that the “scattering matrix” is symmetrical 
and unitary for a lossless junction.'° 


4. Approximate Solution of Integral Equation 


A first approximation to the solution of the integral equation (3-5) is 
obtained when we assume that the non-uniformity of the propagation 
constant has no effect on Q. Thus we put 


Q™(r, 6) = e* (4-1) 


in the integral on the right and obtain an expression for the second approxi- 
mation Q®(v, #), and so on. Here we shall not go beyond Q®)(2, @). 
It is convenient to expand g(v, 6) in a Fourier cosine series 


eo 


g(v, 0) = >. an(v) cos 8 


n=(0 


(4-2) 


a,(v) = = | g(v, 0) cos n6 dé, é = 1;e, = 2, > 0. 





CONFORMAL TRANSFORMATION 113 


The second approximation, obtained by substituting (4-1) in (3-5), may then 
be written as 


0” (09, 0) = P ai aa se > i Yn cos my 
m=0 


; (4-3) 
-| Am(v)e Re !e-voltm ay 


The nth approximation R{" to the reflection coefficient (when the 
electric vector lies in the plane of the bend) is defined in terms of Q™ by 


Limit Q° (0, 6) = & ™" + RG e™ (4-4) 


v-+—00 


R{” is also equal to the integral obtained by replacing Q in (3-6) by Q®, 
We have 


A) 

1) 2 sian —2ik 

” = 0, ©) — —jk2 | aj(rv)e dv, 
— 


@ = RY — ik? Y (4ymen) 
m=( 


«© (-) ; 
: | dt Am (vo) [ dd dm(v)e Rete) !9—v0l 7m 
ae - 


where Ym is given by (2-1) 
The results of this section have the same generality as the integral 
equation (3-5) in that they are not restricted to corners. 


5. Truncated Corner—Magnetic Vector in Plane of Bend 


When the magnetic vector lies in the plane of the bend the reflection 
may be calculated by a similar procedure. The wide dimension a of the 
wave guide now replaces the narrow dimension 6 in Fig. 1. We shall call 
the result of making this change the ‘“‘modified Fig. 1”. We again assume 
the frequency to be such that only the dominant mode is propagated without 
attenuation. In place of equations (1-3, 4, 5) involving Q we have those of 
(1-6, 7, 8) involving P. 

The conformal transformation which carries the modified Fig. 1 into 
Fig. 2 leads to 

2 2n 
ot 
dv" 06? 


P=0 at @6=0 and @=f 


+ [1 + g@, O° P = 





114 BELL SYSTEM TECHNICAL JOURNAL 


where 
xk = 2a/Xo = —il'oa/m, ; (x2 — 1)'2 = ak/b 
5. = m—1-—e= 
Ao = free space wave length, 
and 
fmoa(? + 10) \??/a? = 1 + g(v, 8). (5-3) 


Here fmoa(w) pertains to the modified Fig. 1. Since the expression for 
f’(w) given in Appendix I is proportional to 6 and since the modified trans- 
formation contains a in place of 8, it follows that g(v, 0) for the magnetic 
corner is exactly the same function, given by (2-6), as for the electric corner. 

It is again assumed that the incident wave coming down from the left 
in the modified Fig. 1 is of unit amplitude and of the dominant mode. 
At large distances from the corner 


>= [fe “* + Ree’) sin 8, 
re (5-4) 
P= Tee sin @, >+ x 
which serve to define the coefficients of reflection and transmission. The 
subscript 7 on the reflection and transmission coefficients indicate that 
here we are dealing with a magnetic corner. 
The conversion of the differential equation into the integral equation now 
employs the Green’s function 


G = 2 > 6, sin m6 sin moe '"~ °° 
m=) 
which corresponds to 
oY. @F , » 
; +«£V =0 
av" + 00? 
V=0 at 6=0 andé 


The integral equation for P is found to be 


P(t, 0) = € °"* sin O% 
(5-6) 


2 a +0 n 2 
kK Aewl _¢ * t 01 On 
+ | dv I dO g(0,0)P(v,0) >. 26,' sin mB sin mBe |" °°! 
a) 0 


Qn m=] 
where the parameters are given by (5-2). This is a general equation. 
For the corner of the modified Fig. 1 g(v, @) is given by (2-6). 
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By letting 1%) —>* — ~ we obtain the exact expression 
ix a] r : ogi s - 
Ry = — — dv | dO g(v, 0)P(2, Ae sin 8 (5-7) 

Wwe % 0 


When dealing with the electric corner we saw that Re — 0 as k > 0. 
The presence of c in the denominator of (5-7) suggests the possibility that 
Ry— —1asc—0. For Ry must remain finite and this may perhapscome 
about through P(v, 6) — 0 in the region, say around v = 0, where g(v, 4) is 
appreciably different from zero. This and the fact that P(r, @) must 
contain a unit incident wave suggest that for » < 0 the dominant portion of 
P(v, @) is 2i sin cv which gives Ry = —1. Incidentally, it is apparent that 
the approximations for P(v, 6) given below in (5-8) and (5-10) (and therefore 
also the approximations (5-11) for Ry) fail when ¢ becomes small. 

The first approximation to the solution of the integral equation (5-6) is 


P(x, 6) = e—" sin 0 (5-8) 


When we introduce the coefficients 


b,(v) = | g(v, 8) sin 6 sin nO dé 
0 


Ald 


sin 6g(v, 6) = >. b,(v) sin n0 


n=) 
bi(v) = ao(v) — ae(v)/2, ba(v) = [dn—i(v) — adnyi(v)]/2, n> 1 


we find that the second approximation is 


4 
2) ce . 2 1 el _s ‘ 
P (29 , 0) = € "® sin 0 + «2 >, 6,' sin mb, 
m=) 


5- 
" (5-10) 
f bn(ee “" a 
— 3 
‘The successive approximations to the reflection coefficient are 
F 2 pto 
1) 2 K Liev 
>» = 0, = | dv b,(v)e “” 
2c — 0 
x a + 
3) 2) - 4 1 j . 
@ — R® — it Yo (4c8,) | diy bya (0) (5-11) 
m=) — © 


»+@ 


| dv bm(vo)e 
— 2 
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6. Series for R® When Corner Has No Truncation 

The integrals which appear in the approximations for the reflection 
coefficients are difficult to evaluate in general. This section serves to put 
on record several expressions which have been obtained for R® when the 
corner is not truncated. Corresponding evaluations of R‘ would be 
welcome since the work of Section 7 for small angle corners indicates that 
R® — R® is of the same order as R®. However, I have been unable to 
go much beyond the results shown here. 

As mentioned in the introduction, H. Levine has studied the effect of a 
corner in a wave guide by representing it as an equivalent pi network having 
an inductance for the series element and two equal condensers for the shunt 
elements. Early in 1947 he derived the following expressions (in our 
notation) for the elements corresponding to a simpe & corner :* 


i[¥(P51)- 9-9] 


(kr) cot (82/2) 


Ba/Vo 


Be/ Yo 


il 


where Yo is the characteristic admittance of the straight guide, iB, the 
admittance of one of the two equal shunt condensers, —7B, the admittance 
of the series inductance, ¥(«) the logarithmic derivative of I(x + 1), and 
Br is the total angle of the simple corner (for no truncation we set 8 = 2a). 

When the reflection coefficient for the corner is computed from the 
equivalent network for the case 8 — 0 it is found to lie between the approxi- 
mate value R}? given by (7-3) and the considerably more accurate value 
RY given by (7-5). All three approximations are of the form Ag? + O(6°) 
where A differs from approximation to approximation but is independent of 
B, and O(8°) denotes correction terms of order 8°. Since rR? gives the exact 
value of A, it may be regarded as the standard when the three approxima- 
tions are compared. If this comparison be taken as a guide, it suggests 
that the rather cumbersome expressions (6-2) and (6-5) for R¢” given below 
are not as accurate as the simpler expressions resulting from Levine’s work. 
Dr. Levine has also obtained corresponding results for the general /-corner 
of Fig. 1. It is hoped that his work will be published soon. 

When the corner is not truncated it is convenient, as mentioned above, 
to replace 2a by # so that 7m is the total angle of the bend. For no trunca- 
tion ¢ = 0 and (2-6) becomes 


s 8 
ect we E: + cos 4 wk. (6-1) 


chy — cos@ 


* I am indebted to Dr. Levine for communicating these expressions to me. 








i 
| 
} 
{ 


% 4 
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. From (4-2) and (4-5), or from (3-6), 
’ : —2ike 
‘ P=] «ae | acl, He ” 
(6-2) 


cos 0/ch v 


ala + 1)- 











where we have expanded g(2, 
and integrated termwise. 


For a right angle corner B = 
may be obtained by subtracting the sum of the series corresponding to 
k = 0, namely 


Tin — ik)T (n+ ik) < 
~ al(n — 1)!2— 


(- — 28) 2m (8) n—m 
Axo (2m)!(n — m)! 





—ik r= 


n=] 


6) as given by (6-1) in powers of 
The notation is (a)o = 1, @)n 


-(a+n-— 1). 


1/2, and a more rapidly convergent series 





32-3 U2), (6-3) 
a= 520 
Thus for 8 = 1/2 
43) = —tk [ozs 2 = ¥ ¢ o Bey — 40], 
n=1 
we (6-4) 
A, = rk/sinh zk, A,=A, [] (+m), n>1 


m=1 


The rate of convergence of the more general series (6-2) may be increased 
in a somewhat similar way. It is found that 





: = -8 [a — 26°(1 — Aj) farsa — As) 
-% g (23 + 208° + 28')(1 — As) - | 
J=K+L (6-5) 
ee @.. 1 = bs a o 
K = a nin 1-8 v1 B) 5772 
pa (—28) 2m . (S).- es — = (1/2 sage 8) m : 
aia A= (2m)! sou (n— m)\n Bp > (1/2)mm(m — B) 
where .5772 --- is Euler’s constant, V(x) is the logarithmic derivative of 


II(x) = I'(x + 1), and A, is given by (6-4). 
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wry ° (2) . . . y 
[he results corresponding to Ry, are quite similar. When the corner 
is not truncated 


i) 


C— a 2 e® T(x —_ ic) (n + ic) : (—28)2m(B)n m 
lil " (a + 1)!(2 — 1)!2c m=p (2m) !(n — m)! 


; E — (1 — 4)) - . (2+ 6)(1 — A) (6-6) 








lk 


ll 


<¢ 2 
e (23 + 208° + 28°)(1 — 43) — --- 
270 : 


jul +e «SR 


~ (1 .~ 3) m 
— BU — 8) > (1/2),,m(m — B)(m —- B+ 1) 








in which A, is obtained by replacing ¢ by & in the expression (6-4) for A,. 

The evaluation of the integrals for R¢? and Rf’ for general values of / ’ | 
appears to be difficult although it is possible to obtain approximate expres- 
sions for the case when / is large. 4 


7. Reflection from Small Angle Corners 


The expressions for R® and R® may be evaluated approximately when 
the angle of the corner is small. It turns out that, for ¢ = 0, they are of the 
same order of magnitude and both of them must be considered. Moreover 
R™ for n > 3 differs from R® by terms of the same order as those neglected 
in our approximations so that there is no point in going to the higher values | 
of m. 4 

We first obtain the approximation for Re for a corner with no truncation : 
having the total angle 78. Since 2 is very small (6-1) may be written as 

g(v, 6) = exp [By] — 1 = Be + B¥y"/2! + 0(6*) 


io ij 
(7-1) ; 
log (chv + cos 0) — log (chv — cos @) 


I| 


¢ 


where 0(3*) denotes terms of order 3%. ‘The expression ¢ becomes very large 
near the two points (0, 0) and (0, a) (the coordinates being (2, @)). The +t 





following considerations indicate that this does not invalidate our procedure. 
The remainder, denoted by 0(6%), in (7-1) is less than | By |? exp | By |. 
Near (0, 0) ¢ is approximately equal to 2log(2/r) where r = 2? + @&. 
Consequently the remainder is less than (28 log 2/r)3 (2/r)”’. When the 
expression (7-1) for g(v, @) is set in the integral equation it is seen that all 
terms, and in particular the remainder term (by virtue of the inequality 
just stated), of the double integral converge at (0,0). Hence the contribu- . 


tion of the remainder term is of order 8° even in the worst case when the 
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Green’s function is replaced by —log r._ A similar result holds for the other 
point in question, namely (0, 7). 

Integrating (7-1) from @ = 0 to @ = m and using equations (A2-1, 3) of 
Appendix IT gives 


ag(v) = B{l,(v, v) — To(v, v)] + 0(8%) 


=48 dD ne + 0(63) (7-2) 
n=1,3,5,--- 
where v > 0. We consider only positive values of v since g(v, 8) and the 
a,(v)’s are even functions of v. Thus (4-5) yields 
RP = —ik2e D> n(n? +B)" + 068’) (7-3) 
n=1,3.5,--- 
This is an approximation to the exact value given by the double series in 
(6-2). Comparison of (6-5) and (7-3) when 8 and & approach zero gives, 
incidentally, 


md n/t = 1S 


From (4-2), (7-1) and the expansions (A2-2) of log (chv + cos @) it follows 
that 


Amn(v) = 48m e—™l"l + 0(8?), m= 1, 3, 5, 
(7-4) 
dm(v) = 0(B?) , m= 0(0,2,4,6--- 


Equations (4-5), (A2-4), the relation y,, = m? — k?, and (A2-8) give us the 
answer we seck: 


R® = RY — ik*2¢’ 2. Yn_m” J(m, m, k, Ym, 0,0) + 0(8°) 


m=1,3,5,+ ++ 


—ik2e >> nim? + 0(6°) (7-5) 
m=1,3,5,-++ 
It is not necessary to go to Ry? because it differs from Rf” by only 0(8°). 
When // lies in the plane of the bend the reflection from a small angle 
corner with no truncation may be obtained by much the same procedure. 
For brevity we shall not write down the order of magnitude of the remainder 
terms. From (5-9), (A2-1), and (A2-3) 


bi(v) = ao(v) — ao(v) /2 (7-6) 
Bh — Iz — (Is — 1h)/2) 


I 


> 


= fir +e SF ate = Pe | i teas 


n=3,5,--- m=2,4,--« 
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where we have written J, for J(v, v) and assumed v > 0. Then, using 
(5-11), 


O = Hoek t+2 Dwi +e)7 
n=3.5.-+ 


2 -17.3 , 2-1 (7-7) 
—-2 Dd ale’ -—1)°(n"' +2)" 
n=2,4,+++ 
When we put 
b,(v): = [an—i(v) — dngi(v)]/2, a> (7-8) 
b,(v) = 28[(m — 1)Me-@—Y ll — (nn + 1) Me“ FIl], n = 2,4, 6, --- 
ba(v) = 0(6"), ilies i, 3, 5, tite 
in (5-11) and use the results of Appendix II we obtain 
RY = RP — ‘ce 5, [(n — 1)? J (n — 1,” — 1,8,, 0,0) 
n=2,4.6--- 
+ (n+ 1)? J(n +1, + 1, c, 5, 0, 0) (7-9) 


— 2(n’ — 1) *J(n — 1, + 1, c, bn, 0, 0)] 


The values of the first two J’s, obtained by setting m = n + 1 in (A2-7), 
may be simplified by using 


e+ (m+1+4+ 6)? = 2n+ 1)(n+4+ 38) 


where we have dropped the subcript » from 6,. In order to eliminate 6 


from the denominator we multiply both numerator and denominator by 
n — 6 and use 


(n — 5) (6 + 2n + 2) 


(nw + 13? + c? — 6(m + 2) 
wW—-Ff=-14+c= 


Setting in the value, given by (A2-9), of the last J and separating the 
terms (into those which contain the first power of 6 and those which do not) 
enable us to write the term within the square brackets in (7-9) as 


no “(I (n —1) —1 rf (n+ 1) +1_ 
(in? — 1)? — PL (mn — 1)*{c? + (nm — 1)7} (mn + 1)*{c? + (n + 1)?} 


2n{2(n' + ¢) — «(mW — oY 
+ x(n? — 1)°(n? + c?*) ; ies 


It is found that when (7-10) is put in (7-9), the contribution of the first 
two terms within the square bracket of (7-10) exactly cancels the summation 


i 
i 
| 
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which is taken over 3, 5,7 --- in the expression (7-7) for RY. Moreover, 
if we make use of 
> 4n(n® — 1)? = 1 
n=2,4,6,°+- 
we see that the contribution of the last term within the square brackets of 
(7-10) cancels the remaining terms in a. Only the contribution of the 
first term in (7-10) remains and it gives 
@ = ee" YO win — 1) 782 + 06) (7-11) 
n=2,4.6,--- 

The relative simplicity of this result indicates that there may be another 
method of derivation which avoids the lengthy algebra of our method. 

Recently approximate expressions for the reflection coefficient of gentle 
circular bends have been published’. In our present notation these may be 
written as 


. — k 
2 —2| Sin u cos u — e “™! 
Ry = — ib p; | — — 4k > — 
24 manl,3.5.--° * m' Ym 
. —us 2 
R ig? pv? | SB ™ cos u — e “énle n 
n= — 1a py ~ ~ e 7 | 
81? c? n=2,4.6,--: sO (n? — 1)8 


where {7 is the angle of the bend, p; is the radius of curvature of the center 
line of the guide and u is 2m times the length of the center line in the bend 
divided by the wavelength in the guide: 


u = Brrkp,/b = Brcp;/a 


The first expression for u is to be used in Rg and the second in Ry. If we 
now let 8 — 0, keeping p; fixed, then «> 0. The trigonometric and expo- 
nential terms may be approximated by the first few terms in their power 
series expansions, and part of the series which make their appearance may be 
replaced by their sums given, for example, by equations (4.1-7) and(4.1-8) 
of reference’. After some cancellation, the above expression for Rg and Ry , 
which hold for gentle circular bends, reduce to (7-5) and (7-11), respectively, 
which hold for the sharp corners. In other words, the reflection coefficients 
for both the sharp and the circular bends approach zero as 8 — 0, and 
furthermore their ratio approaches unity. 

We shall merely outline the derivation of the approximation Rf? for a 
truncated corner. Instead of (7-1) we have from (2-6), 


g(v, 0) = exp [ag] — 1 = ay + a%g*/2! + O(a’), (7-12) 


yg = 2 log[chv + cos 6] — log [ch(v — #) — cos 6] — log [ch(v + t) — cos 8} 
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The Fourier coefficients of g(v, 8) may be obtained by using the results of 
Appendix II. Assuming v > 0, m > 0, 

ao(v) = 2a(v — tv) + 2a? {47,(v, v) + Mi(v — tho — 0) 

+ [i(0+ thu+t) — Ale(v — t, v) — 4Io(v + ft, 0) + 21,(0 — th 0 + DO}, 
Gm(v) = 2Zam™{—2(—)m™emmlel + emmle—4l emir tl] (7-13) 


where ¥(v) = 1 when 0 < v < ¢and ¥(v) = 0 whenv > ¢. Substitution 
of the values (A2-3) for 7; and J, gives 


ao(v) = [2a(v — t) + 2a2(v — £)*}Y(2) (7-14) 
+ a> we + en 2ne—2nt — 4(—)"en2ne-nt 
n=1 


+ e -2n|v—t| +. de—n|v—t| —njoti| 4(— )re—n|r—t| “nel 
The second approximation to the reflection coefficient is 


R® = iak™ sin®kt — io®k22-1(2kt — sin 2kt) 
— iko® >> n(n? + kb)! {2—(—)"dem (7-15) 
n=1 


+ [1 — 2(—)"e™** + €-™*IJcos 2k 
+ nk“"[e?n! — (—)"2e—"*|sin 2kt! 
The typical term in the summation (4-5) for Rf? is 


a ina 


+00 
~ik(v+v9)—lv—v9| e 
—_- — dv an(to) [ dv dm (ve eet rea (7-16) 
4 mEm —2 —oO 


When m = 0, € = 1, yo = tk, and ao(v) is 2a(v — t) + Oe?) for0O << v <1 
and is O(a*) for v > ¢. The integral may then be approximated by replacing 
the upper limit « in (A2-14) by ¢. The value of (7-16) for m = 0 is found 
to be, to within O(a’), 


2T1e?P(e?*** — 1) — (3/4)ia®k-*(sin 2kt — 2kt) (7-17) 


When m > 0, €m = 2, Vin = m°* — k*, and the substitution of the value 
(7-13) for am(v) enables us to express (7-16) as the sum of six J’s where J is 
defined by (A2-4). The J’s may be evaluated with the help of (A2-7) and 
(A2-8). Substitution of this value of (7-16) and the value (7-17) for m = 0, 
together with RY given by (7-15), in the expression (4-5) for R® gives 
our final result 


RE? = iaksin? kt + o22-(e~2#*t — 4) (7-18) 


+ ia? [4k-*(2kt — sin 2k) — Bsin 2kt + kD, 270 Aal 


n=] 





phWes® 
So 











= 
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where 
oo 
Be z. n-fe-2nt — 2(—)ne-n4] 
n=1 


A, = cos 2kt — [2cos kt — (—)"e~™*f? 


3 


Equation (7-18) is an approximation, to within terms of order a’, for the 
reflection coefficient of a truncated corner which turns through a small 
angle 2ra. The electric vector lies in the plane of the bend. When? = 0, 
(7-18) reduces to (7-5) by virtue of 2a = 8. 


APPENDIX I 
CONFORMAL TRANSFORMATION OF TRUNCATED CORNER 


We shall use a Schwarz-Christoffel transformation* to carry the guide of 
Fig. 1 into the straight guide of Fig. 2. The first step is to transform the 
interior of Fig. 1 into the upper half of an auxiliary complex plane which we 
shall denote by ¢. Let the points 2; , 32, 23, 24, 25 in Fig. 1 correspond to 
the points — h, h, 1, «©, —1 in the ¢ plane. A suitable transformation 
is then 


ig 
s= D+ E[ (r + h)-*(r — h)-(r — 1)7°(r7 +: 1)7dr_—s (A 1-1) 
0 


where D, E and hare to be determined from the geometry of Fig. 1. Because 
of the symmetry of our transformation about the line joining 2 and 2, 
it follows that z = 2) corresponds to¢ = 0. Hence D= 2. As ¢ travels 
from 1 — eto 1+ e, e being very small and positive, along a semicircular 


indentation above ¢ = 1, z as given by (A1-1) increases by 


l+e nm 
Ra + VY "2 5 (r —1) ‘dr = ‘Ear 


ie 
a ; (1 1) 


while, according to Fig. 1, it increases from « + 10 to « + ib. Hence we 
set the real part of E equal to —2bm"(1 — h*)*. We have tacitly assumed 
the factors in (A1-1) to have their principal values at 7 = 1 + € and also 
thatO < ht <1. Asz goes from 2; to 22, ¢ goes from —h to +h. In this 
range arg(r + 4) = O and arg(r — h) = @. 
Consequently, if sz. — s;, = ¢, then 

h 


Zo — 21) = le or = —Ee*™* (h? na T) "i noe t*)—dr 
—, 


*See, for example, S. A. Schelkunoff, Electromagnetic Waves, New York (1943) 
pp. 184-187. 
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and we see that £ is purely real. Hence 


h 


(= Yr — ie fe = 8 = Ade 
—h 


is an equation from which / may be determined as a function of f. Setting 
7? = h’x, expanding (1 — /®x) in powers of #® and integrating termwise 
leads to 


—4 " ea : " . 
. aa (1 — h’)*h'**F(1, 353 — a;h’) 


2b r(3 — a) 
si. & — 9 9 
= = Td = a) h'**F(2 — a, 1 — a;$ — a; Fh’) (A1-2) 
r(} — a) 
= <7 9 9 f 9 
= tp TCO) pty _ yr, 414651 -— #) 
sin Ta (3 — a) 


where we have used relations from the theory of hypergeometric functions. 
The term 1/sin 7a is the reduced form of an original term containing a 
hypergeometric function which has been evaluated by the binomial theorem. 
The second and third expressions are suited to calculation when h? < 1/2 
and A? > 1/2 respectively. 

Now that the guide of Fig. 1 has been transformed into the upper half of 
the ¢ plane, the next step is to transform this upper half into the straight 
guide of Fig. 2. We want ¢ = —1, i.e. 35, to go intov = —» and ¢ = 1, 
i.e. 33, to go intov = +. Again using the Schwarz-Christoffel formula 
with w = » + 76 (the exterior angles at 1 = + are equal to 7) 


i? 
w=D,+ E, [ (7 +1) (r — 1)" dr (A1-3) 


We take the point 2» in Fig. 1 to correspond to v = 0, @ = Oin Fig. 2. Since 
this corresponds to ¢ = 0, D; must be zero. Also dw/df is real because w 
traverses the walls of the guide of Fig. 2 as ¢ moves along the real axis in the 
¢ plane. Hence £; is real. As ¢ goes from 1 — ¢€ to 1 + e€ around a small 
circular indentation above ¢ = 1, w changes from « to « + ir. Thus 


im = E\2-'"(—iw) or £, = —2 (A1-4) 


When (A1-3) is integrated, (A1-4) inserted, and the result solved for ¢ 
we obtain 


Se 


¢ = tanh w/2 (A1-5) 





as 


aitgeicieagiay Wh Lali 
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The function we require is obtained by differentiating (A1-1) and (A1-3): 
dz _ dz /dw 
dw d¢/ d 

= EP — #Y"(P — ITE — 1) 
= (1 — #)°(0 — h’)*b/x 


I'v + 16) = f'(w) = 


(A1-6) 
- | ch'w/2 a 
mw Lsh}(w — t) shd(w + 2) 
az | a 
r (e-* abs 1)(ert# am 1) 
where 
h = tanh t/2 (A1-7) 
For a 90 degree corner a = 1/4 and 
£ = 1 — das) (A1-8) 
26 
where, in Fig. 1, d = | ss — 29| and dy = | 24 — 26]. In order to obtain 





the relation between ¢, defined by (A1-7), and d/dy various values of h” 
were picked and the corresponding values of ¢ and d/dy (using (A1-2) and 
(A1-8)) computed. Representative values are given in the following table. 


d/do l d/do t 
1.000 0 5796 1.2302 
9041 0633 5385 1.4910 
8565 1417 5000 1.7594 
8292 .2007 4615 2.0634 
7745 3500 3127 2.8872 
.7196 5421 .2804 4.0096 
.6919 6549 .1708 5.987 
.6273 .9624 0959 8.294 


APPENDIX II 
INTEGRALS ASSOCIATED WITH CORNERS OF SMALL ANGLE 


The derivation of the integrals encountered in Sections 7 and 8 will be 
outlined here. The first ones are 


I,(u,v) = . | loz(ch u — cos 0) log (ch v — cos @) dé 
T Jo 
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I.(u, 7) = : loz(ch u — cos @) loz (ch » + cos 6) dé 


(A2-1) 
2 
I;(u, v) = cos 26 loz (ch u — cos @) loz (ch v — cos 0) dé 
T 0 
) . 
I,(u, rv) = — [ cos 20 loz (ch u — cos @) log (ch » + cos 0) dé 
T 0 
Assuming “ and ¢ to be positive and using the expansions 
x 
log(ch u — cos 0) = log(e"/2) — 2>, ne"cos nd 
n=1 
(A2-2) 
P<) 
log(ch u + cos 0) = log(e"/2) — ar. (—)"u“"e-"“cos nO 
n=! 
leads to 
2 
I\(u, v) = log(e"/2) log(e’/2) + 2). nem" 
n=] 
a] 
T(u, vy) = log(e“ 2) log(e” 2) + 3h (—)"n~e-"™ nv 
n=l 
I3(u, v) = —e“log(e"/2) — e**log(e"/2) + Je" 
(A2-3) 


20 
+ 2>- Nn lin oa 2) le nu =“ oh e") 


n=1 


T4(u, v) = —e-“log(e”/2) — e**log (e*/2) — 2e-*-* 


o 


2 
+ 2>> (—)"n 1( yy + 2) lp nu nv (p—2u + e~2”) 


n=l 


When w or v are negative they are to be replaced by their absolute values 
in the expressions (A2-2, 3). 
Now we consider the double integral 


a +2 +o 
J(u, (C85 7,-5) = | dio | d» 
> 4] — 2 


7 
-exp [—u|t% — r| —m|v — s| —icle + vo) —é]7 — wf] 


in which pw, m, c, 6 are real and positive and r and s are real. The double 
integral may be reduced to a single integral by substituting 


P -\v-—vol __ i (8° 4 x”) A fate ro) dx, (A2-5) 
T J-x 


AS ea 








eee 
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interchanging the order of integration, and integrating with respect to 2 
and tv. Assuming s — r > 0, the integral is then evaluated by closing the 
path of integration by an infinite semicircle in the upper half plane and 
calculating the residues of the integrand at the poles i, ¢ + im, —c + in: 


ss 4iume*°-?-"*) dy 
Jt Mh. 6. 65:7, $s) = 
2 a( 


(5? + x*)[p? + (x + cia? + (x — 
7 (b+te)8+(b6—te)r 
46 : 
wy: E + (¢ + 16)*][m? + (¢ — 16)? 


ms+(m—2ic)r 


. m|6> + (¢ + im)*][u? + (2c + im)?} 


(A2-6) 


r—(u+2ic)s 
+ eo - ° 9 9 a ° 9 | 
wld? + (c — in)?|[m? + (2c — in)?} 
Substituting special values for the parameters gives the results required 
in the text. Thus, 
J(m, m,k,y; t,t) = e?** I(m, m, k, y; 0, 0) 
JI(m, m, k,y;—t,t) = &** Jim, m, k, y; 0, 20) 
J(m, m, k,y; —t, 0) = &** J(m, m, k, y; 0, 0) (A2-7) 
2m(5 + 2m) 
(ce? + m*)[c? + (m + 6)?} 
which hold irrespective of any relations between the parameters. The 
derivation of the last result is simplified by setting a = c + im, & = ¢ — im 
and factoring the denominators in (A2-6) so as to obtain terms of the 
form & + 16,a@ + 16. 
When 7? = m® — k* considerable simplification is possible and we obtain 


J(m, m, c, 6; 0, 0) 


7 | 1 m | 
I(m k, y;0,0) = — 
J(m, m, k, y;9, Re E mn? + ze 
ikt yt oe k k kt) (A2-8) 
ve a: e m cos kt — k sin Ri 
J(m, m, k; y, 0, 1) = - 
sis eats k* | Y m> + Rk | 
If we put p = n — 1,m = n+ 1, and set & = WP -—-1—-—C=nr7 -K 
where x” = 1 + c, (A2-6) yields, after some reduction, 
n —1 (n — 1)6 
JI(n -- 1 1,¢,6;0,0) = = — — 
” ee het (cn + 16) + 2(n + 1)(1 — ic)*(n — ic) 
+ ye Lge (A2-9) 


a(n — 1)(1 — tc)?(n + ic) 
= cw —s nd[2(n> + c’) - K (ne — 1)] 
~ KA(n? — 1) x(n? — 1)(n? + c) 
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The form of the final expression has been chosen so as to be suited to the 
use we shall make of it. 
Another double integral which appears in our work is 


I(k, y) = [ | dv a(v%) 
ee (A2-10) 
. i] dv a(v) exp [—tk(v + %) — y |v — 20 |] 


where a(v) is an even function of v and is such that all of the integrals 
encountered converge. We begin our transformation by dividing the 
interval of integration (— ~, ©) for vp into (— », 0) and (0, ~). Making 
the change of variable 1 = —1,v = —’ in the first interval, dropping 
the primes and using a(—v) = a(v) leads to 


I(k, y) = 2 | dio a(vo) / dv a(v)e*'” *"! cos R(v + 2%) ~— (A2-11) 


We now split the interval of integration of v in (A2-11) into the intervals 
(— 2, 0), (0, v0), (v0, ©). In (—~, 0) we change the variable from v to 
—v', drop the prime, and use a(—v) = a(v). By paying attention to the 
sign of v — v we may remove the absolute value sign. By changing the 
order of integration in the double integral arising from the third interval 
(in which 0 < v% < *%, % < uv < &) we may show that it is equal to the 
double integral arising from the second interval. Thus 


1(k, y) = 2 | dt a(c) [ dv a(v)e *” *”® cos R(t — 2) 
0 0 
. ” (A2-12) 
co 4 | di a(w) [ dv a(vje 7°" cos R(t + 2) 
0 0 
When a(v), y and & are real we may write (A2-12) as 
I(k, y) = 2 [ dv a(v)e | 
0 
- , (A2-13) 
+ 4 Real | diy a(r)e a dv a(v) e’*** 
0 0 
and when y = 1k we have 
I(k, ik) = 2 | drya(o) | dv a(v)e 
0 0 
(A2-14) 


) ve 
+ 2[ dy a(t) | dv a(v)[e** + e *tr oy 
° ) 


ise ie in 


































CONFORMAL TRANSFORMATION 


APPENDIX III 


INTEGRAL EQUATION WHEN GUIDES ENTERING AND LEAVING IRREGULARITY 
ARE OF DIFFERENT SIZES 


Here we shall indicate how the integral equation method may be extended 
to cover the case mentioned in the above title. It is supposed that only the 
dominant mode is propagated freely in both guides. 


E in Plane of Irregularity 


Let the notation for the guide carrying the incident wave be the same as 
for the E-corner. 6 denotes the narrow dimension of the guide and the 
quantities & and y», are given by (2-1). Both guides have the same wide 


™ 








A 
' 
b ‘INCIDENT WAVE 
We 





A 

! 

' 
by 
! ' 
' ' 
y —— y 

° 
Fic. 3 


dimension a. The narrow dimension of the guide shown on the right of 
Fig. 3 is b}. We introduce the new quantity 


ky = [(2b1/Xo)? — (b1/a)?]'” (A3-1) 


to correspond to k. Since, by assumption, only the dominant mode is 
freely propagated in both guides both & and ; are real positive quantities 
less than unity. 

Let z = f(w) carry the system of Fig. 3 into a straight guide of width 
in the w = v + 76 plane (see Fig. 2), and let g(v, 0) be defined by 


1+ g(v, 6) = | f’(w) ??. 
The behavior of g(v, @) at infinity is shown by the table 


v dz/dw g(v, 8) 
—~« b/ 0 
+a by/m kik? —1 


where },/b = ki/k has been used. It is convenient to introduce the ap- 
proximation ¢(v) to g(v, 6). g(v) may be chosen at our convenience subject 
only to the conditions that it be differentiable, g(— ©) = 0, and g(~#) = 
kik? — 1. 

When we define G by equation (3-3) so that, as before, it is the Green’s 
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function corresponding to a guide of width 6, we may use equation (3-4) 
to derive the new integral equation 


O(% ; Ao) = Pees fe a | ao | dé 
2 — 0 


, \ P Pa —tkyv ’ A3-2 
-{ g(r, AC(r, 0) — g(r) Tre "1G (0, Oo 5 0, 8) \ ) 
+ Tr F (zo) 
in which 
F(a) eel *1°0 § (29) &( i) leah one N (vo) = eo N* (v9) 
y ™ ne 1 = 1 wp l(a \ > i(ky—k)e ? 
N (v9) ys k(k, k) :. 4 (ve di (A3-3) 
N* (to) = 2 "k(ky oe k) | i (ve ares di 
Here ¢’(v) denotes dg(v)/dv. Equation (A3-2) and 
Limit Q(7, 0) = Tre“ (A3-4? 


vo 


are to be solved for the unknown function Q(2, @) and the unknown quantity 
T,. The method of successive approximations may be used in somewhat 
the same fashion as in the simpler case but we shall not give a general 
discussion. 
The first approximations are found to be 
Ty’? = 1/N-(~), RY? = —N+(—@)/N-(#) (A3-5) 
where the .V’s may be obtained by setting 7% = +~* in equations (A3-3). 
One of the simplest choices for g(v) is to let it be zero for negative values 
of v and to have the value (%) = kjk? — 1 for positive values of 7. Then 
Te = 2kkit+ ky, Re = (k-—k)(R+h)" — (A3-6) 
These are quite similar to the corresponding expressions for a transmission 
line which have been used extensively in wave guide work. 
In working with these formulas, when & is small, it is sometimes convenient 
to use the result 


v2 7 v9 rT 
| di | dé g(x, 0) = wb | ao | d0 | f'(w)/) — (22 — um) (A3-7) 

Vi 0 V1 0 
where the evaluation of the double integral on the right is made easier by 


the fact that it represents the area in the original guide (in the (x, y) plane) 
enclosed by the lines corresponding to v = 2, and v = 2. v2 and % are 


— 
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chosen to be moderately large positive and negative numbers, respectively. 
It turns out that, when &; and & are very small, this is related to the “excess 
capacity” localized at the irregularity whose effect must be added to that 
of the mismatch, indicated by (A3-6). 

When the entering and leaving guides are of the same size it is still possible 
to use the formulas of this appendix. -V~(v) may be replaced by an expres- 
sion which now has for its limiting value 


N-(~) = 1+ i(k/2) [ 6(v) de (A3-8) 


IT in Plane of Irregularity 


Let the figure corresponding to the irregularity be Fig. 3 with 6 and 6, 
replaced by a and a,, respectively. In addition to the quantities c and x 
defined by equations (5-2) we define 


' 2 1/2 
K; = 2a,/Xo, Cy = (xi — 1) (A3-9) 
where we assume « and x; to lie between 1 and 2. Atv = —x P(x, @) 
still consists of the unit incident wave plus the reflected wave given by the 


first of equations (5-4) and g(v, @) is still zero. However, now, atv = =x 


b 


P(v, 0) = Tye “* sin 6 


—~-1=«x (4-2) (A3-10) 


” 
O(%2) = KK 


The integral equation for P(v, 6) and Ty is 
9 


P(t , %) =e ° sin 0 + ~ / dv [ 


a J-x “0 
/ me ils . (A3-11) 
dO} g(v, 0) P(v, 0) — g(v)Tye “" sin O}G(z0 ,  ; 2, 4) ” 
+ Ty sin 0) F y(z0) 
in which 
F y(t) a € —s &(%) &( x) —e *: M (1%) = ih M* (v9) 
M (v9) = x’ (2c) ‘(Gy — c) ' | (ze we dv (A3-12) 


M ‘(t) = x (2c) aC + ¢) s i (ie i(eyte)e dv 


v¢ 
First approximations are 


TY oe 1/M~-(), RY = —M+(— 2 )/M~() (A3-13) 
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which, when we choose §(v) to be zero for v < 0 and kik — 1 for v > 0, 
become 


f= 2lato7, Rey =(c—alartec (A3-14) 


which again agrees with results obtained from transmission line considera- 
tions. When the entering and leaving guides are the same size we may use 


M~(o) = 1+ ix°(2c)7 ‘a 6(2) de (A3-15) 


It seems difficult to give any general rules for the choice of g(v). Since 
for Ry and Ty, , the factor sin @ reduces the effect of the singularities on the 
walls of the transformed guide, the choice ¢(v) = g(v, /2) suggests itself. 
The factor sin @ is not present in the formulas for Rg and 7 and regions 
near the walls are more important. In this case the selection 


ian l ” e(v, 6) dé 


may be useful, especially since it allows us to use the result (A3-7) when k 
and k; become small. 


APPENDIX IV 


VARIATIONAL EXPRESSIONS FOR REFLECTION COEFFICIENTS 


The reflection coefficients are proportional to the stationary values of 
certain forms associated with the integral equations. In order to obtain 
these forms we proceed as follows. It is readily seen that the values of 
x, and x2 which satisfy the symmetrical set of equations 


A141 + Ayox%2 = Dy 
(A4-1) 
@j2X; + Ayo%2 = be 
are the ones which make 
je aux of 2410% 1X2 +- AaoX = 2b 4x1 —= Qhoxe (A4-2) 


stationary when x; and x2 are given small arbitrary increments. This 
stationary value of J is 


Li = —byx, = boxe 


If we take the integral equation to be the analogue of the set of linear 
equations, the reflection coefficient turns out to be proportional to J,. 
In order to set down the actual expressions it is convenient to write r for 
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(v, 6) and dS for the element of area dvd@ so that the integral equation 
(3-5) for Q(v, 6) may be written as 


Oro) = @ + R(2e) | g(r)O(r)G(ro , r) dS (A4-3) 


where the integration extends over the interior of the guide and G(r, r) 
denotes the Green’s function (3-3). 

If the number of equations in the set (A4-1) were increased from two to a 
large number .V, the set of x’s would correspond, say, to the values of Q(r) 
or of g(r)Q(r), and the 6’s would correspond to the values of exp(—ik?»). 
In any event, we take the analogue of J to be 


Je = [ momen — 2e“"| dS 
(A4-4) 
— §2(2e)" al g(nO(r)g(ro)Q(ro)G(ro , r) dSo dS 


where the subscript £ indicates that we are dealing with an electric corner. 
It may be verified,* by giving Q(r) a small variation 6Q(r), that the function 
Q(r) which makes J, stationary is the one which satisfies the integral 
equation (A+-3). Furthermore, when we assume Q(r) to satisfy the integral 
equation, the expression for J reduces to an integral which is proportional 
to the integral (3-6) for the reflection coefficient Ry. More precisely, 
Rx is given by 
Ry = = [Stationary value of J ¢] (A4-5) 
It follows that if, by some means, we have obtained a fairly good approxi- 
mation to Q, we may obtain a better approximation to Rg by computing 
J, and using the formula 


Re = 1k(2r)"JS 


When we use the first approximation exp(—7kv) for Q to compute J¢ it 
turns out that the above formula gives the third approximation, Rf’, to 
the reflection coefficient. 

The magnetic corner may be treated in much the same way. The 
integral equation (5-6) for P(v, 8) becomes, in the notation of this appendix, 


P(ro) = € “sin 0) + K (2r) * | sronatn ,r) dS (A4-6) 
in which the vin dS = dvd@ is integrated from — ~ to + * and @ from 0 tox, 


* See Courant and Hilbert, Methoden der Mathematischen Physik, Julius Springer, 
Berlin (1931), page 176, where a similar problem is treated. 
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as before, and G(r , r) now denotes the Green’s function (5-5). We define 
Ju by 


Jun = [ rove — 2e*" sin Old 
(A4-7) 


—«2(2n)-! I] g(r) P(r)g(ro) P(ro)G(ro , 7) dSo a. 


Jy is stationary with respect to small variations in P(r) when P(r) 
satisfies the integral equation (A4-6). Furthermore, from the integral 
(5-7) for Ru , 


Ry = ix*(rc) [Stationary value of Jy] (A4-8) 


which may be used in the same way as equation (A4-5) for Re. 

J. Schwinger has used variational methods with considerable success to 
deal with obstacles in wave guides.* However, his variational equations 
differ somewhat from those given here. Some light on the relation between 
Schwinger’s equations and the present one may be obtained by returning 
to the simple algebraic equations (A4-1) and (A4-2). A rough analogue 
of the expression required to be stationary in Schwinger’s theory is 


(a1 + 24 2X 1X2 a dooXs) /(byx, + box2)* (A4-9) 


The essential point here is that the stationary value of the expression 
corresponding to (A4-9) gives the value of an impedance or combination 
of impedances appearing in some equivalent circuit. Expression (A4-9) 
may be obtained by expressing J, defined by (A4-2), as a function of x; 
and y = %2/x,. J is still to be made stationary but now it is a function of 
x, and y. Solving 0J/dx,; = 0 for x; and setting this value of x; in J gives 
the following function of y 


—(b1 + bey)? (air + 2airey + doy"), 


which is the stationary value of J with respect to variations in x, when y is 
held constant. This function is still required to be stationary with respect 
to y. The same is true of its reciprocal which becomes (A4-9) when both 
numerator and denominator are multiplied by x} and the definition of y 
used. When (A4-1) is replaced by a larger number of equations similar 
considerations lead to a generalized form of (A4-9). The expression required 
to be stationary by Schwinger is obtained when the sums in the general- 
ized form are replaced by integrals. 

* An account of the method together with applications is given in ‘‘Notes on Lectures 


by Julian Schwinger: Discontinuities in Waveguides’ by David S. Saxon. An account 
is also given by John W. Miles." 
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A Set of Second-Order Differential Equations Associated with 
Reflections in Rectangular Wave Guides—Application 
to Guide Connected to Horn* 


By S. O. RICE 


In dealing with corners and similar irregularities in rectangular wave guides 
it is sometimes helpful to transform the system, conformally, into a straight 
guide. Propagation in the straight guide may then be studied by an integral 
equation method, as is done in a companion paper, or by a more general method 
based upon a certain set of ordinary differential equations. Here the second 
method is developed and applied to determine the reflection produced at the junc- 
tion of a straight guide and a sectoral horn—a problem the first method is unable 
to handle. The WAB approximation for a single second-order differential 
equation is extended to a set of equations and approximate expressions for the 
reflection coefficient are derived. 


N A companion paper! the disturbance produced by a corner in a rec- 
tangular wave guide is examined by transforming the system, con- 
formally, into a straight guide. Although the medium in the straight guide 
is no longer uniform, an integral equation may be set up and approximate 
solutions obtained. 

In that paper the wave guide is assumed to have the same cross-section 
at +x asat —x. When this is not so, a conformal transformation may 
still be used to transform the system into a straight guide provided one 
dimension of the original cross-section is constant. However, now some 
advantage appears to be gained by replacing the integral equation by a set 
of differential equations. Since two cases appear, corresponding to £ and / 
corners, there are two sets of equations to be considered. 

These two sets of equations are studied in the present paper. After their 
derivation in Sections 1 and 2 several remarks are made in Section 3 con- 
cerning their solution, special emphasis being laid on the problem of deter- 
mining the reflection coefficient. In the remainder of the paper the general 
theory is applied to a system formed by joining a rectangular wave guide 
to a horn (with plane sides) flared in one direction. The reflection coeffi- 
cients for sectoral horns flared in the planes of the electric and magnetic 
intensity, respectively, are given approximately by equations (6-1) and (7-1). 
These approximations assume the angle of flare to be small so that, as it 
turns out, only the first equations of the respective sets need be considered. 

As was mentioned in the companion paper, Robert Piloty has recently 
made use of conformal transformations in wave guide problems. In his 

* Presented at the Second Symposium on Applied Mathematics, Cambridge, Mass., 


July 29, 1948. 
‘See list of references at end of paper. 
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method the propagation function g(v, @) is derived graphically from the 
geometry of the wave guide irregularities and the result used in one or the 
other of two sets of differential equations which are equivalent to those 
derived below. Piloty’s work is scheduled to appear soon in the Zeifschrift 
fiir angewandle Physik under the title ““Ausbreitung el.-magn. Wellen in 
inhomogenen Rechteckrohren.” 


1. Differential Equations when Electric Vector is in (x, y) Plane 


The partial differential equation to be solved is, from equation (2-3) of 
the companion paper’, 


vO . #0 J 
wi + a + [1 + g(r, |KO = 0 (1-1) 
Ov 0g 
where 
) 
ou = Qat@d = Oand@d=-Fr 
1+ g(v,0) = 1+ D> ancosnd = | f'(v + 6) Pr /P (1-2) 
n=0 


k = [(2b/d)” — (b/a)’} , Ao = free space wavelength 


In (1-2), 3 = x« + iy = f(v + 76) is the transformation which carries the 
wave guide system in the (x, y) plane into the straight guide of width @ = r 
in the (7, 6) plane. For the sake of simplicity we shall always assume that 
far to the left the system becomes a straight wave guide of dimensions 
a, 6 (6 < a) such that only the dominant mode is propagated without 
attenuation. This insures that the a,’s (which are functions of v) will 
approach zero as v—> —x. The dimension (of our system) normal to the 
(x, y) plane is a throughout. 

Since the normal derivative of Q vanishes on the walls at 6 = Oand 6@= 7 
we assume 

Q = Fy + F, cos 0+ F2 cos 20+ ---, (1-3) 


where F, , F2, +--+ are functions of v, and substitute it together with the 
Fourier series (1-2) for 1 + g(v, @) in (1-1). 

The equations obtained by setting the coefficients of the resulting cosine 
series to zero are 


” an k = . ‘ 
Io + (1 + ao)k Fo + D) 2. ark, = 0) (1-4) 
oa n=1 
Fr + [1 + ao + dom/2)R? — m Fn + amk’ Fo (1-5) 


RB a) es 
+ 9) jy Rien + Ontm) Fn = (0) 
2 n=l 
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where m = 1, 2, 3,---, F,. = @F,,/dv®, and the prime on indicates that 
the term x = mis to be omitted. In grouping the terms we have assumed 
that Fy is the major part of Q. 

The principal problem is to solve equations (1-4) and (1-5) when the 
fundamental mode F4 is of the form 


Fo _ e~ th v fh Rye” v 1—> —© 


(1-6) 

Fo = T,z(v), v— +o 
in which Rex is a constant and 7T,(v) represents a wave traveling towards 
v= x, Atv= +x F,, Fy.,--- have the form of waves traveling (or 


being attenuated) away from the region around v = 0. As before, we shall 
be mainly interested in determining the reflection coefficient R. 

It is assumed that only the dominant mode is propagated without attenu- 
ation in the straight wave guide far to the left and hence Ff, , F2,--> all 
become zero as v—> — ®. 

2. Differential Equations when Magnetic Vector is in (x, y) Plane 
The partial differential equation is now given by equation (5-1) of the 
companion paper! 
er . eFP sai 
a2? + 36? + [1 + g(r, O)|x P = 0 (2-1) 
2 2 


where the dimension of the system normal to the (x, y) plane is now 8, a is 
the dimension (in the (x, y) plane) of the straight guide at the far left and 
P=0Oat@0=Oand0=fr 

2) 
1+ g(v, 0) = 1+ : dy, cos nO (2-2) 

n=l 

x = 2a/Xy , Xo = free space wavelength 
c= (x — 1)!” 


Since P = Oat 6 = Oand @ 


II 


wT We assume 


> F,, sin n6 (2-3) 


n=1 


P 


II 


where the F’s are functions of v to be determined by the equations 


Fy + (1 + do — 2/2) — WF + 5 20 (ana — angs)Fn = 0 (2-4) 


Pe + Ot 4 ~ etd) ~~ 0, + 5 er 


2 © 

K . 

n 2 iy (a\m n| — Amin) Fa = 0 
2 n=2 


4 
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in which m = 2, 3, 4,--- and the primes on F, and >> have the same 
significance as in (1-4) and (1-5). 
The principal problem here is to solve equations (2-4) and (2-5) simul- 
taneously subject to 
PF, = eter + Rue’, v— —o 
(2-6) 
Fy = Tn (v), ica ++ nw 
which again corresponds to a unit wave in the dominant mode incident from 
the left. 7,(v) and the remaining F’s correspond to outward traveling 
waves as before. /2, F3, --- all approach zero asv— — ~. 


3. Remarks on Solving the Equations of Sections 1 and 2 for the Reflection 
Coefficient 
Suppose that we have a system in which the wave propagation is governed 
by the single differential equation 


“= hy a (3-1) 


where h = h(v) isa positive imaginary function of v, twice differentiable and 
such that h — ic, c being a constant;as v—» —*. We desire the solution 
of (3-1) which, together with its first derivative, is continuous everywhere 
and at + satisfies the conditions 


y _ enter a Re'’, v—~ —@ (3-2) 
y + (ht h’/(2h))y > 0, v > @ (3-3) 


The constant R (the reflection coefficient) is to be determined. Condition 
(3-3), in which the primes denote differentiation with respect to v} is sug- 
gested by the fact that we want y to represent a wave traveling in the positive 
v direction (the factor exp (tw/) is suppressed). In writing (3-3) we have 
assumed that / is such that for large values of v the two solutions of (3-1) 
are asymptotically proportional to* 


y=h . (3-4) 
& = ¢(v) = icv + | (h — ic) do. (3-5) 


Physical considerations suggest that solutions satisfying (3-2) and (3-3) 
exist in most cases of practical importance. However, if the function / is 
picked arbitrarily the corresponding solutions may be incapable of satisfying 

* S. A. Schelkunoff? mentions that this approximation, sometimes designated by 
“WKB’’, goes back to Liouville. The ideas we shall use are quite similar to those in 
Schelkunoff’s paper. 
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the conditions. For example, if 4 = tc/(1 + exp v) then h — ic exp (—?v) 
as v — 2%, and the solutions of (3-1) behave like Bessel functions of order 
zero and argument c exp (-v). It may be verified that these solutions do 
not satisfy (3-3). Again, condition (3-3) may be satisfied without y having 
much :esemblance to an outgoing wave atv = 2%. Thus if # — ta/v as 
v — x, vinc.eases like v" whe en’ — n — a’ = 0. When0 <a <1/2 both 
values of ” lie between 0 and 1, and both solutions satisfy (3-3). Despite 
these sho tcomings it still seems best to etain (3-3) to specify the behavio- 
of yatv = a. 

It should be mentioned that P. S. Epstein*® has obtained the reflected 
wave by transforming the hypergeometric differential equation into the 
form (3-1). This method has been extended by K. Rawer’ who gives a 1G 
number of references in which the approximation (3-4) is used to study 








° ° . . . . . ” a 
propagation in a medium having a variable dielectric ‘‘constant’. An 1,1 
interesting paper on the general subject of reflection in non-uniform trans- i 
mission lines has been written by L. R. Walker and N. Wax°. ; 
1. When most of the reflection occurs in a short interval, say near v = 0’ "i 
R may be obtained by numerical integration of (3-1). One method is to 4 
start at 7 = 0 with the initial conditions y = 1, vy’ = 0 and work outwards tg 
in both directions. Let V,(v) denote this solution and ¥,(v) the solution 
obtained by starting with y = 0, y’ = 1. The general solution is 
y = Ci¥a(v) + CoV (2). (3-6) 
C; and C2 are to be determined by the conditions 
y = (constant) ke , wv > r% (3-7) 7 
y = (ic/h)' "fet + Re] , v<y (3-8) 
where 7 and v are large negative and positive values, respectively, of 2. iq 
These conditions lead to equations for C; , C2, R: I 
[y+ Oy], ve = Q 
ly’ — Oy + (ich)! ? e-*],.., = 0 (3-9) ‘ 
[y’ + Oty — 2(ich)' *Ret],.., = 0 
in which € is given by (3-5) and 
6 =h+th'/(2h). (3-10) 
The required value of R is obtained by letting 1. > —*, v2— » in the } 


4 
expressions, which follow from (3-9), .' 
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y = Co/C, = —[(Ya + OV a)/(V¥s + OY s)] ome 
P= [y'/ylomv, = (Va + VV s)/(Va + 1¥ Jomo (3-11) 


‘1 
R= ((@* + 1)/(@ — P)J.-», exp | ~2ir - 2 | (hk — ic) | 
-00 
where the arguments of Y,(v) and Y,(v) have been omitted for brevity. 

If 4 should change from a positive imaginary quantity to a positive real 
quantity in (2 , v2) and remain greater than some fixed positive number for 
v > vgit may be shown that | R | = 1 (y and Tare real and Im 6* = Im @, 
Real 6+ = —Real @ at v = 2). This complete reflection is to be expected 
from physical consideration. 


2. An exact expression for the reflection coefficient which holds when / 
ee satisfies the conditions following (3-1) (in particular it must not pass through 
zero anywhere in — x <v< &) is 





‘i R = Mic)? [. e * y(v) = h? dv (3-12) 


where & is given by (3-5). Before this integral for R may be evaluated 
y(v), and hence R itself, must be known. Nevertheless, when R is small a 
useful approximation may be obtained by using the IVAB approximation 





y(v) = (ic/h)' e* (3-13) 
Thus 
a i ‘a3 P i 
‘ R=j[ « h de" di : 
? _— (3-14) 
1 (” -2i 5. fdKY 1-320 K 
= ie ies -) — -K lv 
eae E (x) 4 | 
i in which A = —/’. 


The expression (3-12) for R is obtained by letting 7%) — — ~ in the integral 
equation 
‘ef a +o 2 
4 a, } —to </ \723 d } 

y(to) = (ic/ho)'e *° — Galto, t)y(v)h’ —. h* do, 
‘ Ba dv? 
‘ 1 ui |e a” v < to > sa 
Galto,t1) = —3ho'h*; , , (3-15) 

je * q > to 
evo 


&—i= i h dv, ho = h(r), &) = &(v9). 


v 


Ga(v , 2) is the approximate Green’s function suggested by (3-13). The 
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integral equation may be obtained from the differential equation (3-1) and 
the boundary conditions (3-2) and (3-3) by the one-dimensional analogue 
of the method used in Section 3 of the companion paper'. If we multiply 
both sides of 


2 
dy 
dv? 


—i'y = s(v) (3-16) 


(where s(v) has been added for generality) by Ga(vo , v), integrate twice by 
parts over the intervals (7: , 1% — €), (v% + €, v2) with e > O and 1 < 
Y < v2, and finally let ¢ — 0 we obtain 


y(t) = [ : G.(t0, v) | s) —_ y(o)i} a h : dv 
7 2 


(3-17) 
+ Galt , v1) Ly" ~— Vem, _— Gal to ’ v2) [y’ + g° _— . 


Equation (3-15) follows when we put s(z) = 0 and let 1, > —%, 12> , 
It will be recognized that (3-17) and (3-15) are closely related to integral 
equations occurring in the work of R. E. Langer® and E. C. Titchmarsh’. 

When ih has, for example, one or more simple zeros in —~* < vy < @ 
the integral in (3-15) contains a factor which becomes infinite and the 
integral equation fails. However, we shall not concern ourselves with this 
case beyond remarking that it involves results obtained by H. Jeffreys”, 
Langer’, Furry! and others. 


3. So far we have been considering the solution of only one equation 
whereas we really require the solution of a set of equations. If it is apparent 
that most of the disturbance is given by the first equation of the set it may 
be possible to proceed by successive approximations, each of the remaining 
equations being of the form (3-16) with s(v) determined by the solution of 
the first equation. 

Another method of dealing with a system of V equations is that of numeri- 
cal integration. As a contribution towards obtaining the boundary condi- 
tions at large positive and negative values of v we shall state a generalized 
form of the WKB solution. Although this solution is related to the general 
results obtained by Birkhoff!*, Langer’, and Newell!® concerning the asymp- 
totic forms assumed by the solutions of a system of ordinary linear differen- 
tial equations of the first order, it is worth mentioning explicitly. 

Let the mth equation of the set be 

N 
jm = 2, AmaYn, m=1,2,-°-,N (3-18) 


n=1 


where the A,,,’s are relatively slowly varying functions of v (see equations 
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(3-22) for a more precise statement of the assumptions) and the dots denote 
differentiation with respect to v. We shall reserve primes to denote trans- 
position of matrices. It is supposed that Amn = Anm (equations(2-4) 
plus (2-5) satisfy this condition and (1-4) plus (1-5) may be made to do so 
by setting PF, = 2'/2F 9). 

The solution of (3-18) is approximately 


N 
Ym = 2, Smeleede + & df] (3-19) 


f=1 


where the dz are the 2N constants of integration and 


N 
7 Sint = z. AaaSut 


n=l 


e: 2, Sut = 1 (3-20) 


v 
& =| ge dy 
v3¢ 


serve to determine gy , & , and Sn¢ (the iast to within a plus or minus sign). 
We assume the V roots ¢}, ¢2, «+ yy of the determinantal equation arising 
from the first of equations (3-20) to be unequal, and denote by ¢¢ that square 
root of yg? which has a positive real part or, if the real part be zero, which has 
a positive imaginary part. 3, is any convenient constant. 

The approximation (3-19) may be obtained by setting the assumed form 


in (3-18). The result is a set of NV equations of which the mth is 
fim £ Dime E me + Ee = 2) Amn Bn (3-21) 
We also assume 
Lol Kl e715 | Fm] K| me | K| gme? | (3-22) 
Sm = Smo + Smi + Sm2 + °°: 
where gm, and its first two derivatives satisfy inequalities of the type 


| gmo | >> | gmi| >>| gme] oe 


The first and second order terms in (3-21) give, respectively, 


Smo ¢ apes bs Aun Sno _ 0 
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Lm? ae 7 Aan far — +28mo09 + Lmo¢. (3-24) a 


The vanishing of the determinant of the coefficients of the gno’s in (3-23) 







: are B 2 te 
determines .V values of ¢’, the (th being g7. Once ¢? is selected the gno’s : 
are determined to within a common multiplying factor (which may depend 4 
upon v). This factor is then fixed to within a multiplying constant by the i” 






necessary and sufficient condition that (3-24) be consistent", namely, 


z Nm (28mo¥ + £moe) _ 0 (3-25) 


m 










where /i,, is any solution of the transposed system 


Nim g = Du Aum hy, = 0. ‘beg 






Because Ann = Am, we may take hf, to be gmo. Equation (3-25) may 
then be integrated and leads to the second of equations (3-20) when we set 
the constant of integration equal to unity and identify ¢ and gmo with ¢¢ 
and S,,¢, respectively. The first equation in (3-20) follows directly from j 
(3-23). i 

Since equations (3-20) do not completely satisfy (3-24) (gm remains to be 
determined) our WAS solution for a set of equations is not, in a sense, as ; 
good an approximation as it is for a single equation. Nevertheless it still NW 
represents, just as in the case .V = 1, the leading part of the asymptotic 
form approached as the 4,,,,’s vary more and more slowly with 2. 

In matrix form, the WAB approximation to the solution of 


























S rt 


y = Sed + Se *d (3-27) 





where y and d~ are column matrices, .4 and S square matrices, and exp ; 
(+2) a diagonal matrix having exp (+£,) as the (th term in its principal | 
diagonal. The element in the mth row and (th column of S is S,,.¢ whence, 
from (3-20), 








(3-28) 











== 4 





where the primes denote transposition of elements, J is the unit matrix of 
order V and @ is the diagonal matrix having ¢, as the fth term in its principal 4 
diagonal. That the non-diagonal terms of S’S are zero follows from the iid 
first of equations (3-20) and from gy, ¥ gy if f # k. 
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Because of the second of equations (3-20), at least one of the S,.:’s becomes 
infinite as g- passes through zero. Hence we expect the approximate 
solutions to be valid only over intervals in which none of the ge’s become 
zero. When we hold fast to a given solution y; , y2,--- yn of (3-18) as v 
passes through a value which makes one of the ¢gy’s zero, we expect the set 
of constants dz to be replaced by a new set (Stokes phenomenon). What 
is the relation between the new set and the old set? It may possibly be 
much the same as for the case V = 1 (see Jeffreys’’, Langer’ and Furry"). 


4. When the matrix A is such that the approximation (3-27) remains 
valid over the entire interval — x < v < © (unfortunately this restriction 
prevents us from applying the following results to the horn of Section 4) 
the matrix analogue of (3-15) is the integral equation 


% 


y(t) = Soe = f* + 450 [ e = 1S’ + 26S’ + OS'ly(v) dv (3-29) 


— OO 


+ 4S / c= =[S’ — 26S’ — &S"ly(v) do 


9 


where the subscript zero on S and & indicates that they are to be evaluated 


at v = v%. The column matrix y(v) giving the solution of (3-29) is that 
solution of j = Ay which satisfies the conditions 
y= Se*ft+ Séf, wor—e (3-30) 
y = Se *gt, v—> © (3-31) 


where the column matrix f+ (corresponding to the incident wave) is given 
and f~ (corresponding to reflected wave) and gt (corresponding to the trans- 
mitted wave) are to be determined. The elements of f*, f-, and g* are 
constants. It is further supposed that S satisfies the conditions 
S'S — S'S = 0, v—> +o 

which are certainly met if the elements of A approach constants at +. 
In the wave guide problem we assume ¢y to approach the limit 6; as v 
approaches —«. For this case it is convenient to define the fth element 
in the diagonal matrix = as 


v 
f= bert [| (ve — be de 
L— 00 
In any event we have 


v 
=e —-eb= | ® dv 
vo 
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Although the WKB approximation has the same form as (3-30) in the 
region where 1 is finite, we regard (3-30) and (3-31) as being the exact limiting 
forms of y. Hence, gt may differ from ft. 

Letting 1% — — © in (3-29) and comparing the result with (3-30) gives the 
exact result 


f=} : e =[S’ — 26S’ — OS']y(v) dv (3-32) 


which leads to an approximation for the reflected wave when y(v) is known 
approximately. 
The integral equation (3-29) may be obtained by premultiplying both 

sides of j = Ay by the transpose of the approximate Green’s matrix 

(—3Se* = So , v < 0% 

G,(% , v7) = pa 

\—$Se°""" So, v> %. 
and integrating by parts twice. It is seen that each column of Ga(v , 2) 
is an approximate solution of j = Ay, in which the columnar constants of 
integration are the columns of Sj , and represents a wave traveling away 
from v in both directions. Ga (vp , v) is continuous at » = v% and 


Ox - : . 
|2 Galv m | a 2 Gav . | = So Dy Sb i 
ov v=ro+0 Ov v=ry—0 


Thus the mth column of Ga(v , v) gives the approximate values of y,(v), 
yo(v), «++ , ¥n(v), subject to the conditions that all these and all of their first 
derivatives are continuous at v9 = % except ¥,(v) which has the jump 
Jn(te + 0) — Jn(vo — 0) = 1. 

The presence of 


26S’ + oS’ = 5’ — os’SS! 
= (S'S — S’S)S™ 


in (3-29) and (3-32) makes the N variable case somewhat different from the 
case N = 1. 


5. When Zp, and Y», are slowly varying functions of v the approximate 
solution of the transmission line equations 


dVm _ > Pa 
dv - smn n 
: (3-32) 
dS _ = VV 
dv — im mn n 
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where Zmy = Znm and Vinn = YVamis, as in (3-19), 


N 
Vn = >, Sme le de + € © dé] 
f=] 


N 
Im = 2, Tut leeds — &* di]. (3-33) 
f=1 
Here & is the integral of g as given by (3-20), and gy is determined by 
setting the determinant of the matrix ¢’J — ZY to zero. Wheng;, is known, 
Sme and T,¢ are determined (to within a plus or minus sign which may be 
absorbed by the constants d7 of integration) by the relations 


N 
et Saat =. oe > Lan Tne 
n=1 
N 
et j = = Y an Smt (3-34) 
n=1 


N 
be SmeTme = 1 
n=1 
The last condition, which arises from the condition that the equations for 
the second-order terms be consistent, may be regarded as a generalization 
of Slater’s® result for the case V = 1. 


4. Transformation for Wave Guide Plus Horn 


The system to which we shall apply some of the preceding equations con- 
sists of a straight wave guide starting at « = — and running to x = 0 
where it is connected to a sectoral horn. The horn is flared in the (x, y) 
plane only. The dimension of the system normal to the (x, y) plane is 
constant and equal to a or 6 according to whether the electric or magnetic 
vector is in the plane of the horn. 

One might expect that the field in this system may also (in addition to 
our method) be determined by an alternating procedure of the type described 
by Poritsky and Blewett!® using the equations obtained by Barrow and 
Chu" for transmission in the horn. However, we shall not investigate this 
possibility as we are primarily interested in using the system as an example 
to which we may apply the foregoing equations. 

If the total angle of the horn is 2a7, and if the sides of the straight guide 
are at y = Oand y = 3, (assuming the electric vector to be in the plane of the 
horn), the equation of the lower side, i.e., the continuation of the side y = 0, 
of the horn is y = —.x tanaz and that of the upper side is y = 6 + x tan az. 
If zs = x + iy and w = v + 76 then the Schwarz-Christoffel transformation 
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z = f(w) which carries the guide plus horn in the z plane into the straight 
guide with walls at 6 = 0, @ = m in the w plane may be obtained from 


“ = (1 — &")"b/e (4-1) 
dw 


This gives, upon setting 


dz 


F = [fe + 70) ? = [1 — 2” cos 20 + el ib /r , 
dw’ 


the relation (4-2) 
o» 4v 
1+ g(v, 0) = [1 — 2c" cos 26+ e'|* 
from which the a,,’s may be obtained in accordance with (1-2). 
5. Expressions for the a,,’s for Horn 


The Fourier coefficients of 1 + g(v, 8) appearing in (1-2) and (2-2) are 
the same. It may be shown from (4-2) that 


(e"F(—a, —a;l;e = 5: 2S 
1+a= 411+ 2)/Mite , .=0 (5-1) 
{ F(—a, Os ile e’’) , v<O 
and 
(2e*"*"(—a),F(—a,r —aj;r+1;e)/r!, v>O0 


fav 


(—a),(1 + do)ro/(1 +a) , v=0 (5-2) 


2 
Qe = {2 
2e (—a),F(—a,r—a;rt+ 1; e')/r!, v<0 


where the F’s denote hypergeometric functions, r = 1, 2, --- and we have 
used the notation 


(B)o = 1, (8), = BiB + 1)---(B+r-—1) (5-3) 


When x is odd, a, = 0 because of symmetry about @ = 7/2. The expres- 
sions for v > 0 in (5-1) and (5-2) may be verified by expanding the two 
factors in 


ed 216 


1+ g(x,4) = eI —e" "1 — 7)" 


by the binominal theorem and picking out the terms containing e”"”.. When 
v < 0 we use the relation 1 + g(v, 0) = e'""[1 + g(—2, 6)], and when v = 0 
we may sum the hypergeometric series. 

Differentiation of (5-1) and (5-2) leads to 

dae*” F(—a,1—a;l;e “ ; v>0 

2a(1 + do)v—o , v= 0 (5-4) 
dae” FL — a, + aj;2;e", 1 <0 


d | 
ae (1 ot ao) — 








———— 
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tea "ta Pr * Ma. a; 1;¢°"), ?>O0 


? 
- (1 + a) = 2 4p 4v\2a—1 po 4v 
a loa e"(1 — e")" Fla, a; 1; e”), 1 <0 


(5-5) 
where in obtaining (5-5) use was made of Euler’s transformation 
F(a, b;¢c;x) = (1 — x)@*>F(c — a,c — b3 0; x) 
It is seen that d(1 + do)/dv is continuous at v = 0 but the second deriva- 
e “okey 2 1 
tive becomes infinite as 7“ 
When 1 + ap and a2 are expressed as the customary integrals defining the 


Fourier coefficients it is seen that one of the coefficients occurring in equation 
(2-4) for F; is given by 


II 


> 7? ss as 
1 + do — a2/2 - i] (1 — 2” cos 20 + e*”)* sin’ 6 dé 
® 0 (5-6) 


(e” + 1)°* F(—a, 4; 2; sech’ v) 


II 


Atv = 0,1 + ad) — a/2 and its first and second derivatives are continu- 
ous, their values being 
I'(2 + 2a) 2aI'(2 + 2a) 
ritar(2+a)’ Fitaelr(2+a)’ 
da(2e? + 2a + 1)r(1 + 2a) 
rl + a)P(2 + a) : 


(5-7) 


respectively. These may be obtained by differentiating the integral in 
(5-6) and setting » = 0. 
A second expression for 1 + ay — d2/2 follows from (5-1) and (5-2): 


fav . op» 
e“"|F(—a, —a; 1; e*") + ae?’F(—a, 1 — a; 2; vas) F 


v>0 
1+ ap — ao/2 =< 


| F(—a, —a; 1; e*”) + ae’’"F(—a, 1 — a; 2; e*”) 


v< 0 
(5-8) 


6. Approximation to Reflection Coefficient of Horn, Electric Vector in (x, y) 
Plane 


When the flare angle 2am of the horn is very small the reflection coefficient 
may be shown to be 


i - 
R; = i + O(a?) (6-1) 








i50 BELL SYSTEM TECHNICAL JOURNAL 


where O(a’) denotes correction terms of the order a. This result is based 
upon the fact that when terms of order a” are neglected the set of dif- 
ferential equations (1-4) and (1-5) reduce to the single equation 


Fo + (1+ av)k*Fo = 0 (6-2) 
where, from (5-1, 4, 5), 
v>0O v<0 
1 + ap eta 1 
d 
= (1+ do) 4ae*” 0 


9 


7? 
as (1+ do) 16a%e*"(1 — e#”)?" = 16a%e**(1 — et?) ?e-1 
dv" 


The reflection coefficient (6-1) is the one corresponding to the differential 
equation (6-2) and may be computed by setting 


(1+ a) = -—PW=K (6-3) 


in the integrals (3-14). 

The expression (6-1) for Re may be obtained quickly (but the procedure 
is not trustworthy) by assuming that the principal contribution to the first 
integral in (3-14) comes from the region close to v = 0, say in -~e <0 < €, 
where the second derivative of 4~'” is infinite but integrable. When the 
integration is performed approximately by replacing the second derivative 
by the first, (3-14) gives 


1f,4d is 
Re au = h J (6-4) 


€ 


1 jd ,/° ta 
2ik EX + ao) | . Qk 


where ¢€ is assumed to be so small that 1 + ag is effectively unity and 
d(1 + ao)/dv changes from 0 at — € to ta at + €. 

A more careful investigation based on the second integral in (3-14) also 
leads to the value (6-1) for Ry. It further suggests that possibly most of 
the correction term, denoted by O(@’) in (6-1), is given by 


Ut 


lt 


2 ce) Y 


a 2t-2av 1 eé . 
nie patie 2.  £Sife) — 2/2 iCilx 5-5 
vik dy e ( ree + r [Si(x) — 2/2 + iCi(x)] (6-5) 


with x = k/a and 2 = ixlexp (2av) — 1]. Si(x) and Ci(x) denote the 
integral sine and cosine functions. Incidentally, the rather curious result 
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2 © 1 a 
n=1 m=1 mn(m + n) oi > ‘2 

turned up in the investigation of the orders of magnitude of the various 

terms. 


7. Approximation to Reflection Coefficient of Horn, Magnetic Vector in (x, y) 
Plane 


The work of this section is quite similar to that in Section 6 except that 
here we enter into more of the details. We shall show that when a is small 
the reflection coefficient appearing in equation (2-6) is 


Ra = ~~ + Oa). (7-1) 


2c8 
From (2-4) the analogue of the differential equation (3-1) is 
FY + [1 + ao — a2/2) — 1]Fi = 0 (7-2) 
and the K appearing in the second of equations (3-14) is now 


= —/? = (1 + ay — a2/2) — 1 (7-3) 


The largest terms in the expression (5-8) for 1 + ao — a@2/2 yield, to within 
terms of O(@), 


K = w(e*" + ae’) — 1, »v>O0 

K = (4a e@° — Jae’) (7-4) 
K = «(160% + dae-*”) 

K = #®(1 + ae’*) —-1= C+ rae” , 1< 0 

K = 2axe” (7-5) 
K = 4axe?” 


where the dots denote differentiation with respect to v and c? = x? — 1. 
We have retained the @? in K as given by (7-4) because at this stage we do not 
know whether it may be neglected or not. 

When v < 0, the definition (3-5) of € and (7-5) yield 


g = ico + if (a — c)dv 
(7-6) 


I 


icv + ic | (1 + °c? ae”)' — 1| dv = icv + O(a) 
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and we have 


= O(a’) 
which may be neglected. The other integral suggested by (3-14) is 
a) - 0 
2 p---3/2 je 2icv —3 2 20 
| e °K Adv= | € ¢ take dv 
oo — © 
= Jax’ c°/(1 — ic) 


When v > 0, 


ee (e+ ae?) — 1}! ae 
0 


> 


] [ (x _, 1)! dv + O(a), 
0 


i 
= = [y — tan’ x —c + tan cl] + O(a), 
£Q 


2 4av } 2\—1 
(x e*" — 1)’, 2a dv = x1 +x) dx 
In the integrals containing exp (— 2v) as a factor, & may be taken to be icv 
since the integrand becomes negligibly small by the time icv differs signifi- 
cantly from (7-9). We have 


l 


a 


@ 
2t ,-—5/2 yf-2 2t7 2 sav . 4.3 dav 
as 4 KK’ di -{ e “(xe —1) ‘xk a(4e” — 
0 


@ 

2 -—5 4 24, 8av 
[ ex «x al16e" dv 
/0 


a 
ss ~“Leettan tells 9 
ail sa | e i{r—tan~!r tan ake {: 4 fe x °) dx 
e 


where the integrals containing e~?” and e~*” have been neglected since their 
contribution is O(a”). When a becomes exceedingly small the exponential 
term oscillates rapidly and the last line of (7-10) is likewise O(a”). This may 
be verified by integrating by parts, starting with 

exp V dx = iax (1 + x*)d(exp Y), 


Y = —i(x — tan'x) /a 


The last integral which must be considered is 
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ry 


-_ [ etigeh ~17 ate + 4a “ab 
0 


= 16x'a° | e* 
0 


(7-11) 


@ 
Qicv-2v —3 2 
+ | e¢-  ¢ K 4a dv 
0 


oo 
ak teat eae tothe 2 2-8, . 
= Sa | en en ee + 2ak ¢ / (1 + ic) 
c 


= O(a’) + Par? co? /(1 + ic). 


That the integral having « as the variable of integration is O(a”) may be shown 
as in (7-10). 
When we combine our results in accordance with (3-14) we obtain 
1 r —2tr 5 p-—5/2 peo 1 p-—3/2 EP 
Ry =< e “l3sK °K? —1K “" K] dt 


21 J-x 


2 3 
aK ¢ 1 1 9 
-—- ies - | + Oa) 
47 , — 1 1+ tc 
ia/(2c*) + O(a) 
which is (7-1). 

If, instead of discarding (7-10) because it is O(a’), we retain it and the 
corresponding integral in (7-11) (in the hope that they represent most of the 
difference between the approximate value (7-1) for Ry and the true value) 
we obtain the approximation 


ina. [ ge ete eee 4 de: (TAD) 
2) 4 
in which the integral may be evaluated by numerical integration. 

The approximations (6-1) and (7-1) for the reflection coefficients may also 
be obtained from an equation given by N. H. Frank.'* However, care 
must be taken to suitably define the wave guide characteristic impedance 
which appears in his expression. 


8. Speculation on the Reflection Obtained from Horn Flared in Both Directions 


All the work from Section 4 onward applies only to a horn flared in one 
plane. Nevertheless, it is interesting to speculate on how close an estimate 
of the reflection from a three-dimensional horn may be obtained by super- 
posing the two reflection coefficients (6-1) and (7-1). It must be kept in 
mind that the flare angles (the a’s) may be different in the two directions, 
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that & is given by (1-2) and ¢ by (2-2), and finally the difference (not the 
sum) of Rg and Ry must be taken. In (6-1) Rg is the reflection coefficient 
of the component of the magnetic vector normal to the (x, y) plane (which 
is proportional to Q), while in (7-1) Ry is the reflection coefficient of the 
transverse electric vector (which is proportional to P) and there is a difference 
in sign just as in the case of voltage and current reflection coefficients. If 
a > band Xo is the wavelength in free space, the superposition gives the 
following expression for the reflection coefficient of the electric vector: 


R= Ry — Rz 


Pia all : (8-1) 
= 5 [(2a/Xo) — 1] *(@x/[(2a/do)” — 1] — aag/b) 


where 27a, and 2ma,z are the total horn angles in the planes of H and £, 
respectively. Of course this approximation can be expected to hold only 
when ag and ag, are small. 


9. Numerical Calculations—Ry for 60° Horn 


The value of Ry , the reflection coefficient when the magnetic vector lies 
in the plane of the flare, was computed on the assumption that only the 
dominant mode need be considered.* Thus, instead of the system of 
equations (2-4) and (2-5), only their simplified version, namely the single 
second order differential equation (7-2), was used. This equation may be 
written as 

dF, 


—-+ KF, = 0 (9-1) 
dv 


where, according to (5-6), 
K = —f = (1+ a — a,/2) — 1 (9-2) 
1+ ado — do/2 = (e" + 1)*F(—a, 1/2; 2; sech? v). 


The problem was to obtain the Ry appearing in that solution F, of (9-1) 
which satisfies the boundary conditions (2-6). 
No computations for Rg were made. 
In the first method of calculation the integrals in the approximation 
(3-14), namely 
oo 2 
Ry = . | e igi : - K ms dv (9-3) 
21 J-« dv 
‘ 


2& = icv + 2i (K' — c) do, 


i OO 


* I am indebted to Miss M. Darville for carrying out the computations of this section. 
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were evaluated by Simpson’s rule. The second derivative of K~'/* was 
computed from the even order central differences of K'*."* For a= 1/6, 
corresponding to an angle of x/3 between the two sides of the horn, calcula- 
tions at two representative wave lengths led to the table 


No Cc K = 1+¢ Ru (9-4) 
1.549a 8173 1.6680 — .0420 + 7.0724 
1.610a .7376 1.5441 — .0551 + 1.0878 


An idea of the variation of K may be obtained from its values at —@, 
—.6, 0, .6, 1.8, 3.6 which are approximately .67, .76, .98, 1.62, 4.56, 17.4, 
respectively. The range of integration was —3 < v < 4.4. 

The second method of computation used the formulas (3-11) with F; 
playing the role of y. The differential equation (9-1) was integrated by the 
Kutte-Runge method, the interval between successive values of v being 


0.2. For c = .8173 the values obtained were 
vy Ve Y r Ru (9-5) 
—.6 6 —.202 — 1.981 —.142 — 7.794 —.0167 + 7.0658 


—1.2 1.2 —.218 — 11.004 —.049 — 7.696 —.0525 + 1.0754 
—1.8 1.8 —.225 — 1.989 +.086 — 1.716 —.0512 + 7.0753 
—2.4 2.4 ~.220 — 71.000 .136 — 71.842 —.0424 + 1.0722 


In order to gain an idea of the meaning of these values of v it should be 
recalled that w = v+ 76 and the walls of the guide are at @ = 0 and @ = 7. 
An interval of length r = 3.14 --- in the v direction therefore corresponds 
roughly to a distance equal to the width of the guide. The above table 
indicates that, loosely speaking, most of the reflection occurs close to the 
junction of the horn and wave guide. 

The last value of Ry in (9-5) agrees quite well with the value —.0420 + 
1.0724 obtained from the approximate expression (9-3). It appears that the 
method leading to (9-5) is superior to the one based on (9-3) since, in theory, 
it may be made as accurate (insofar as the single equation (9-1) may replace 
the set of equations (2-4, 5)) as desired. Moreover, less actual work seems 
to be required. 

The approximation (7-1) yields, for ¢ = .8173, 

la i(§) eo 
he = aa” gq 
which is considerably in error, as we might expect, sincea = 1/6 is not small. 
However, if we use the approximation (7-13) and evaluate the integral 
by Simpson’s rule we obtain 
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Ry = 1.153 — (.061 + 1.077) 
= —.061 + 7.076 


which is in better agreement with the earlier values of Ry . 

No similar computations have been made to test the corresponding 
approximation for Ry obtained when the correction term (6-5) is added to 
the leading term in (6-1). However, it appears that fora = 1/6 and the 
representative value k = .38, (6-5) is only about one sixth as large as ta/ (2k) 
and hence is relatively unimportant. 
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Abstracts of Technical Articles by Bell System Authors 


Pulse Echo Measurements on Telephone and Television Facilities.| L. G. 
ABRAHAM, A. W. LesBert, J. B. Macoio, and J. T. Scuorr. Pulse echo 
measurements have been used on telephone and television facilities since 
1940 to locate impedance irregularities and control quality in manufacture 
and installation. These sets send a pulse into a line and observe on an 
oscilloscope the echoes returned from irregularities. The shape and width 
of the pulse, the rate at which it is repeated and the pulse magnitude are 
important in determining the accuracy of the results and the requirements 
of the measuring apparatus. The “‘coaxial pulse echo set”’ is used for factory 
and field testing of coaxial cables. The ‘“‘Lookator’’ was developed for use 
on much narrower band systems such as spiral-four field cable and open 
wire lines. 

Television Network Facilities’ L.G. ABRAHAM and H. 1. Romnes. This 
paper describes television network facilities which are needed to connect 
studios and other pickup points to transmitters in the same and in distant 
cities, and discusses their transmission characteristics. Short-haul tele- 
vision circuits may be by microwave radio or over wire circuits. Long- 
haul television connections may be by radio relay or over coaxial systems 
of the type originally developed for carrier telephone circuits. Transmission 
requirements include adequate frequency band, accurate gain and phase 
equalization, and freedom from interference resulting from excessive noise, 
crosstalk, or modulation. Radio and wire systems are under development 
to provide extensive high-quality television networks. 

A Carrier Telephone System for Rural Service.’ J. M. Barstow. The 
M1 carrier telephone system was designed for the purpose of extending tele- 
phone service into areas served by rural power lines, but not served by co- 
existing rural telephone lines. To the local office operator and to a carrier 
subscriber the service provided is the same, so far as procedures involved 
in establishing a connection are concerned, as a voice-frequency line. 

At the office end of the system a telephone wire line extends from the 
office to a point near the power line. Here is located a converter (called 
common terminal) which converts the voice-frequency signal to be trans- 
mitted to the subscriber to an amplitude-modulated double-sideband carrier 
signal. This signal is coupled to the power line through a coupling unit 

1 Trans., A. 1. E. E., vol. 66, 1947 (pp. 541-548). 


2 Transactions, A. 1. E. E., vol. 66, 1947 (pp. 459-464). 
3 Trans. A. 1. E. E., vol. 66, 1947 (pp. 501-507). 
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and high-voltage capacitor. At the subscriber’s location the signal is taken 
off by similar means and led by separate wires to the subscriber premises, 
where it is reconverted to voice frequency by means of a subscriber terminal. 
A signal transmitted from the subscriber to the central office goes through 
similar conversions. 

The usual number of parties per two-way channel may be assigned ac- 
cording to local custom, and divided-code or full-code ringing is provided. 
Equipment is available for five two-way channels over a single-power line 
employing frequencies in the range 150 to 425 kilocycles. A sixth channel 
has been discontinued because of radio interference. 

A description is given of the manner in which the power line should be 
treated in order to reduce reflection effects. The power line treatment does 
not affect its capabilities in regard to power transmission. 

A pplication of Rural Carrier Telephone System’ FE. H. B. BARTELINK, 
L. E. Cook, F. A. Cowan,* and G. R. Messmer. This paper deals with 
the application of a carrier system developed primarily for providing rural 
telephone service over power distribution circuits in areas where this means 
of extending telephone service may be more attractive than other avail- 
able methods. The modifications required in the power circuits to permit 
carrier frequency transmission are described, including the effect of these 
modifications on the operation of the power system. Construction fea- 
tures also are discussed. The use of the rural carrier telephone system over 
open wire telephone pairs is discussed briefly. 

An Improved Cable Carrier System H.S. Biack, F. A. Brooks, A. J. 
WrerandI.G. Witson. A new 12-channel cable carrier system is described 
which is suitable for transcontinental communications. Important fea- 
tures are negative feedback amplifiers of improved design, new arrangements 
for accurate equalization of the cable loss, and automatic thermistor regu- 
lators which continuously control the transmission of each system. 

Joint Use of Pole Lines for Rural Power and Telephone Services® J. W. 
CAMPBELL,* L. W. Hitt, L. M. Moore, and H. J. Scuoiz. The use of 
poles to carry both power and communication circuits is not new, having 
been employed before 1890. There are today more than 6,000,000 poles 
used jointly by power and telephone organizations in the United States. 
The great bulk of these poles are located in urban areas where the voltages of 
the power circuits concerned are generally less than 5,000 volts and the 
span lengths between poles generally do not exceed about 150 feet. 

As power and telephone lines were extended into rural territory, new 


4 Trans. A. I.E. E., vol. 66, 1947 (pp. 511-517). 
5 Trans. A. I. E. E., vol. 66, 1947 (pp. 741-746). 
§ Trans. A. 1. E. E., vol. 66, 1947 (pp. 519-524). 
* Of Bell Tel. Labs. 
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problems were encountered in the application of joint construction because 
of the use of longer spans and higher voltages for the power circuits and the 
increased noise induction in the necessarily longer exposures. However, 
progress in the art through cooperation of the telephone industry with the 
Edison Electric Institute and the Rural Electrification Administration has 
brought about the developments reviewed in this paper which now make 
long span higher voltage rural joint use feasible where conditions are favor- 
able. 

Atomic Energy.’ Kart K. Darrow. (The 1947 Norman Wait Harris 
Lectures at Northwestern University.) This little book, which reproduces 
four lectures substantially as they wer? given, is at once a very readable and a 
very accurate account of enough of the facts of nuclear physics to convey a 
good understanding of the atomic bomb and the possibilities of atomic 
power. The scientific accuracy of the presentation is instanced by the 
author’s apologies for his title; he emphasizes that in reality his subject is 
nuclear energy, but that on the day of Hiroshima somebody wrote of an 
atomic bomb and the misusage spread like a chain reaction. 

The role of electrons, protons and neutrons in atom building is told in a 
simple and entertaining style (but with a degree of ornamentation that may 
disturb some readers), and the discussion of rest mass and the Einstein rela- 
tion between mass and energy is pointed up by well-chosen numerical il- 
lustrations beginning with the lightest composite nuclei. The role of fast- 
particle bombardment in increasing and decreasing the size of nuclei is also 
explained. The reader thus acquires a clear understanding of the basic 
phenomena for which nuclear fission is famed. ‘The text is augmented by 
well-chosen cloud chamber photographs. 

Though the author’s treatment is accurate, his style and marshalling of 
facts are very readable. This is well illustrated by the closing paragraph 
of the third lecture, which follows immediately upon the author’s develop- 
ment of the idea of the chain reaction: 

Here is the climax of my lectures, and here is where you should be frightened; and. 
if I had an orchestral accompaniment, here is where the orchestra would have mounted 
to a tumultuous fortissimo, with the drums rolling and the trumpets blaring and the 
tuba groaning and the strings in a frenzy, and whatever else a Richard Wagner could 
contrive to cause a sense of Gétterdimmerung; for, let there be no doubt of it, this is 
something that could bring on the twilight of civilization. But at this crucial junc- 
ture I have only words to serve me, and all the words are spoiled. We speak of an 
awful headache, a dreadful cold, a frightful bore, and an appalling storm; and now 
when something comes along that is really awful and dreadful and frightful and 
appalling, all these words have been devaluated and have no terrorin them. I have 
to fall back on the saying, of unknown origin and dubious value, that the strongest 
emphasis is understatement. Let then this picture, with its circles and its symbols 


and its numbers, be considered an emphatic understatement of the most terrific thing 
yet known to man. 


eS. Publishedby John Wiley & Sons,Inc., New York, and Chapman & Hall, Ltd., London. 
80 pages. $2.00. 
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The book will be welcomed particularly by those who at one time or 
another have had a general acquaintance with radio-activity, cosmic rays 
and the results of cloud chamber research, but whose vocational activities 
have forced such special knowledge well into the nebulous regions of their 
memories. 

Network Theory Comes of Age.® KR. L. Diretzo_tp. The third decade in 
the growth of modern network theory, the decade of maturity, is considered 
in this review of the advances in network theory evolved over the past ten 
years. New types of networks developed during the war are included. 

Thermistors as Components Open Product Design Horizons. K.P. 
DowELL. These thermally sensitive resistors with high negative tempera- 
ture coefficients have come a long way since they were laboratory curiosities 
and are now available in a wide range of types with diverse and stable char- 
acteristics. You may be able to transfer to your own problems some of the 
unusual design ideas described here. 

Gas Pressure for Telephone Cables... R. C. Giese. Communication 
cables consist of a number of electric conductors insulated from one another 
and encased in a metal sheath. This encasement is subjected to numerous 
hazards, such as those caused by electrolysis, crystallization, various kinds 
of mechanical damage, and lightning burns. Any damage to the sheath 
which will permit water to enter the cable will decrease the effectiveness of 
the insulation material and thus cause an impairment or an interruption to 
the service. The entrance of moisture through small openings in the sheath 
can be materially retarded when the space inside the cable, not occupied by 
the conductors or insulation, is filled with a gas maintained under controlled 
conditions. Nitrogen is the gas usually used for this purpose because it 
is inert and does not combine chemically with the conductors or insulation. 
In addition the use of the gas provides a method of locating openings in the 
sheath by means of a pressure gradient, which is a material aid in cable main- 
tenance. 

Rural Radiotelephone Experiment at Cheyenne Wells, Colo." J. HAROLD 
Moore, Pau K. Sey_er and S. B. Wricur. The first rural party-line 
telephone service by radio installations operating on the subscribers’ prem- 
ises was inaugurated August 20, 1946. This paper describes the equipment 

“used, how it operates, and the results obtained during the preliminary testing 
and the initial period of regular operation. Radio is one of several new 
methods which the Bell System is exploring in its program for extension of 
telephone service in rural areas. It is expected that experience gained in 

8 Electrical Engineering, September 1948 (pp. 895-899). 

» Elec’l. Mfg., August 1948 (pp. 84-91). 


® Transactions, A. I. E. E., vol. 66, 1947 (pp. 471-478). 
" Trans. A. I. E. E., vol. 66, 1947 (pp. 525-528). 
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this experiment will aid in developing a standard rural radiotelephone 
system. 

Effect of Passive Modes in Traveling-Wave Tubes." J. R. Pierce. As 
the beam current in a traveling-wave tube is increased, the local fields due 
to the bunched beam become appreciable compared with the fields propa- 
gating along the circuit. The effect is to reduce gain, to increase the electron 
speed for optimum gain, to introduce a lower limit to the range of electron 
speeds for which gain is obtained, and to change the initial loss. 

New Test Equipment and Testing Methods for Cable Carrier Systems." 
W.H. Tipp, S. Rosen and H. A. WeNk. ‘Three portable test sets developed 
for the improved cable carrier telephone system are described: A high sen- 
sitivity selective transmission measuring set covering 10 to 150 kc, a decade 
oscillator for frequencies from 2 to 79 kc, and a tube test set. 

A New Microwave Television System.’ J. F. Wentz and K. D. Smirn. 
A microwave point-to-point radio system is described which is designed for 
the transmission of television programs. This system is intended to sup- 
plement wire facilities for local distribution of television signals from pickup 
points to studios or from studios to broadcast transmitters and to long dis- 
tance network terminals. The circuits and equipment are described in 
detail. Performance obtained in tests during 1946 is given. 

2 Proc. 1. R. E., August 1948 (pp. 993-997). 


'8 Trans. A. 1. E. E., vol. 66, 1947 (pp. 726-730). 
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